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PREFACE 


DunING the session 1954-55, I was invited by the authorities of the 


-Maharaja Sayajirao University of Baroda to deliver a course of six lectures 


on some selected topics in Design and Analysis of Experiments for the 
benefit of the post-graduate students in Statistics. The invitation was 
repeated during the session 1955-56 also. The present book covers more 
or less the substance of these twelve lectures. I have received numerous 
requests from friends and students for copies of these lectures in the 
recent past and this demand has encouraged me to print these lectures 


in the present form. 


forms a fascinating branch of Statistics and 
though primarily it originated from agricultural experiments, it is finding 
more and more applications in various pther fields. In India Dr. R. C. Bose 
started research in the mathematical ‘onstruction of designs in 1938 and 
he has influenced a large number of workers in India and America. I 
had the privilege of attending some of his lectures in 1943. Most 
Universities offer advance courses in Design and Analysis of Experi- 
ments and it is hoped that this book will be found useful to the students.“ 
I have not included any experimental material for two reasons. Firstly, 
I have not personally conducted experiments in any experimental station. 
Secondly, there are many excellent books dealing with this aspect. While 
teaching the subject, I have felt the need for solving a large number of 
problems which only can give a deeper insight into things. Unfortunately 
such a collection of problems is not easily accessible to the students. 
Kendall in his book, Exercises in Theoretical Statistics (Charles Griffin 
& Co. Ltd., 1954), has removed a long-felt want of students but he has 
also not included exercises on experimental designs as he had ‘to draw 
a line somewhere.’ The present book contains more than 150 problems 
and I believe that a genuine attempt at solving the problems will give the 
students a firmer grasp over the subject. 


Design of Experiments 


No, originality is claimed except perhaps in the presentation of the 
material. I am deeply grateful to numerous friends and students for 


encouraging me to undertake this work. I shall greatly appreciate if any 
reader will bring to my notice mistakes, misprints or any defect in the 


book. 


M. C. CHAKRABARTI 


Bombay 
15 July 1962. 
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CHAPTER I 


THEORY OF LINEAR ESTIMATION 


1. Notations. Throughout this book, letters in thick-faced type will stand 
for matrices. In order to ecconomise space, the single column matrix. 
E 


a 
ds 
d 
will be denoted by (4: 45...4,). Ip will be used to denote the identity 
matrix of order p. The p x q null matrix will be denoted by O,, and when 
no confusion is likely, by O. The p x g ‘matrix all of whose elements are 
unity Will be denoted by Ej, Clearly 


Ej, X Ej — QE, (1.1.1) 


The square matrix of order p having ay, do, ..4 a, along the principal diagonal 
and zero elsewhere will be denoted by diag (a, ay, ...,a,), The letter E will 
stand for mathematical expectation and matrix V will stand for a dispersion 


< 


matrix. 


2. Estimable Parametric Function and Condition for Estimability. Let y;, 
Unt be n independent stochastic variates having a common variance o? 


and expectations given by 


Yı 3; 319 +++ lim 9 
E(y) =E Yo = | 4g dag - Com 0. =A0 (1.2.1) 
Yn Lan Gna +++ lnm 8, 


* 
where the ays are known constants and 0; 0, ... Om are unknown para- 
meters. A linear function i 


b' 0 = 0,0, + ba 0, + ... + bm Om (1.2.2) 
is said to be estimable if there is a linear function 
€ y = C141 + Ca Ya F -a H Cn Yn (1.2.8) 


such that 
E(c'y) —b'8 (1.2.4) 
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irrespective of the values of the parameters. Evidently then, 


c’A =D’ or A’c=b (1.2.5) 
and a necessary and sufficient condition for b' 0 to be estimable is that 
Rank A’ = Rank (A’, b) (1.2.6) 
(1.2.6) is equivalent to 
Rank A’A = Rank (A'A, b) (1.2.7) 


For, if (1.2.6) holds 
Rank AA < Rank (A'A, b) = Rank (A’, b) [ 
= Rank A'A 
Hence (1.2.7). On the other hand, if (1.2.7) holds, we can find x = (a, 25, 
+++) Em} Such that (A'A) x — b. Writing z =A x, we can see that A/z =b 
is solvable and hence (1.2.6). Oftentimes (1.2.7) is easier of application 
than (1.2.6). 

When b'0 is estimable, we can find a unique linear estimator (1.2.8) which 
(i) is unbiased and (ii) possesses the least variance among unbiased cstima- 
tors. This will be called the best unbiased linear estimator of the estimable 
linear parametric function (1.2.2). Since (1.2.2) is estimable we have at least 
one linear function (1.2.3) satisfying (1.2.4) and hence (1.2.5). This leads to 
(1.2.6) and (1.2.7. Now E (c’y —b’0)?=o* c'c. Hence we have to 
minimise c'c subject to (1.2.5). Let A — (24, Ag, ..., Am} denote Lagrangian 
multipliers ; differentiating c'c — 2A’ (A’c — b) with respect to c, we get 

c=AA (1.2.8) 
b=A’AA (1.2.9) 
On account of (1.2.7), we can find A from (1.2.9) and hence c from (1.2.8). 
If A'A is singular (1.2.9) will not give a unique A but no matter which solution 


of (1.2.9) we take c is the same. Let AM and A? be two distinct solutions 
of (1.2.9) and let c(D = AA, c(2 = A 42. Then 
1 


O = AA (40 — A?) 
(c) — cl)? (c) — c) = (AM — AD) A'A (AM — 2) =O 
zoe = c 


A On, , 
Os 3 < Rank (A’, b) 


That the minimum is actually obtained can be demonstrated in the following 
manner. Let us assume for the sake of simplicity that the first r columns of 
A are independent, the rest depending on them and A, is the matrix of these 
r columns and A, the matrix of the remaining columns. Then A = (A, A). 


Let b, = (bi, ba, ..., b,}. A solution of (1.2.9) is 
b 
A , =1 1 
(A5 Ay) (s ] : 


mer, 1 
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Hence (1.2.8) can be written as 


7 Sb. 
, c = (Ay, As) (A4 A1) * [a] = A (AA) by... (1.2.10) 


and c'c = b^ (A,A;) + b, (1.2.11) 
Now for every c satisfying (1.2.5) 

e'c — c' [E— A; (A^, Ay) A4] c + b'i (A,A,)! b, 

— c' [I— A, (AA) A4]' [I — A; (AA) TA] e + DY (AA) bi 
< b'i (AA) bi. 

Hence the minimum is achieved for c given by (1.2.8). 

In the general case where ?,th, ith, ...7,th columns of A are independent, 
the rest depending on them and A, is the matrix of these independent 
columns, b; = fb; , bip, «++» bi, }, c and c'c are given by (1.2.10) and (1.2.11). 
Also it is evident that bi, bi. 2:5 b; can be chosen arbitrarily and any other 
constituent of b, say bj, in order that (1.2.2) may be estimable, must be 


linearly related to b; b s; in the same way as jth column of A is 


is, Us oc 
related to the ith, Mino. ith columns. 
A linear function 
ly =hy + lY + DEI (1.2.12) 
will be said to belong to error if > 
EQy)—0 | (1.2.18) 
irrespective of the values of 0,, Os ..., Om A necessary and sufficient 
condition for this is 
Le Q) (1.2.14) 
If œ Op; +++» %m denote the column vectors of A, by (1.2.8) 


€ = 404 F Ag Oe + «+--+ Am 05, Where A= Du Ao, aas 24] 


is given by (1.2.9). Thus c lies in the vector space generated by 0t, a, ..., m 
which is also called Estimation Space. The vector 1 correspondidg to a linear 
function I’ y belonging to error is by (1.2.14) orthogonal to the estimation 


space. ‘The vector space orthogonal to the estimation space is called the 


Error Space. 
3. Method of Least Squares and Markoff's Theorem. As estimate of ® we 
usually take @ which minimises 

(y — A6)' (y — A8) a 
Differentiating (1.8.1) with respect to 0 and equating the derivative to zero 
e get: the normal equations > 


wt 


A'y — A'A8 (1.8.9) 
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The above system of equations is always solvable for Rank A’ = Rank 
(A’, A’y) and hence Rank (A'A) = Rank (AVA, A'y). The solution is unique 
if A'A is non-singular. Let b'0 be an estimable parametric function, its hest 
unbiased linear estimate is from (1.2.8), (1.8.2) and (1.2.9) 
c'y = A'y —- XA' A8 — b'6 (1.8.8) 
This provides an easy method of obtaining the best unbiased linear estimate 
of an estimable linear function of the parameters. It is also evident that 
(i) E(A'yj— A'A0 (1.8.4) 
(ii) V(A'y) E [A' (y — A0) (y —A0) A] =c2 AVA ^ (etu) 
(wit) c'y — A' (A'y), i.e. the best unbiased linear estimate of an estimable 
linear function of the parameters is: a linear function of the elements of A’y. 
(iv) If in (1.8.2), we can formally write 6= D A’y, where D is a square 
matrix of order m, then in (1.2.9), we can write A = Db and the variance of 


b 8, the best estimate of the estimable linear function b'6, is A'b = c?b'D'b. 
One way of choosing D is to take it equal to 


(&,A;)7! O 
(0) o f 


(w) When A’A is singular, 6 is not unique. But no matter which 


solution of (1.3.2) we take A @ remains the same; just as in (1.2.10), we 
can thus take 


A6 — A, (A^, Aj)71 A^, y (1.8.6) 
We can find a non-singular matrix G of order r such that 
(A’, A,)7! = G’G (1.8.7) 
Hence 
Rank A,(A' Aj)! A', = Rank (A, G' GA’,) = Rank (G A^) 
= Rank A^, =r (1.8.8) 
Rank [I — A, (A^, A) A4] 24 — y (1.8.9) 


Also I — A, (A'5 Aj)? A’, is symmetric and idempotent and hence its non- 
zero eigenvalues are all +1. Hence 


E [y — A8y (y — A6] =E (y — A0) [I — A (A,A)-1 A4] (y — A0) 


= (n —r)o? (1.3.10) 
As 


E (y' y) — nc? +0 A'A0 


and 
yy —(y — A8y (y — A6) 6 A'A0 
E (6 A'A8) «rct 0 A/A0 (1.8.11) 
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(1.8.8) and (1.8.10) constitute Markoff's Theorem, though Markoff’s result is 
slightly less general than what we have proved above. 
he expressions 

0 A’ AG (1.8.12) 

(y—A8) (y—A9) (1.8.13) 

are usually known respectively as the ‘sum of squares due to regression when 
01, O2 ---» Om are fitted’ and the ‘residual sum of squares’. 


4. Tests of Hypothesis regarding Parameters. We now assume that Y1, Yo, ..., Yn 
are independent normal vari&tes with a common standard deviation c and 
expectations given by (1.2.1) Now 
(y— A80) (y —A8) - (y — A8) (y—A8)--(6—0) A’ A (68—0) 
=(y— A60) [I— A, (45 A,) 1 AL] (y — A0) + 
(y—A6)’ A (A^ A,) A^ (y—A 90) (1.4.1) 
Remembering (1.8.8) and (1.8.9) and by Gochran’s rule the first and second 
constityents on the right hand side are distributed independently as x? o? 
with n—r and r drgrees of freedom respectively. Hence for testing the 
hypothesis 0 = O, we may apply the F-test where 
- Lo y’ A, (A5 A) LA^ y/r 
(y—A8)y (y-A6)/(n—r) 


A 9 A ylr 8 
(y' y—6' A’ y) | (n—7) 


with d. f. r and n—r. 6’ A’ y can be easily obtained by pre-multiplying the 


(1.4.2) 


left hand side of (1.8.2) by 8. Let 
Z= [us oy Sa HZ = (2525 tm Ze = ni Samy os Bp fe 
The rows of G A’, where G is given by (1.8.7) are orthonormal. For 
G A’, (GA) —GA',A,G'—G (A, Aj) (A,A) Gc 1—G G-1—T 


Henee we can find further n—r orthonormal vectors L, in an infinity of 


GA‘ : 
ways, such that ( 3 | forms a complete orthonormal system. Writing 


[T coe 
zZz-(y—A9) (y—4A9) 


y—A9) A (A^ A) A’, (y—A 0)=[G A'1 (y—A 9)]' [GA] (y—A 9)] 
—Zz,—5 kaka ae (1.4.8) 


( 
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(y—A 6)’ [I-A; (A^, A,)7* A4] (y—A0)=2'n Z, = 22, tate be Em 
(1.4.4) 


24, Z5. ..., 2, are independent normal with expectations zero and variance o°. 
From (1.2.10), the best unbiased linear estimate of the estimable parametric 
function b’ @ is b'; (A’, A,)~1 A’, y=b’ 6+b’, G’ GA’, (y—A 6)==b’ 0+ 
ba G' z, and depends on z}, Z» ..., 2; and not on Spt» +++) 2,4, Hence to test 
b' 0=0, we may apply t-test where 

b^; (AA) A^ y[ /b'; (A, Ay)? b, 


VI y-8 A' y) | (n—r) ` 
with (n—r) degrees of freedom. e 


nen (1.4.5) 


5. Test Involving Several Linear Functions of Parameters. Let b'; (8). i=1, 
2, ...,kbek estimable parametric functions and let the first k’ of them be 
independent, the rest depending on them. Let b;—(b;, bios ..., Dg, «+s Dim} 
and ba—(55, bie .-., bir}. Then by (1.2.10) big (A', Aj)? A‘ y will be the 
best unbiased linear estimate of b, 0. The dispersion matrix of the first X’ 
of these estimates is the non-singular matrix 
Bi (A A)! BY 
where B1— (bi, b’s,, Diah. 
When b; 6—0, i=l, 2, ..., i 
Y Ay (A^, A4)? BY [B, (AA)! B'1]-* B, (A4 A) t A^, y 

= (y — A0) A; (AGA) 7 BY, [By (4A) 1 B',]-! B, (AA) A^ (y — A9) 

—Z,GB*,[B;(A,A,)! B/,]-1 B, G'z, (1.5.1) 


is a quadratic form in 2,25, ...,2, of rank kh’. 


A 


G B^, [B,(A',A,) ! B^,]-! B, G” is besides indempotent. Hence (1.5.1) is 
distributed as y?c? with hk’ degrees of freedom. Hence for testing b';0 = O, 
i = 1,2, ..., I, we can apply the F-test where 


p y' A, (44,A,) ! BY, [B, (A5A,) ! BY] B, (AA) ? A’, y/k' 


= (1.5.2) 
(y'y — 8 A'y) [n — r) 


with d.f. k’ and n — r. 


Let B = (bij, bra an; bia}. The dispersion matrix of the best estimates of 


b’,0(¢=1,2, ...,k) is B (AA) 1 B’. Let v= fv. Ure vl. Then, if we 
write down the equations 


B(A5A,) A, y = B(A^A;)-! B v, (1.5.8) 


the expression V B(A,A,) 1A’ y will have the same value no matter 
which solution of (1.5.8) we take. The numerator of (1.5.2) is proportional to 


the value of v' B (A5A,)71 A^ y corresponding to Oya = Uppa m.m Ug = 0. 
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Hence obtain any solution v of (1.5.8), pre-multiply the left-hand side of 
(1.5.8) by v'. Then (1.5.2) can be written in the form 


e po VBA Any 
(y'y — 9' A'y) | (n — 7) 
with d.f. n, = k', ny = n — T. 
6. Tests of Subhypothesis. Method of Fitting of Constants. 
Let 0, = ( 0, 0, .-., 6, } and 6; = (89,5 91m e Dm) 


Let K:N 
A 0 = (A; Ay) and let Rank 4, — s « t, (1.6.1) 
2 
Let 6, (0,) satisfy 
A's (y—A 0.) = A's As 6, (65) (1.6,2) 
Then ( 6, (8) 
6 Alyy =A‘ A (1.6,8) 
6: 
` = From (1.8.2) 6, 
A'jy — ASA rh (1.6.4) 
6, ° 
From (1.6.8) and (1.6.4) we get ^ 
ô, — 6, (09) 
O= ALA - (1.6.5) 
0, — 0; 


Now (6—6)’ A'A (6-8) 


E 6, ri 8, (82) E 8,(0,)—6; | H A | | n 
0, — 0; (0) 6, —0, 
6,—6,(0,) UM [sd T hod eu [er 
ere: 6,—9, o 5 


(1.6.6) 
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From (1.4.1) and (1.6.6) we get 


A E 6,—8,(006.) | ' ( 6,—6,(0 
(y—A6)‘(y—A®) =(y—A8)'(y—A6) + s 1 2) ADA. y " 2) 
6,— 0, 6,—0, 


6,(8.)— OOR 
A p | 4 [80 | da 
ur MI t 


Let the first s columns of A; be independent, the rest depending on these s 
columns and let the matrix of the first s columns be denoted by As. Then 


us e EE] 3 COEN) A A (88)—8)) 


o 
—(y—A9)' A; (A'; A;) 7 A’; (y—A0) (1.6.8) 
This is a quadratic form of rank s. The first constituent on the right hand 


side of (1.6.7) is a quadratic form of rank (n —7). The rank of the second 
form is at least r — s. The second constituent is 


(y— A9/ [A,(A’, A,) TA — AS(A' 5 A;)) 2A’, ] (y— A0 1.6.9) 
since A, (A’; Aj) 14', — A; (A’, A;) 1 A’, is subject to n — r + 8 
independent restrictions. 

A’, [A4 (A4A;)^! A', — A; (A'5sA;)1A,] =O (1.6.10) 

T [A (A A,)1A*, — A, (AA) A] = O (1.6.11) 
where t = {t,, t» A and AT =O, the rank of (1.6.9) is at most (r — s). 
Hence the rank of (1.6.9) is exactly r — s. The constituents of (1.6.7) are, 
therefore, independently distributed as y?c? with n — r, r — s and s degree of 
freedom respectively. Ths first constituent is independent of 0. The second 
constituent depends on 8, and not on 0, and the third both on 0, and 6,. 
Hence for testing 0, = 05, we apply the F-test where 


6, —6 e 6, — 6, (02 
| IS OM APR | KN AEAN, 
F= 9, — 03 J 9, — 6; 
(y'y — & AY) I (n —7) (1.6.12) 


with d. f. ny =r — 8; ng =n —r 


When 6, — O, the numerator of (1.6.12) is proportional to 


$ o Sorin E 3 
5; | AA | a | — 8r (0) A, A,ô, (0) 
6, 6; 


—6'A'y —6,'(O) A,y (1.6.13) 


` matrix of y and E (y) be gi 
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= Sum of squares due to regression when 0, 0, ..., Om are fitted — Sum of 
squares due to regression when only 0;, 0, ..., 0, are fitted. 
In this ease z 


(6 A'y —6, (0) A'sy)) (r — 8) 


= - (1.6.14) 
(y y —9' A' y) / (n—7) 
with d. f. ny =T — S; n, =N —7. 
Itis evident that no test is available when r — s. 
Often- 


7. A Particular Test Involving Equality of Some of the Parameters. 
times, we are required to test 0,1, = 0,19 = -o = Om- Let A, consist of a set 


of 7’ independent columns of 


an dg ap ap F daka Tos m 
Gay Qag ee Goy Gotya F ores Toc 08m 


sss. An Ana F Aniti To F nm 
the other columns depending on these columns. Then 
(y — AO)’ (y — A9) = (y — A9) [1 — Ay (AGA) A] (y —- A0) 
+ (y — AG) [Ay (AD Ay) AS — Ag(A'5A,) A's] (y — A9) 
+(y— A0)' Ag (A'5A9) 7 A's (Y — A9) (1.7.2) 
The constituents on the right hand side of (1.7.2) are independently distribu- 
ted as x? o? with n —r,r—r' and 7’ degrees of freedom respectively. When 
Oi = O42 = = @,,, the second constituent is simply 
y' [Ay (A5A1) A — As (A's Ag) 1 A'e] Y 


Hence the appropriate F-test is 


a (y! Ay (ADA) A5 Y TY! A, (A'5Ag) TA y) rr) Gee 


(y'y = 6! A’ y)/ (n a, 7) 
with d fn erc rim ent ;' denoting the rank of (1.7.1). 
8. Linear Estimation with Correlated Variables. Let V be the dispersion 
i ven by (1.2.1). We assume that V is non-singular, 
n-singular H such that H' H = V-1, The case of 
duced to the case of independent variables by 
H'1iyory — H'z. 


We can, therefore, find a nor 
correlated variables can be re| 
means of the transformation Z — 


B ; tric function b’ 0 in the same way as 
an estimable linear parame < is ano 
We define and sufficient condition as the availability 


in secti and obtain a necessary 
of aee that A/c — b. This gives (1.2.6) or the equivalent condition 
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(1.6.7.). The variance of the estimator c' yis c' V c and in order that this 
may be a minimum subject to A' c — b, 


"We —AUore-— VA (1.8.1) 
D= AV AX (1.8.2) 
where 


eX. O) | 
o IJ 


< Rank (A’, b) = Rank (A^) = Rank (A' H^) = Rank (HA) (H A) = Rank 
(A' V71A) : 


-. Rank (A’ V-1A) = Rank (A' V~? A, b) 


Rank (A'V-1A) < Rank (A' V-! A, b) = Rank (A/,b) ( 


and (1.8.2) is solvable. No matter which solution of (1.8.2), we take c in 
(1.8.1) is unique. For let A and A? be distinct solutions of (1.8.2). Then 


O — A' V71 A (AM — 2%) 
c(0 = V-1AA0), c9 = VIANO 
(c? = cy vi (c? = cà») = (4 = 2?) AL Var A(A® —2?)) =0 


2) 


rec 
GANGGA I) 15 
A particular solution of (1.8.2) is A = Ex Y o 1) 3 


Hence from what we have proved above 

€ —V-! A,(A^5, V7! Aj)71 b, (1.8.8) 
For all c such that A' c — b, we have 

ce Ve=c'[ V — A, (A^, V= A)? A45] c +b’, (A; V71 A)T b, 
—CH'[I—H'? A(A',; V1 4)71 A, H1] H cb’ (Ah V7: A3)! b, 

=T T +b’, (A, V-1A,)71 b’ 

Where T is then x 1 matrix 
[I —H' A, (A^, V~ Aj?1 A, HO] Hc 
"o C Vez b^ (A V^1 Aj) 1 b, 

and the minimum is actually assumed for (1.8.8). 


It is easy to see that the best linear unbiased estimate of b’ @ is 


^ 


b'6 (1.8.4) 
where 6 is the value of @ which will minimise 
(y— A 0) V- (y — A 0) (1.8.5) 
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Differentiating (1.8.5) with respect to @ and equating the derivative to O, 
we get 
o A'V-1y =A’ V-2A6 ` (1.8.6) 
Now from (1.8.1), (1.8.2) and (1.8.6) 

ey SAX A'V-Ly 2A AI V1 A = b'6 (1.8.7) 


. Also from (1.8.1), (1.8.2) and (1.8.3) 


V(c'y)2c' Vc 
2 =’ AVI AA=A'b (1.8.8) 
"mpbu(A VA)! b; (1.8.9) 


A linear function I’ y will be said to belong to error when E (I’ y) = O irres- 


pective of the values of the parameters. Thus if!’ y belongs to error, A'1— O. 
There are n — r independent linear functions of the y’s belong to error. The 
best unbiased estimates are c' y where c = V-1 AA. Covariance between 
c'y and] yisi Vc —l'Ay — O. Thus the best unbiased linear estimates 
ith linear functions belonging to error. 


are uncorrelated w 
^ ^ 
0 we take, H A 0 remains the same. 


© 

In (1.8.6), no matter which solutions 
Hence 
9 A'V18 (1.8.10) 


remains the same. But since 
Ab) VOX (y Ae) =v Vy WE Verg (1.8.11) 
5 ; 3 i 
it also remains the same no matter which solution 6 of (1.8.6) we take. A 
partieular solution of (1.8.6) is 


4 —1 A. po NA V-1 
d 2 a 4 (1.8.12) 


Using (1.8.12) 
(y— A9) V^! (y — ^9) 
Ho [I — HA, (A,;A,) 1 ASH]H !(y — A0) (1.8.13) 


with zero means and unit 


—(y-—A9/ 


Hy A 0) are independent normal variates 


variance and 
1—HA, (A5 Aj)! AH" 


is à symmetric indempotent matrix of rank n — r. Hence 


B(y— A9) V^ (y —^98)— n—rT (1.8.14) 
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EXERCISES TO CHAPTER 1 

1. Given three independent stochastic variates Yy Yo, Ys, having a common 
variance c? such that E(y,) = 0, + 0), Elya) = 0, + 0, and E(ys) = 0, + 0, 
show that 1,0, + 1,0, + 1,0, is estimable if and only if 1, = l; + 1,. Find an 
unbiased estimate of c?. 

2. Let y; i= 1, 2,..., 8 be independent stochostic variates having a- 
common variance c? and expectations given by E(y,) = 0, + 05, E(y;) = 0, + 05, 
Elya) = 0s + bo Elya) = 9, + 05, Elya) = 9, + 07, E(y,) = 05 + 07, E(y;) = 
05 +0, and E(y), = 0, 4-04. Show that 0, — 0, is.estimable and find at least 
four unbiased estimates of 0; — 0,. Find also the best unbiased linear esti- 
mate and its variance. Establish further that 0, + 0, is not estimable. 
Obtain an unbiased estimate of c?. 

8. Assume that y; i = 1, 2, ..., 6 are independent and have a common 
variance o? Given E(y,) = Elya) = 0, Elys) = Elya) = 6, Elys) = E(ye) = 
0,--0,, obtain the best unbiased linear estimates of 0, and 9,. Find four 
orthogonal linear functions each belonging to error. Hence represent 

6 o 


S yi as a sum of squares of six mutually orthogonal linear functions of the y's. 
i-1 

4. Itis known that the period of oscillation ¢ of a pendulum is propor- 
tional to the square root of its length l. In order to obtain the constant of 
proportionality, an experimenter chooses lengths lj, la ..., lų of the pendulum 
and makes 7, Ng ..., ny observations respectively corresponding to these 
lengths. Assuming that the observations are equally reliable, obtain, the best 
unbiased linear estimate of the constant of proportionality and an estimate 
of its variance. 

5. Given o< 4; <m anda; Æ «; for i + k and y; are independent variates 
having variance c? and E(y;) = 0, cos «; + 0, sin ap 4 = 1, 2, ..., n obtain the 
best unbiased linear estimates of 0; and 0,, their variances and the covariance 


between them, Show that 21 (y; cos aj — y; cos ai) (a) is an unbiased 
2 
jmd 


estimate of 0, and comparing its variance with the variance of the best 
unbiased linear estimate of 0, obtain the trigonometric inequality 


geld» | [an] SAS aan) 


i=l 


where the summation Z' extends over all values of j except i. 
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6. E(yy) = vin à =1, 2 aki) = L, 2, n Making usual assumptions 
about y; obtain separate tests for ' 
(i) ua = uam = Pk 
and (ii) py = us = "= Ue = 0- 
7. y's are independent normal variates having a common variance o?. 
(i) Given E(y;) = u, i = 1, 2, «7, obtain a test for p. = o. 
(ii) Given E(yj) = pry t= 1, 2, u ny and E(yj) =p, $ =n +1, san5 
n, n obtain a test for py = pa. 
8. E(yg) =u +u t 5;* = 1,2, .4,5;j =1,2, .,v. Making usual 


assumptions about yi;'s 
(i) Show that 1, mi + 1g tao +1, % d$ estimable if and only if 


1-1. +1,=0 


(ii) Obtain a test for Tı = Ta = + = Tp = 0. 
9. E(yy)=e+a +748 (aj — z) i=l, 2 DIF = 152,50 055 
^ p 4 


are known constants, a = > x; | pq and yijs are independent normal 


fol jel 
having a common variance o?, Obtain the best unbiased linear estimate of 


B and a test for p = o. Show that 1, 7, + lg ++ + loto is estimable if and 
only if 1, + 1, + + + 1, = O. Obtain the best unbiased linear estimate of 
and its variance. Indicate a-test for the hypothesis +; =T=... =1)=0. 


717773 
10. Given known constants aj; (t=, 2, .. k;j =1, 2, -.., mj) 
Th (E : 
$;— S ty [noa = x nj (m + na cbe + n) gus to be independent 
jel i-1 
normal variates having a common variance c? and E(y;) = «; + B(aj— a 
(a) Obtain a test for the hypothesis B = Bs = «+: = Pr 
(b) Assuming 0; = Ba = = By = B (Say), test 
(G) ai = a + Bm (Ti — 2) 
(ii) B = Pm 


11. Ely) —«4-B(v; — a), (6 — l, 2 e) a's'being constants, @ = (a + 
wa +... + a@,)/n and y's are independent normal having a common variance 


2 
G^. 
(i) Obtain best unbiased linear estimates for « and 6. 


(ii) Derive separate tests for & = 0 and B = o. 
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12. Ji Yos --- Yn are independent normal having a common variance c? and 
Ely) =a Bian — 23) + Baltai — 22) +.» + B, (rp; — a), as being known 
constants, 2, = (a + e 4- ... +a,)/7,t = 1,2, ..., p obtain best unbiased 
linear estimates of «, Bi, Bo, ..., 8, and find their dispersion matrix. Obtain 
tests for (1) 8, = f; and (ii) B, = B, =... =. 


18. Ji Yos -Yn are independent normal having a common variance c? 


and Ely) = «o Xo (v;) + o% X, (a) +... + ap X, (2), p <n and Xo =1,, 


X; (x) for j — o. 1, 2, ..., p are orthogonal functions so that for j Æ k 


n 

SS (2;) Xy (mj) = 0, Tis 25. ..., Vp being a set of constants no two of which 
iei 
are equal. Obtain best unbiased linear estimates of aj, i = 0, 1, 2. ..., p, and 
their variances. Derive a test for a, = o. 


14. Let yy Yo, ---, Yn be a random sample of size n drawn without replace- 
ment from a finite population consisting of the elements Y,, Y, .., Yy. Let 
y = (Y,4 Yat. + Yw)/N,o? = (YG --Y2-4-..--Y& = NYI) [NG Show 
that (i) E(y;) = Y, i = 1, 2, .... n (ii) Dispersion matrix V of y = (us Yasan) 


"m INE! 1 z 
is [o? /(N — 1)] (N I, — E,,,) and V~! = F (+ IRE Em | . Find the 


best unbiased linear estimate of Y and an estimate of its variance. 


15. Yu Yz -+ y, are stochastic variates having dispersion matrix V = diag 
(L/w, lwa, ..., 1/w,)-and E(y)-— «--8a;, (i=1, 2, ..., n) the z's being 
constants. Find the best unbiased linear estimates of « and B, their variances 
and the covariance between them. Work out the particular case a; = i, 
w Nt TN, 2; m. 


16. Let a finite population consisting of N = N, +N, + ...+-N;, elements 
be divided into b strata, the ¿i th stratum consisting of N; elements 7 = 1, 
2, .., k. Let Yi, (6— 1, 2, .., k; j= 1, 2, ..., NG) denote the value of a 
characteristic of the j th element of i th stratum. 


Ni Ni k k 
Let u; = x YN; o;— ` (Yg — Yiy|N,Y = Si NiYil SUN. 
j=l j=l i=1 i=1 4 
Let y — 4 Yir Yao o Vins Yar Yan -s Yang c Yki Yeas s Yin, } be a random 


sample from this population, the sampling being without replacement and 
there being n; elements from the ith stratum and letu = {tas sy ey 
Show that 


(i) E(y)— diag (E, 5, E, e E,, 1) u 
(i) VY) = diag [N in, — E, s Jot, 1) «5 (Ney, — Es) ot] (N, —1)] 
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nah > N,—1 f 1 Nz—1 
=] 1 ' 
(iii) V— (y) = diag [ Ne Un "xa T NS 


Find the best unbiased linear estimates of 
k 


(NG pa Nata + +» NR WISAN; and obtain an unbiased estimate of its 
. i=l 

variance. d 
17. We are interested in estimating the population mean of a character y. 
For increasing the precision of the estimate, we utilise an auxillary variate 
æ which is highly correlated with y. Suppose that the finite population of y 
is divided into k classes according to the value æ; of æ, i = 1, 2,... k, there 
being NN; values of y corresponding to 2 = v; N=N,+N,+...+N, A 
random sample of size n is drawn, from the population, the number of units 
of y having the value æ; being n; Assume that (i) Elya) —a-- Bao 
(ii) Var le =a) = yt; (iii) the sampling within a class is carried out, 
without replacement. Find the best unbiased linear estimate of N(« + Ba), 

k 

m DFN a;|N. Obtain the variance of this estimate and also an unbiased 


i+] 
estimate of y. 
ig. A finite population is divided into M clusters, each cluster contain- 
ing N elements. Let a sample of mclusters be chosen at random and from each 
selected cluster let n unitsbe chosen at random without replacement. Let Vg 
denote the value of a characteristic for the j th individual chosen in the ith 
selected cluster (i = 1, 2, > m; j =1, 2, ... n). Denoting by c?, the square 
of the population standard deviation, p, the population intra class correlation 
coefficient and by y ={Y11 -++ Yani -Ymr Yma «-Umn}» the observed sample, 


show that 


fuir nt PAP Pa 
Oyo RO ee ais Pede ant 
AA vj [BBB A 


where 


V, = al, + bEnn Va = CF nm 
M a -- nb + nc (m — 2) 
A=(1/@)In = —|-z-- (a--mb — nc) (a+ nb — ne + nem) 


B= — BEqn| (a + nb — ne) (a + nb — ne + ncm), a = c? (1— p), b =p c? 
ande — — c? [1 -En-—160]/N (M —1). 
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Obtain the best unbiased linear estimate of the population mean and obtain 
its variance. Derive also an unbiased estimate of this variance from the 
sample values. : 


19. Letz 2... v, be an ordered sample from a rectangular population 


with centre p and range >. The standard deviation c is evidently (2/12). 
Show that 


(i) Ela) = u + ca; where a; = 4/8 (2i — n — 1) ] (n. + 1), 
(it) V = e? (wy), wj =12i(n —J+1)[(m+1P2m+2),i<f 


ejes MO Wks A) 
Sh Fest Wy. X9 
(i8) [1262/(n+-1) (n+ 2) V-1 = QOS B= uo. 0 


2 


Sey, 


Find the best unbaised linear estimates of y, and c and show that the two 
estimates are uncorrelated. 


20. Let z,4,..7, bea random sample of size n from the rectangular 
„Population having centre u and range = 24/8 c and let X, Xa... Xp be the 
censored sample of size kin which Ajis the greatest observation. The (n—k) 
censored observations are known to be greater than X, Let the k observa- « 
tions be arranged in increasing order of magnitude, the r th smallest value 
being Xpy. Thus j 

Xy, < Xan < Xin 

Show that V—!, the inverse of the dispersion matrix of Xp Xam Kaya is 


2-1 Oe OM Ome 50 


(B'S ed COA AD aS al eke eee, a ENS 
12 02 O° 0r 20 — quor 


Obtain the best unbiased linear estimate y. and c from the censored sample 
and find the variances of these estimates and the covariance between them. 


21. Develop the theory of Linear Estimation when the parameters are 
subject to linear restrictions. 


CHAPTER II 
GENERAL STRUCTURE OF ANALYSIS OF DESIGNS 
1. General Theory of Analysis of. Experimental Designs. Let there be v 


treatments and b blocks having hy, ks ; ky plots respectively and let the 
ith treatment be replicated 7; times (7 = 1, 2....v). Evidently 


Tag Ta e Hro = kr ke-t a. H ly =n (2.1.1) 


Let nj denote the number of times the ith treatment occurs in the jth 
DIGI (TE y Dive ege E19 t b). Then 


N — [n] (2.1.2) 

will Be called the incidence matrix of the design. Now 
EN = (ky, ho, <., hy) (2.1.8) 
NE, = {rp Tos ge To} : (2.1.4) 


Let the total yields of the blocks be denoted, by B,, By, ..., By respectively 
and the total yields of the treatments by 7, T, ..., T, respectively. Let 


B = (B, Bs, ...,B,} (2.1.5) 
T-—[T,7,..4,7) (2.1.6) 
G=E,,B=E,,T (2.1.7) 


We shall assume that the yields of the plots are independently and normally 
distributed with a common variance o? and expectations given hy the sum 
of a general cffect yu, the effect of the block containing the plot and the effect 
of the treatment applied to the plot. Let the effect of the blocks and the 
effect of the treatments be denoted respectively by 


H a= {an QU 2e ay} (2.1.8) 
c= (ts Aste To) (2.1.9) 
In this case the normal equations (1.3.2) are 

G nly, E diag (hy, ..., 15), E, diag (74... SE) 


B | =| diag (ki, ..., 4) Ej, diag (A, .., ky), N' (2.1.10) 


IIR? E? 


ai L diag, (74, ...,7,) Eg. N, diag (r4, ... , 7) 


. 
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The rank of the first matrix on the tight hand side of (2.1.10) is the same as 
that of 


[s (Es, kas a.a, kg), N' 


‘ (2.1.11) 
N,  diag(ri,..., Tol 
and the rank of (2.1.11) is the same as the rank of either 
diag (ki, ka, a., ky) N’ 
(2.1.12) 
0 [0] : 
or 
D N’ 
f (2.1.13) 
O diag (n, Ta, ..., D») 


where 


5 


n H 1 1 
C = diag (r, Te ee To) — Ndiag 6 B ip. z) N’ (2.1.14) 
1 b 


* 


1 
ae zil N (2.1.15) 


1 
D 
1 

tn 


D = diag (kik, ... 1) — N' diag ES 
1 u 


(2.1.12) and (2.1.13) can be obtained by pre-multiplication and post-multipli- 
cation of (2.1.11) by the non-siñgular matrices 


L, o D 


bs QE 


—N diag i nee a Su LC — diag E Mas zd NI, 
respectively. Hence we get the result : 
b +-rank C =v + rank D (2.1.16) 
C and D are obviously symmetric matrices and 
CE, =O, E, G =0,DE,, =0, E,D=0. (2.1.17) 
The rank of the first matrix on the right hand side of (2.1.10) with an addi- 


oO f ^ 
tional column E ] where L = lis lase; ly } is the same as that of 


2.1.18 
L O, GL D ( ) 
In order that (2.1.18) may have the same rank as (2.1.12), it is necessary that 
L+l +... +4, =0. Hence a necessary condition for the estimability of 
L « is that E,, L — O. Such a function of the treatment effects will be 
called a treatment contrast, When rank G — 7—1, every treatment contrast 


I diag (ky ky... 5), NI, O ] 


D 
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will be estimable. A similar consideration will show that a linear function 
of the block effects in order to be estimable must be a block contrast and 
when rank G = v—1, by (2.1.16) rank D = b—1 and every block contrast is 
estimable. Incidentally, it may observed that if rank G =v — t, a set of 
(t — 1) independent treatment contrasts L' c which are not estimable can be 


obtained from 


(€ ] L=0 (2.1.19) 


A necessary and sufficient condition for every block contrast and 
treatment contrast to be estimable is that rank G — v — 1. When this 
condition is satisfied, the design is said to be connected. Pre-multiplying 


(2.1.10) by 
1 1 
[ O, —N diag [z Pee Ej 1 


we get 
e E N di m Sb i 7 B 91.6 
Q—T-—Ndiag DER b ate iy J (2.1.20) 
= Ct (2.1.21) 
Evidently * s 
E (Q)=Gz, V(0)—cG E,Q—O. (2.1.22) 
c ES 1 È 1 1 
Similarly pre-multiplying (2:1: 10) by [o. I, — N diag [ ke ] ] 
1 v 
we get 
z P's 1 Ik os 
P — B — N' diag [o ae 3 T (2.1.28) 
— D (2.1.24) 
Obviously 
E (P) = Da, V (P) = ? D, E; P= O (2.1.25) 


Q and P as defined in (2.1.20) and (2.1.23) are known respectively as 
adjusted treatment yields and adjusted block yields. We can also easily 
verify that in order that adjusted treatment and adjusted block yields may 
be mutually orthogonal, we must have 


o one 1 1 ; 
C diag rot rA I N=0 (2.1.26) 
4 1 1 
or N diag jg a Te D=0 (2.1.27) 
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Thus in order that the adjusted block yields may be orthogonal to the 
adjusted treatment yields, the design is either (i) not connected or (ii) the 
incidence matrix N is such that ny/r; is constant for all j. This also leads 
to the condition ^g; = constant for all Jj. Hence consistent with 
conditions of the design no Mij can be zero in this case, Defining an 
incomplete block design as one in which at least one of the blocks does not 
contain all the treatments, we can assert that in a connected incomplete 


block design all the adjusted block yields can not ke orthogonal to all the 
adjusted treatment yields. 


In a large variety of designs, we have k —d—.o— ky, = k (say) and 
each ny is o or 1. In this case from. ( 2.1.14) we have 


: 1 
C —diag (rj, 7 n) — TaN Ni (2.1.28) 

e 

Equation (2.1.21) now can be re-written as 
ANE T RTI py UN 
Mr [x3] Lr LE a 2 Mica «c fug 
m =e rers Tp E (2.1.29) 
b 


Where A,, denotes the number of times the $ and s' th treatments occur 
together in a block. A particular Ket of solutions of (2.1.10) is 
^ , 1 1) SN TOR 1 

=Gjn, a=diag baa. | pop ^. SS 
1. —G[n, a —diag (= m Ey; p—diag LE =| 
Hence the sum of Squares due to regression when H: 045 Oy, 
are fitted is 


Ginte B $^ T 


: TET 1 D 
— B' diag E29 *) B+ Q (2.1.81) 


The associated degrees of freedom is b +o— t. Whent = 
set of solutions of the new normal equations are’ 


^ ^ / 1 1 G^ © 
e (0) = G/n, a (0) = diag E 5 rJ B TM v (0) (2.1.32) 


Hence the sum of Squares due to regression when only u, 
fitted is 


O, a particular 


Ws ho, .... O, are 


: 1 1 
B’ diag E: s x) B (2.1.83) 


The associated degrees of freedom is b. Hence if r= 9. 


; the difference 
between (2.1.31) and (2.1.33) 


| 
J 


Date JR. En eir. 


Acc. No.L& QS... 
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v 


s=1 


4 
*. Qs 


(2.1.84) 


is distributed as Z? c? with v — t degrees of freedom. The analysis of variance 
for testing treatment and block effects may be given in the following 


summary table. 


TABLE 1 


Analysis of Variance Table (Intrablock Analysis) 


a 


Source DF ss SS DF Source 
Blocks b Blocks 
(ignoring B G (Eliminat- 
en ba x SANUS Z + as ANG 
treat- c kj n ing treat- 
ments) dm ments) 
[US 3 o v x) 3 ^ 3 
Treatments 5 "e D @ Treatments 
(climinat- | v—¢ za ay oe ala! (ignoring 
ing blocks) FEX zm blocks) 
Error in —b — o-- 1, T > n—b—ov--i| Error 
Total Total 


T Obtained by subtraction. 


2. Two-way of Elimination of Heterogeneity. In some designs the position 
within the block is important as a source of variation and efficiency is gained 
by eliminating the positional effect. Consider a two-way design with u rows 
and w^ columns. Let there be v treatments, the i th treatment being repli- 
cated r; times. Let l}; denote the number of times the i th treatment occurs 
in the j th row and mj, denote the number of times the i th treatment 


eio: 
CHA 
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j' th column having treatment i (say) be u -+ aj + B; -+ 2, where a; — effect of 
Jj th row, Bj — effect of j' th column and 7; — effect of i th treatment. Let 


(dem CD cu 31:518: Ga; Boy 7... Bur}, T= (855 Tas 7.) ` 
The normal equations are 
G uu' l, Ww Eim uwE,, Ep diag (Ti Tg, " 
! 2-992520) 
= S w Eu, "UIN kN L' ci 
a u Es Erus ul M' B 
T diag (r,, 75, L, MUR diag (74, 72, = 
753 To) Eyi» RRA) 
me (2.2.1) 


The rank of the first matrix on the right hand side of (2-2-1) is the same 
as the rank of 


wl, E, L' 
E, ul, M' (2:2:2) 
L, M, diag (71, 7, ... , 45) 
(2.2.2) has the rank u +w’ — 1 + rank F which is the rank of 
TOPE a LU 
u 
L DIN ha ec 2.2.8 
O, Tu — g Ewu MIT Pe diag (ry LL n) GREN 
O, O, F 
where 
di m ; STI l MM 4 -LE Li 
F = diag (7, Ta ..., qr w T "RT uu 44 (2.2.4) 


The result is obvious for (2.2.8) can be obtained by pre-multiplying (2.2.2) 
by the non-singular matrix 


1 
e UE. [9) [6] 
Ki 
1 1 
9, Pig EE alent (2.2.5) 
1 E, 
== L + JE ww J — My IL 
Ub Www u E 


Note that F as given in (2.2.4) is symmetric and 
ME, =L Eww LE,,, =M E, Ei, F=0, F Ej =O (2.2.6) 
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It is evident that in order that a linear function of the treatment effects be 
estimable, it must be a treatment contrast. When rank F — v — 1, every 
treatment contrast is estimable. Pre-multiplying (2.2.1) by 


L Ea ;— DL ,— M 15 we get 
au w m 


Q—F ¢t (2.2.7) 
where 
GIO ge een iret (2.2.8) 
Wu w uu 
E, Q=0, E(Q)—Fcand V(Q)—c? F (2.2.9) 


The row totals and column totals are orthogonal to Q. Also every contrast 
. of row totals is orthogonal to every contrast of column totals. Hence we 
have the following analysis of variance table. 


TABLE 2 


Analysis of variance for two-way elimination of heterogeneity 


ə 


Source DF ss 
Asis E DW Le 
si 4 2 
Rows u — 1 R; E Me 
ul wu 
j=l 
Ry 
; 5 G2 
Columns w — 1 Gia ey 
2 u uu 
Jr=1 
v 
Treatments "3 
v—t as 
(adjusted) = i 
8-1 
Error uw —u—wu' +1—v+t ii 
Total uu'—1 Ly? —G*wuw' 


+ Obtained from subtraction. 
It is easy to verify that when 


1 1 
E,, = L' diag i ELE xj M (2.2.10) 
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the estimable row contrasts and the estimable column contrasts are mutuaily 
orthogonal. In this case, estimable row contrasts are also orthogonal to esti- 
mable treatment contrasts when 


1 1 
[c l,— L' diag Ux Marte a 4 L’'’=0O (22.111) 
1 [J 


and the estimable column contrasts are orthogonal to the estimable treatment 
contrasts when 


[c L, — M' diag [e =| x M'-0 (2.2.12) 
1 t 


The estimable row and column contrasts are given respectively by 


^ ^ = a 1 1 
w I, — L' diag lL Gal eae Lie (2.2.13) 
1 v 


1 1 
and E L,—M' diag [= reads 3) x] B (2.2.14) 
1 v 
and their best unbiased linear estimates by 
: 1 M 
R — L' diag E Jedes ] T (2.2.15) 
KI, To J 
3 1 1 
and GQ — M' diag [ We =| T (2.2.16) 
. 1 D 


3. Analysis with Recovery of Interblock Information. The yield of the 
treatment in the j th plot of the i th block has in its expectation the term a; 
the block effect, which was regarded as constant in the previous analysis, 
If a; is itself a random variable which is normally distributed about zero as 
mean and standard deviation oa, the analysis becomes slightly changed. In 
what follows, we shall assume that there are b blocks of k plots each and v 
treatments are assigned to these bk plots, the s th treatment being replicated 
r, times ( $—1, 2, ..., v )- No treatment is assigned more than once to a block 
and the number of times s and s' treatments occur together in a block is Asse 
Let n,; denote 1 if the s th treatment occurs in the i th block and o therewise. 
Let y;; denote the yield of the j th plot in the ; th block (i = 1, 2; 
1, 2, ...,k). For fixed ay, da, ..., at Vj S are independently 
tributed with a common variance oj and expectations give by u + a -- ^0) 
where w represents a general effect, %;, the additional effect of i th block and 
^j, the effect of the treatment applied to the j th plot of the ; th block; 
Zp %, -..,%,are themselves normally distributed about zero as mean and 
Standard derivation op. 


very O3 j= 
and normally dis- 
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Making use of the following results on conditional expectation 
ER (#|z)=E(a@) 
Vara = E Var (w| 2) + Var E (æ|z) 
Cov (w, y) — E Coy (a, y | 2) + Cov [E ta] 2), E (y | z) 


we get easily 
E (Yy) = u + 58) 
V (yy) — ei + ei 
f E Jose emu 
Cov (Yir Mop) = Q 0 otherwise 


Hence if 


Y = (Vae Yi» e Yus Yous e Your s Yis Yos s Mu) 
2 


RU 1 E me uot 
"E imd FIG E I EFIE EIS 
V = diag (of I, + os Eps -o 0; Ip + 62 Ej) (2.8.4) 
E (y) = (Ea o «s a) (5) =A pu (2.8.5) 


Where t = (Tis Tas ---, To) &, is an bk X 1 matrix consisting of zeros and 7, V's 


€, X, = ry QV Oy, 0, SAS 
w Eq, 
j te — w’ 
cm nde TING (nia Ey. -s My Ey) 
AN MES 
w — 1€ E 
wa, i (na Eis 2 Na Ey) 
w G 
w 
wR - x Bi 
AUVI y = 
, 
w 
w Q.. S Bg 
: k 


(2.8.6) 
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bkw’ Tw Taw. ES Ty to" 
» 


where B; — 28 ng Bik, B; and Q, having the same meanings as in (2.1.5) 
ied 


and (2.1.20) w = l/c} and w' = 1/(6? + ko3). Since 


Ey; Ey, 0 
an e h 
A = es , in order that rank A’ = rank M^ s. 
a’. a'o i, 
ly tlet- +1, =0. Conversely, if l -+l +... +1,=0 
Ey, 0 o g Ey 
rank a, lh | =rank ag dl =v =rank an 
a, ls ae le a'e 
Mae ae a, h ay 


Hence every treatment contrast is estimable. Now 


w—w' Xig(£0 — zw") w — w' 
71 [ E k ] pins k SON [ Jamin 


| 
cobre dE x eee 


is 


wG +- NN (2.8.10) 
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and is nothing but a’, | Vie [ Gy ts... 0, ] , the rank 


a's 
7 
a, 


of which is v. Hence a particular solution of the normal equations 


E 


G w bkw’ , w E,, diag (ry, Tos =: s To) m 
i = à zm n | (2.2.11) 
wQ + x NB w diag {7,,...,7y) Ey, WC + E NN’ T 


^ w’ 71 co^ 3] 
y = 0, t= [wc +- NN’ ] [we + NB] (2.3.12) 
The best unbiased linear estimate of a treatment contrast L’ q is 
w' 71 w 
Saar) L' E + T NN’ wQ + p NB (2.3.13) 
k 
and the variance of this estimate is 


i E) -1 
L' (se Nn”) L (2.8.14) 


Making use of the transformation y = L’~!x where V^! = LL’ and re- 
membering (1-6-13) we can easily deduce that when v = 0 


M m .G ` 
EE NB Je Em (2:3-15) 


has 72? distribution with (v— 1) degrees of freedom. 


4. Expectations of Sums of Squares in Intrablock Analysis with the Model of 
Section 3. r 


Erom (2:8:11) 


D 


E (G) = Mu $i ns 


s—1 


e dd ies PAK kl 2 
E(G@2/bk) = V(G/Vbk) + [E(G/ Vek) Y = -y 7 { vii Pee ia is 
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j A b k 
E [Total Sum of Squares] = Æ IS D Vu es | 


D 


i=l j=1 
= (bk — 1) 0? + k (b — 1) o2 + S rae EXE 
(2.4.2) 
E (T,) —7, (9. ba), E(T2) = 7, (ey + o3) +72 (u + 2,)2 
a Tie . dor nd Ey. X 
B| SSL Te | = © AH Ua +(e bod + EnEn aab 
gl 
(2.4.8) 
b i b 
2 G2 By s 
B > r pe | = OH We? + hod) + at tu oes trys)? 
i=l i=l 
21, Ta E 
E (2.4.4) 


^ linear function e’y = (e'i e', ..., €y where e = (en, ej ... , ĉn) belonging 


to error in intrablock analysis also belongs to error in the analysis with re- 


covery of interblock information. For these linear functions, we have 
Q, 4 = 

also e'; E,, =0. 

7 : 1 2 2 z 2 

V(e' y) — e' diag (sj I, +03 Ey, ... o? Ij, + 5; Ex,)e =c e'e (2.4.3) 


Tf eU" y and e?" y are orthogonal and both belong to error. 
Coy (e y, eO y) = e% diag (o? T, + 63 Ey, ... , 02 I, + c? Ej) e? 


2 all (2 
= ej e ef =0 + (2.4.6) 


Hence 


E [Error Sum of Squares in Intrablock Analysis ] 


= (bk—b—v +t) o? (2.4.7) 


In the ease of connected desigris, we have the follow 


ing table of expecta- 
tions. 


2t 
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( 


319. LX E VES " 319 1 
uno Weug) THR net Tru 
8? (1—4) y + le (1 — qq) 8? (T= 9) 4:2 — qq) 
39opqeaug) LOL lo( I--a—q—yq) 101 +a—q— 44) TA Cporquaqug) dosag 
SAN ngar aE 


Suður Bony, 


Jos (7 "a )— 
VNN aht (204+ 1e) (Tg) 


n) sxoo[gr 


994008 
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5s MES (2.4.8) 


E ( Vi) = (bk —b — v + 1) c1, E pP LG 


E (V,) == (0—1) 0} + (bk—v) a3 


: V, (b—1) V: art ; 
Z p2 CE a | Gan 
1 A v—k y BE 2 D 
E (gm (EM —€— 93) o ko (2410) 


When the experiment can be arranged in r complete replications, 
E E l (say ). The expected value of the yield in the g th plot of the 
= 


j th block of the i th replication for fixed $;/sis given by 
E (Yj) =v Hot BG bia) (2.4.11) 
ug p qr D aa b edge 


where p; denotes the effect of the i th replication, B; denotes the effect of the 
j th block of the i th replication and tjg denotes the effect of the treatment 
applied to the qth plot of the j th block of the ith replication. We assume 
that y's are independent normal having mean zero and variance o? and Yija S 
for given B,jsare independently distributed with mean given by (2.4.11) and 


variance o2. This leads to 
E ( Yijg = 9 + ei * Tin 


V (Yig) = i + 98 


UN NE EUIS t, 
be Ld cif) —7,) —J'.qq 
Cov ( Yijgs Vigor) = a otherwise 2.4.12) 


With the above set-up, the following expectations of the constituent sums of 


squares in the intrablock analysis can be easily obtained. 
= 
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D 


iid (ay =H] ee 


*]8}0} puwan 
juounwo4 YZ $ Jo TEJO, 
'uonvordoz y} ? 243 jo qv30j, ! 


52 (1 — 4) (y — a) + Yo (« — q) (D) 


$9 (1—@) 4+ 19 — qq) é we L— NG monL 
8 
52 (1 +a — q— 41) “a — "w(Lka-q—30) — r--a—q4— qq aoan 
1 =° 
——. r na SEA Er 1 e 49 4 T (po3snfpeun) 
E "ur) x z] * + 39 (3—2)-- Io(r—a) 29i € ats sjusun gag, 


i (poqsn{pe) 
(9) “iT (4 — q) 4—14 suoivordoa 
urs. sxooper 
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I=} 
19 a 
- 4 suomnvord. 
E T < T L suonvordoyy 
(SS) SSK SS a mog 


[== Si ae 9 a ee 


susisagq aqvajosagy papauuop 


fo 2owmiinA fo sishpuy ur so4onbg fo suns Jo suouvjoodan 


CHG wmv 


o 
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o 


~. Est (sta = ee (2.2.18) - 
Tw 
5. Total and Partial Confounding: Let there be. k treatments and let us 
divide them into m sets os k cach. We assign to a block all the treatments 
of a set and there are r replications. In this case from (2.1.14) 


C — diag (age x Eg, rl— ZAR seen E En. (2.5.1) 
L k k k 

L’ = is estimable if rank G — rank (C, L). Rank € =m(k—1) and every 
intraset contrast is estimable. All interset contrasts are not estimable and 
they are said to be confounded with the blocks. In the notation of Section I, 
the adjusted treatment sum of squares is 
1 A 

^ ^ 959 
Tar BNNB (2.5.2) 


1 
nd - 


and has m (k — 1) degrees of freedom. The unadjusted Block Sum of Squares 
carry mr — 1 degree of freedor:: Subtracting these two from the total sum of 
squares we get the Error Sum of Squares carrying m (k~~ 1) (r — 1) degrees 
of freedom. In order that we get an estimate of error r must be greater 
than unity. 


We may, however, not coufound the sanie treatment contrasts in all the 7 
replications. Let us take up the case where the experiment is laid out in 
such a manner that any intraset treatment contrast in the 7 th replication is 
in the i' th replication either confounded (built up of interset contrasts) or 
estimable (built up of intraset contrasts). Let the v =mk treatments be 
numbered as 1, 2, ... v; let the total yield of the s th treatment be denoted 
by T, the treatment applied to the g th plot of the j th bloek of the ; th 
replication by Tuj and the corresponding yield by yj, Let yir denote the 
yield of the s th treatment in the 7 th replication and Bj; the total of the 
jth block of the ith replication. Let L' t=), Ti-l ta- ++ H late ADS 
estimable in the Ist, 2nd, ..., 71, th replications and confound in (7, + 1) —th, 
(ry - 2) th, ...,. th replications. Let further E (y;;,) = v. + pi tBu + Lega and 
V (Yaim) = °; Cov (Yar Yirjrar) = 0- ; 

Then 


L/c — lj, baa te t ho So 


mo k 
= > X lija a) (2.5.8) 


j=l qe 
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where (2) ij; + lja + --- + ljr =0, for i = l, 2, ..., T, 


and (4) lj; = lja =... = li, = dg 
m 
Si dj =0, i =r +l, n F2, ar 
j=1 
n £e UD k 14 : 
Now — [TEES Vi 
T. da Jur T * Fils) 
1 Ge Uu Wa Qnm k T 
E— 7 log. = = l ay =— Le 
"n 2j zx 5 ija Yija "n 2j > 5 üa Ta) 5n 
i=l j=1 q-1 i=l j=1 g=1 
—L'« (2.5.4) 
Also 
Tı m k r m k 
1 1 
m zi ` lija Yigg = n 1 ` 5 lja Yija — 
2 i=1 j-1 q-1 isi j=1 q-1 
j r m k 
ks ` > ` lija Vis 
irl  jei oi 
t: m 
LTi +T +.. +47, 1 
= rs > 53 di Bi (2.5.5) 


Ti 1 isr +1 j=1 


is a linear function of block totals and treatment totals. Hence it is the best 
unbiased linear estimate of L' v. The variance of this estimate is 


| o? Ti v c v 
cs 
En Ta = Bg 2.5.6 
nm S > : A E ( ) 
t=1 j=1 j=1 


Had L’ « been estimable in all the replication, its variance would have been 


c? 2 c? 2 
= ` S g-- S I (2.5.7) 
i=1 j=1 j=1 


Hence if L’ q is partially confounded, the variance of the estimate is E 

g 1 
times the variance that would be obtained if it had been unconfounded in 
all the replication, It is also to be noted that the best unbiased linear esti- 
mate of an estimable treatment contrast is obtained from the corresponding 
linear function of the yields from those replications of the experiment in 
which it is unconfounded. 
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EXERCISES TO CHAPTER II 


Y. It is proposed to test seven treatments A, B,C, D, E, F & G in 
blocks of three plots according to one of the following plans : 


Plan Block 1 Block 2 Block3 Block 4 Block 5 Block 6 Block 7 Block 8 

aa Era oy ne ovi. fex ae La 
T.U A,B C B,F,D C,D,G D,A,E E,C,F F,G,A G,E,B 

LES AS BI C. B,C,D C,D,A D,A,B D,F,G F,G,E GED E,D,F 

III A,B,C A,C,D A,D,E A,E,F A,F,G A,G,B i 


Which plan would you recommend and why? ^ 


2. From the incidence matrix of a design given below, obtain (i) the 
estimable treatment contrasts and (ii) the degrees of freedom associated 
with the adjusted treatment and adjusted block sum of squares. 


>O Ommo 
On Ome: 


O O m H j 
Hm oo 


0 0 0 


8, Ina connected design, obtain an expression for the average variance 
of the estimates of all elementary treatment contrasts in terms of the non- 
zero eigen values of C. 

4. Suppose that the elementary treatment contrast +, — «,, is estimable. 


Show that the variance of its best unbiased linear estimate will lie between 
20" mas ANA 26?/Amin where Amar and Amin are the maximum and minimum 


non-zero eigen values of C. 

5. Show thata necessary and sufficient condition in order that all elemen- 
tary treatment contrasts may be estimated with the same precision is that C 
has (v — 1) equal non-zero eigen values. 

6. If, in any design, every treatment is replicated the same number of 
times, the total relative loss of information is one less than the average 
number of blocks per replication. 


Show further that the result is valid for designs in which heterogeneity 
is eliminated in two directions. 


7. An incomplete block design is to be used for two-way elimination of 


heterogeneity. State the characteristics of such a design, if all treatment 
contrasts are to be estimated with the same accuracy. 
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8. Using the notation of section 2, show that a necessary and sufficient 
condition that estimable treatment contrasts are orthogonal to both estimable 
row contrasts and estimable column contrasts is 
F diag ( 1/rj, 1/rs, -.., 1/ry) L — O and F diag (1/r,, 1/72, ..., 1/r,) M = O. 


9. If km treatments are divided into m sets of k each and if treatments of 
a set are assigned to k-plot blocks and if there be r replications, show that 
the design is such that the adjusted block effects and adjusted treatment 
effects are mutually orthogonal. 


10. Below is given the plan of a design, there being 12 treatment and 
8 replications. The numbers in the body of the table indicate treatments. 
Identify treatment contrasts which are partially confounded with the blocks 
and indicate the loss of information on them due to confounding. 


Replication 1 Replicatioa 2 Replication 3 


Block 1 Block 2 Block 1 Block 2 Block 1 Block 2 


2,8,5,3,9,12 1,4,6,7,10,11 1,4,6,7,9,12 2,8,5,8.10,11 1,4,5,8,10,11 2,3,6,7,9,12 


11. Form an analysis of variance appropriate to the design whose in- 
cidence matrix N = 2E,, and compare it with ipse of a design whose incidence 
mat:ix is Ej. 


12. Show that in each of the following designs, the rank of C is v — 1 and 
there is only one non-zero eigen value. Hence find the variance of best 
unbiased linear estimate of an elementary treatment contrast 

) (i) N=E, (ii) G=al — (ajo) E, 

18. Defining Q, — T, — Q, — (G@/bk) and other symbols as in section 1 ob- 

tain the following in the analysis with recovery of interblock information : 


G) V(Q), (H) V(Q,), (iii) Cov (Qu Qur) — (iv) Cov (Q Qs) 
and (v) Cov (Q’,. Q’,). 

Show that in the set up of section 3, u + z, is an estimable parametric 
function. 


14. Q and B having the meaning of section I, establish that in the analysis 
with recovery of interblock information the dispersion matrix Q and B is 


j [ [o 0 ] [ fe) o 
o? 3E 
O diag (ki, ks, .... hy) O diag ( K?, k3, ..., k? Y 


15. Four fertilisers 4, B, C, D each at two levels — single dose and a 
double dose, and a control ( no fertiliser) are tested in r block of 10 plots each. 
Eight plots are randomly assigned to the fertilisers and two to the control 
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in each block. Find estimates of the treatment effects and give the analysis 
of variance. Obtain the variances of estimated comparisons between any 
two treatment effects. Break up the sum of squares and degrees of freedom 
due to treatments into (i) between fertilisers (single dose), (ii) between 
fertilisers (double dose), (iii) between doses. 


16. Show that the best unbiased linear estimate of an estimable treatment 
contrast L' c is given by d’ Q=d, Q,+d, Q+.. +d, Q, where C d=L and 


the variance of this estimate is o? d’ Cd = c? d' L. If C = al — E Ev 


show that the elementary contrast z,— 7, has for its best estimate (Q,—Q,)/a 
with variance 2 o?/a. 


17. Show that the equations for obtaining t, in the combined inter and 
intra block analysis can be obtained from the corresponding equations in 
intrablock analysis if we replace Q,, r, and ày by P, E, and A,, where 
P, — wQ, + w' Q,, E, = r, [w+ w'/( k —1)], As = (w — w')àsy and set in 
addition S 1, TT, =0 


s=1 


18. In the intrablock analysis, find the expected values of (7) adjusted 
treatment sum of squares and (ii) adjusted block sum of squares. 


19. v=s, X s, treatments vy, (i = 1, 2, ...,5,; j — 1, 2,..., 59) are arran- 
ged in b blocks with k plots each such that every treatment is replicated r 
times and v; and ov; occur together in X, k, plots where kı — 0 or 1 
according as i = ù or i #7’ and k, = 0 or 1 according as j — j' or j #7’. 
From the relation ( Ej, N) N’ = E,» ( N N’), establish 


(5$; — 1) Mo + ($4 — 1) Aga + (51 —1) (53 —1) X44 — r(k— 1) 


Write down the normal equations in the analysis of such designs and obtain 
variances of best unbiased linear estimates of comparisons of the form 


E Sa $1 sı 
(i) > gc 5 05, P zd (di) S DUE s vg. j xh 
j=1 j=l i=1 ici 
and (tit) vy -+ viy — vip — Vij 


20. Let there be N = n, + na -- ..., +n, plots and let k treatments be 
assigned at random to these N plots, n; plots having the 7 th treatment. The 
yields are independent normal having a common variance c? and expected 
values given by the effect of the treatment applied to the plot. Test the 
equality of the treatment effects and show that the test is most sensitive 
when m = mg =<., Ng = N/k. 
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21. Show that in a connected design, Q; + r; G[n, i =1, 2, ...,v are linearly 
independent. Hence establish that 
(i) Q + rr'/n is non-singular 
(i) (GQ+rr'jn) r 2 Ej 
where r= {fi Ta, ---) Tofo n = Ty F Ta F oes To 

22. Show that in a connected design 

2 = Gn, * —(G-Err/n) -Q 
a = diag (IJI, Ika, ..., IJ) B — (GÍn) By — 
diag (Ika, ..., 1/k,), N'(c-Err'n)-1Q 

will satisfy the normal equationsin intrablock analysis. Show further that 
the dispersion matrix of f is o? [ (c + r r'/n) - 1 — E,,/n] and the variance 
of the best estimate of a treatment contrast L'r is c? L’ (c++ a =) L 


28. (i) When T, and T, are independent, show that the Tinea func- 
tion of T, and T, having least pate is 


1 [s +o A n) / Cs on tg) 


(ii) Inaconnected design of equal block-size k where further N N' is non- 
singular, show that the intrablock, interblock and combined inter and intra- 
block estimates of the treatment contrast L' t are respectively 


LI RE 
v [e+] Q, L'(NN') NB and 


L' [sex ww ] 7" E T NB | 


Use (1) to obtain the linear combination of the first two estimates having 
least variance and show that it is in general different from the third esti- 
mate. Obtain situations where these two are equal. 


24. Show that 


G —N 
i k, = oy IP Bone? 
OET, N' diag (I, ..., ko) pt Dues u^) 
D w ^ LL 1 
(ii) [wa- AN] P rn, n xo 
When r, — 75 =... — 7, =T, show that the average variance of all ele- 


mentary treatment contrasts with recovery of interblock information is 


xb! c 1 
2 [m wa +NN BI | e 
Tace ( b. T NN) WT fe 1) 


CHAPTER III 


STANDARD DESIGNS 


1. Randomised Block Design. Here we have bk plots divided into b blocks 
of k plots each and each of the v = k treatments is assigned at random 
to one of these plots in each block. Evidently for this design 


N=Ey (8.1.1) 
1 1 
= [1-5] ,D=k [1-5 | (3.1.2) 
G 
Q=T— En (8.1.8.) 
Ray I a bw! 
m rut Lice (8.1.4) 


Condition (2.1.26) is satisfied gud the adjusted treatment and block yields 
are orthogonal. 
TABLE 3.1.1 


Analysis of Variance Table for Randomised Block Design 


Source DF SS E(SS) 


1 G? T 7 
Blocks b —1 B m Wb ad gue pc WIE XE 


l T & Dora x IE 
Treatments o—1 $53 CPP (v — 1)c? + bc = ES 


F Jd ; 1 2 G2 r 
Error (5—1)v—1) yy — b TT a B'B + bk , (b.—3)(v = 1)o? 


G? 
Total bv—1 YN ae ans 


=i 


Fn 
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2. Latin Square Design. Yt is an arrangement of u? plots of u rows and u 
columns in which u treatments are assigned in such a way that each treat- 
ment occurs once and only once in each row and in each column. Adopting 
the notations of Section 2 of Chapter 2, we have 

pia eee 
rank F =rank [uI — Eu] = v —1 and every treatment contrast is 
estimable. Similarly it can be seen that all block contrasts and all column 


contrasts are estimable. The best estimates of estimable row, column and 


treatment contrasts 
G G G 
Ro ous C= u Ba and T— -y Ea e (8-2-1) 


are mutually orthogonal. 


Equation (2:2:7) now becomes 


G 
eS -[a-x]e 00 een 
And & solution of (3:2:2) can be taken as 
A 1 
ie (8-2-3) 


u 


Sometimes when the experiment is over, „one is interested in whether 
the additional column classification was worthwhile or not. By equating 
E(SS) for columns, Treatments and Error in Table 8-2-1 with those for 


"Treatments and Error in Table 8-1-1, we get 


2 Ü 1 Uu l 
u(u—l)o tup [1-22] Brut [1-72 | t= u(u—l)ekg 


il 
purt [ M Ba | T 


2 u — 1)? Ep + Column SS in Latin $ 
Hence Est opp = ( y 5 um) atin Square 


Comparison of (8.2.4) with Ep gives one an idea as to whether column 
classification has led to increased precision or not. 


3. Graeco-Latin and Higher Order Squares. Consider m set of symbols 
Ay, Ajo, IOS} Au 
Ay, lo, ++ Aou 
Ami Amo +++ + Amu 
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` T UN 
c) 

(sn) (r—n) Tq = 
n n u 
2 "a — — anali n) 20 DAS 
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where m < u. Let there be u? cells arranged in u rows and u columns. Then 
an, arrangement of these m w symbols in the u?-plots will be called an ortho- 
gonal Latin Square of order m if i 
(i) in every cell there is a combination of one A,-symbol, one 
A,-symbol, ..., one A,,-symbol, the particular combination in the 
cell in the i th row and j th column being donated by 


Arij) Aati = Amt 
(à) Aian Aas) 
E Aog) = Aien) 
Aitu Arua) ec A uu) 


constitutes a Latin Square, t — 1, 2, --- 
-symbol and an A, -symbol (¢ 47’) 


n Ainu 


m 


iii) any combination of an A; 
M 
appears in one and only one cell. 


1n the special case m — 2, the arrangement is known as Graeco-Latin Square. 
When u is a power of a prime, such arrangements are always possible 
provided m < (u — 1). Let yj be the yield corresponding to the (ij)-th cell, 
Ty the effect of Ay (t = 1, 2, -> M5 i=1, 1, ..., U), % effect of the i th row, 
B; effect of j th column and y, the general effect. Using the same notations 
as in the previous section and denoting by Ti; the total yield corresponding 


to Ay, we have the following set-up 


E(yy)=u + a tB + Tai F adj) d- em) (8.8.1) 
and the following normal equations : 
b b m u 
G —wyd usa +e Dh pusi Si te 
i=l j=1 t=1i=1 
R; = upt uod $8 $E $8 dL 550 
C; UA + Sla tub 3B „J= l, 2, s U 
u 
Tu= up + Sa Sib tute te Dy S Ty it 
t=1 i/=1 
£=1, 2, =» th; i=l, 2, ... W (8.8.2) 


The number of independent equations in (3-8-2) is (m+ 2) (u — 1) +1. 
A set of solutions is given by 

ee NAN Ta ee 
u u? u u? u u? 


R> 
| 
| 
3 
E 
Í 
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The sum of squares due to regression when all constants are fitted with 


d. f.=(m+2) (w—1)+1is 


i=l | E j=l t=1 i=l 


Easy calculations will give the following sum of squares due to regression 
when mm; (i= 1, 1, ..., v) are dropped 


G u f u m-—1 u 
2 RI G? (C5 G S T? G2 
=, 4 L— 1 = —— i E e 
E m I as u a uos ME us 
i=1 LU j—1 t=1 i=l 
(3.3.5 
with d. f. =(m+1) (w—1)-+1. The difference between 
(3:8:4) and (3.3.5) 
u 
IT G? 
A TEL eT (8.8.6) 
I wu w^ 
is distributed as y? c? with d.f. =u — 1 when Tm = ... = Tiny —0. 
TABLE 3.3.1 


Analysis of Variance Table for Orthogonal Latin 
Square Design of m” Order 


Source DF SS E(SS) 


u 


Ri E ye 
Rows u—1 S SoS Giu? (u—1) +u [5 - Sx] 
ici 
Columns  w—1 zi Cjju —G*ju? — (wu —l)c? +u [3s - IE +] 
j=l 


A, symbols u—1 Si Ta ju —@iu  (u —1)?-ru | Zrt- CE ] 


EI 


i=l 


A, symbols u — 1 Sn ju — Gu?  (u —1)c? +u [sc eae “| 


i-1 


Me 2 n 25 
2 
A,, symbols u—1 Di Tzu — Gu? (w—1)o? +u [Se ee a? 


isl 
Error (i —1)(w —1—m) he (u — 1)(u — 1 —m)o? 
cabe uy PE QUA HESSE ES 


Total ue? —1 $e ze 
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4. Cross-over Design. From the point of view of analysis, this design is 
nothing but p replications of an sxs Latin Square Design. Adopting the 
notations of Section 2 of Chapter 2, we have in this case v = $, U = $, U' = sp, 
lj = p, my =l, 

L=p Ess M= Estem) 

G A 3p G 

= E,,t= — — gp En 


F=sp I—p Ev 0 =T- sp stp 


a t G 
j wae - Rank F = s — 1 and hence every treatment contrast is 


TO =—— 5 


sp sp 
estimable. The analysis of variance table can now be obtained from Table 


2.2.1. 


5. Balanced Incomplete Block Design. (B IB D). An arrangement of v 
treatments in b blocks of I plots each, (k < v), is known as BIB D if every 
variety occurs once and only once in r blocks and any two varieties occur 
together in blocks. Let N be the incidence matrix of this design. Then 


Ey, NE; = (bk) = Ej NE = (vr) (8.5.1) 
NN'E,, = [7 4- (o — 1)] Ea = N (N'E,) = (rk) Ec (8.5.2) 
From (8.5.1) and (8.5.2) we get 
bk — vr N (8.5.3) 
xv —)=r(k— 1) (8.5.4) 
We have also 
det NN — [r A (o DI Aa Ska FO (3.5.5) 


Hence rank N — rank NN' =0 <b. 'Thus 
b<vorrSk (3.5.6) 


Ge — » —k)z0 


b—v-—(r—k)z0 


Also 


b>v+r—k (8.5.7) 


When b =v, the B IBD is said to be symmetrical. In this case 


det NN’ = [det NP. =7°(r — DES 
pad (8.5.8) 
det N= + r(r—»)* 


hence when v is even (7 — 2) must be a perfect square. 


det N is an integer, 
For symmetrical BIBD 
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1 A 
(r —2) Lies s. | 


N = (r—2)I LA Ez, (NN)! — N'7! N^ = 


N= [ x- 2 s. | I(r A) 


Postmultiplying both sides by N we get 
N'N —(r —a) I + Ep = NN’ (8.5.9) 


Hence in the case of a symmetrically balanced incomplete block design, 
any two blocks have A varieties in common. 


If the BI B D is resolvable (i.e. consists of r sets gÉ = =a tues each, — ? being 
an integer and each set of blocks contains every dm once), the inequality 
(8.5.7.) can be further improved. In this case N consists of r sets of Z columns 
each, where a set of columns is such that 1 occurs once and only once in 


each row of the set. By adding the 1st, 2nd, Ae —1)th columns to the 


ib? th columns of a set, we obtain a column consisting of 1 only. As there 
T 
are r sets evidently 
v = Rank N «b —(r—1) 
zv ri I (8.5.10) 
(i) Analysis utilising intrablock information only. 


A 
Of ei (8.5.11) 


k 
is satisfied by * = ADS Q (8.5.12) 


Rank C = v — 1 and sum of squares due to treatments is 


k " 
— Q'Q. The best unbiased linear estimate of +, — qy is Ed (Q; — Q,) 
y d , 2k | 2o* kr 
with variance c WOES mx 
1 
2 
GRE uritur) m ERES z . The corresponding 
r ` va(k—1) r 1 
METUIT 
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variance in the case of a randomised block design is 2c7/r ( assuming c? to 
be the same in both the cases). Hence the efficiency of B I B D is 


1 
En 
vA k x 
Jedem c us <1 (8.5.13.) 
v 


(ii) Analysis with recovery of intrablock information. 


ug. wo Seated) QU Em, 


wa |) k , (w—v’) 
[se xw] 7 UA + uw (r — 2) ES w'rk 5. | 


w’ ni w' 
ee [se 2 Nw | OEC NB 
e 


k w' (w — w) 4G ] 
EE ag a + NB E, 8.5.14 
= wor d-w'(r—2) [wo ey | a rk ; ( ) 


Best unbiased linear estimate of a treatment contrast L'r is 


k wi 
a ee eB aa ne | (8.5.15) 
122 -+ :w'(r—2) EL k 


and its variance by (2.8.14) is 


BARN t 
won -+ w (r —A) 


LL (8.5.16) 


The variance of the best estimate of an elementary contrast 7, — Ts, is 


2k 2k(o — 1) 
mun He w (r — A) = woo — 1) + w'(r —A)(v — 1) 
2 Kerl) (8.5.17) 


i. v(k — 1)w + (o — k) w' 


kv —1) is called the effective variance and can be approxi- 
v(k — 1)w + (v — kw" 


mated by E, [1 + (2 — 
nis (8.5.18) 


TET ERIS ee 
v= (eiw + (o — o 


k)v] where E, is the intrablock error mean square and 
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The best estimate of the estimable parametric function u + T is 
^ T , 
= = i = (e k)T, — (v — 1)B;, + (k -De | 5 


ENG 
x [7 “hy w, | (8.5.19) 


Denoting by v, and v, the error sum of squares and adjusted block sum of 
squares in the intrablock analysis, we can estimate w and w' by 


bk —b—v+1 f v k 
2 andan 1 amen Sa 


[5-54 Or UU ]! 
3 jane 


SESINING When, however, the B I B D is resolvable, for estimating w and 
w’ we make use of the error sum of squares and adjusted blocks within 
replications sum of squares and denoting these mean squares by E, and 1 Ena) 
respectively, we get easily 
E, (a) "EL E, 

Est. (w — t.) = ——— q^ 8.5.20 

( E, E, (a) —E, | (8:7120) 
r |E, (a)—B, | 
vr(k—1)E,(a) +k(b—v—r+1)E, 


Est. v (8.5.21) 


When the BIBD is symmetrical, the adjusted block sum of squares can 
be calculated directly. For pis expression is 
b 


v Hj E 2 2 
c F sm. == 


i=1 s=1 


k—X we 2 ; k V Aa oy, I 
= 12 — SB. + EUER >p; eNe E) 


For S15; = (NB) (NB) = B'N'NB=B'[r—1)I 4) Bu] B’ 
= (r — à) EB; + AG? (8.5.23) 
Now W,=(v—k)T,—(v0—1)B;+(k—1)G 
e Siwro o E ST? —2(v—k) (0—1) SB; 
+ (o1 SE,- v(k—1)G (8.5.24) 


*. From (3.5.22) and (3.5.24). 
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) a 1 ie 

Adjusted Block SS — — ANCE Ss WE (3.5.29) 
s=1 

Sun of squares due to regression when 7, = 4, =... = pem) 
mr w’ wG? 
is c'(w Q 7 NB) — ME 

wor +w (r—}) : "NP EE 
o aga S Taty W, — Ty We (3.5.26) 

s=1 > E 
and this has a 7? distribution with d.f. = v — 1 
Also to (bk —b —v+1 ) E, has a y? distribution with d.f. = b%— 
—b —v + 1 independent of (3.5.26). 
Hence 
wu, + cw'(r— A) = ; TRU. 
fo ono) 2) Ue W, FD) 
s=1 
(kw + (0 — k)u' e 
v(k —1)w + (v — k) w! i i PREIRESTTAMS r 
Sa 23 (T, +y W, — T, +9 W,) (3.5.97) 
s=1 
has F-distribution witltn, =v — 1, n, =bk —b —v 416 
Hence 
» (T, +y W, — T, +v W)? EM 


s=1 


E, [1 +(v—k)v] 
can be regarded as approximately distributed as an F-ratio with n, — v—1. 
n — 0k — b — v +1. : 


6. Youden Square. Given a member of subsets Si S3 ..., S, of a given 
set S. These subsets will possess a System of Distinct Representatives ( SDR) 
if we can find a set 

D=[ ay, dy, ..., an] 


such that dj, as ..., a, are all distinct and a; belongs to S; (i —1, Dosey 
n). We can easily prove by induction that S, So- Sh will possess a 
SDR if every t of them contain between them at least { distinct elements 
of S, (£ — 1, 2. ..., n). Consider a symmetrically balanced imcomplete block 
design arranged in a two-way table in which the blocks are columns. The 
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sets S,, $,,..., S, constituting the columns are subsets of the treatments 
1, 2,..., v. Between any £ of them there will be ¢ k elements each one of 
which can be repetition of the same element at most k times. Hence there 
must at least be ¢ distinct elements. We can, therefore, find an S D R of 
these sets which will be nothing but 1, 2,..., v arranged in some order. 
Let these constitute the first row. Consider the residual sets Sj, Sy, ... , S7. 
Any t of them will contain i(k — 1) elements, none being a repetition of the 
same element more than (k — 1) times. Hence there are at least £ distinct 
elements and the subsets possess a SDR which is nothing but 1, 2, ..., 
v arranged in some order. Let these constitute the second row. By repeat- 
ing the argument we can show that the symmetrically balanced imcomplete 
block design can by interchanging the constituents within the blocks, be 
converted into a rectangular arrangement in which every row will contain 
the treatments 1, 2, ..., v in some order. Such an arrangement will be called 


a Youden Square. 


(i) Analysis without recovery of interblock information. 


In the notation of Section 2 of Chapter 2 


F— DUUM TUR ECL E,,. Hence a solution of (2.2.7) is 
T T 
A T 
Umm Q (3.6.1) 


Comparison with (8.5.12) shows that the analysis is exactly similar to the 
analysis of the corresponding symmetrical B IB D; only the sum of squares 
due to replications (rows) is separated out from the error sum of squares 
in the Analysis of Variance Table of the corresponding B I B D. After the 
experiment is over, one may be interested in knowing whether this additional 
feature had been worthwhile or not. An estimate of the per plot variance 
of the corresponding B I B D can be obtained from 


Sum of Squares due to replications + (v — 1) (r — 1) By (2.6.2) 
v(r — 1) 

where Ey is the error mean square in the Analysis of Variance Table for 
the Youden Square. If Ey is substantially less than (3.6.2), the Youden 
Square arrangement will lead to increased precision. In comparing the 
relative magnitudes of these two estimates of per plot variance, the fact 
that the error sum of squares in the Youden Square is based on a lesser 
number of degrees of freedom should be taken into account. 
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(ii) Analysis with recovery of interblock information, | s 
. Denoting by yj, the yield of the plot in the j th column of the i th 
"row, we have i 
E| (va | Bis Bas --- > By) =u t pit BiH ta) 
V| (uyl £s» Bo, see s Bo) =o? 
where p; is the effect of i th row (replication), Taj the effect. of the treatment 
applied to (i7)th plot and Bj's are block (column) effects which are indepen- 
dent normal with mean zero and variance c2: Hence 
E(yg) -EE(v;| Boas B.) =p +9; TT 
V(yg) «EV (yg| Bi.. Bo) + v [E (yy | 81, ++ Bo) ] 
i — + o3 ! 
Cov (yip guy) =E Cov (yg, yos | Bi» «++ Be) + Cov [E (y; | By ..., Be), 
i E(yiy | Bay +5 8,)] j 
! =o3 ifi =i’, j +7’ and 0 otherwise 
| The analysis now follows as in section 2 of chapter 2 i 
| 
From ( 8.5.18 ) and Table we have 
i E,—E, 
: v(r—1)E, > 
The best unbiased linear estimate of U + 7, is 


Est y = 


= (2, +0 W,) 


The variance of the best estimate of +, — 7, is estimated by 
2E, 
T 


[1-(» — k) v] 


7. Lattice Design. Let us have Kk? treatments arranged in k rows and k 
columns. Then if we assign the treatments occurring in rows and columns 
to different blocks, we have 2k blocks of k plots each. Again if m — 2, 
(m—2< k—1) mutually orthogonal Latin Squares exist, we can super- 
impose these Latin Squares on the k® treatments arranged in k rows and 
k columns and assign to blocks treatments corresponding to letters of these 
squares, Such an arrangement will be called an m-ple Lattice. We shall further 
suppose that the whole arrangement is replicated n times. Thus v = k?, 


1 
b —nmk, r=m'n and in C (Cy); Ca — rn ( 1 7 } and c, when s Æ s' 


is either mjk or o aecording as the treatments s and s' occur together in 
a block or not. Let Sy (+, ), Sa (Ts), S; (<,) denote respectively the sums of 
the effects of the treatments (including the s th) occurring in a row, in a 
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column and having the same letter in the (i—2) th Latin Square as the s th, 
$— 8. 4,...,m. The same excluding z, will be represented by S', (Ta). S'a (a) 
and S'i (7). i =8, 4,..., m. 


Then the normal equations utilising intrablock information only are 


@=r [ 1- q) è- em 338408) (8.7.1) 
i=l 
= r$, — (n[k) 3} Si) r (8.7.2) 
i=] 
S; (Q) =n (m — 1) S;(2,) (8.7.8) 


From (8.7.1) and (8.7.2) we have 


m 

= Í Q, + $15,(Q,)/E(m — d jr (8.7.4) 
i=l 

Thus Vi & — Ty) — 2c? (k + 1)/kr or 262 [k + m/m — 1]/kr according as treat- 

ments s and s' occur or do not occur together in a block. 

The following normal equations combining inter and intrablock information 

can be obtained from (8.7.1) by following the method of Ex. 17 of chapter 2. 


P, =r [w + w'/ (k — 1] = aS ees s SiR) — (87.5) 


i-1 


= r$, — (w —w)n $18, (2, )/k (8.7.6) 
i=1 


whence S; (P,) = n (m — 1) w + w') S; (a) (8.7.7) 


T= [^ + (w — w’) S S, (P,)/k (m —1 w + «) | [ro (8.7.8) 
i=1 
From (8.7.8), V(t, — 2,,) is 
(2]r«o) {1 + (m — 1) (t — to')/E(m — 1 w + w')} or 


(2/rw) {1 + m(w —w')/k (m — 1) w +w')} according as the treatments s and s' 
do or do not occur together in a block. 


8. Partially Balanced Incomplete Block Designs (PBIBD). An incomplete 
block design is said to be partially balanced when 

(i) There are v treatments arranged in b blocks each containing k plots with 
different treatments assigned to each ; 
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(8) each treatment occurs in r blocks ; 


( iii) with respect to any treatment the remaining can be divided into m 
groups containing n, na, ..., Tim treatments, such that the treatments of the 
$ th group occurs with the given treatment ^; times. The treatments of the 
i th group are said to be t-associates of the given treatment. The numbers 


Ts fig... Tl; das Ag... Àm ATE independent of the treatment with which 
we start. Some of the Ws may be equal. 


(tv) if the treatment A is i-associate of B, then treatment B is i-associate 

of A. If A and B are t-associates, then the number of treatments common 

. to the j-associates of 4 and J'-associates of B is Djj, and is independent of 
the pair of treatments we start with, Also Dji = phy. 


Let N be the incidence matrix of the design. 
kb = vt 


Then (8.5.1.) is valid and we have 


(8.8.1) 
From condition (iii), we have 


Ma kna +... +n —p9—] (8.8.2) 


o 


The matrix NN’ consists ofr in the principal diagonal and n, times A, n, 


times Aa ..., nm times Am in each row and column. Th 


us 
NN’ E,, = L yo] Ea =N (NEn) = (rb) E, 
i-1 
S din A=r(k—1 ) (8.8.8) 
i=1 


The parameters b, o, r, k, Aq Aas s.s Arms Nis Nas ..., n, are known as primary 
parameters. From (3.8.1), (8.8.2) and (9.8.8), it is evident that there are 
2m + 1 independent primary parameters, Á 


is contained in the ^t; variates which 


4 — 1 t-associates of A, there ‘are 
exactly pi; varieties which are j-associates of B. 


Therefore, $7 pb =n — L If ij, 


among the n; varieties which are 
jei 


j-associates of A, there are exactly. Di varieties Which are at the same — 
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a m 
time j'-associates of B. Thus > Pi =n; Thus 
> HE 
m 
Sin —n—1lifj-i 
ml 
and = n; ifj7 #4 (8.8.5) 


Consider the group consisting of the i th associates of a variety and the 
group consisting of the 7 th associates of the same variety. The number. 


` of varieties common between j'-associates of variety of the Ist group and 


the varieties of the second group is Dire Hence the number of j'-associates 
that can be formed by taking a variety of the group and a variety of the 
second group is on the one hand n; pi, and on the other nj Phr. 


Hence 
4 nip = nj Diy, (8.8.6) 


Thus the m symmetric matrices 


Pi =[Pi], Pa = [r$ ] | -— E, = [pF] (8.8.7) 
A m (m?—1) . : À : 
contain only LETS M independent. elements. Consider the j’-associate 


classes of the n; i-associates of a treatment 4. How many times will a 
particular treatment B belonging to the j-associate class of A occur in the 
above set? How many times will A itself occur? Let B occur z times and 
for the sake of definiteness, let it occur in the j'-associate classes of 4i, 
Aj, ,..., dis Then Ai, 4i, ..., Ai, will occur in the j'-associate class of B 
and hence z — number of common varieties between the j'-associate class 
of B and the i-associate class of A i.e. Dh. In the j'-associate classes of 
4i, Ai, e., Ai, A cannot occur unless ji — i, in which case it will 
occur n; times. Let $;(7,) denote summation over the i-associates of 
Ta Then the aboye result is equivalent to 


m. 
Sy =S: (S) = Sing, Si )iti zi 
j=i 


m 
Ad SD = pi, S; (4) ifj =i (8.8.8) 


j=1 
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(i) Analysis without recovery of interblock information. Equation (2.1.29) 
can be written in this case as 


1 1 y 
Q =r 1— y er DS) 
cm 


"S (Q)mr(1— =) Sy O0 — E HS SAN) (34.9) 
i=l 


Remembering (3.7.8), and imposing the condition 2, + S, (2,) +... + 
Sm (3) =0 S 


k$;(Q,) = ai Sy (Fe) + ai So (Fy) Hooe H lim Sm (2, ) 
i=1,2,. m (3.8.10) 
m 
where aj = Nn; — Savi »t Aj 
l=1 
m 
aj —r(k —1) E Xin; — S app 5j =1,2,...,m. ° (8.8.11) 
l=1 


From (3.8.10), we can write 


S;(4,) = fa $,(Q,) + fia $2 (Qs ) aE sess + fim Sin (Qs ) 


i=1,1,..., m. (8.8.12) 
and finally à 
r(k—1)s —kQ, + SS R) SN fy (8.7.18) 
j=1 i=1 


Hence if the +, — Ty, is estimable and if s and s'-th treatments are i-associates, 
the best estimate of +, — Tẹ, is 


2> 


a 


^ 
aT (8.8.14) 


202 = E 
THES) [ k— Sin ^] (8.8.15) 


121 
(ii) Analysis with recovery of interblock information, 


with variance 


Following the procedure in Ex. 17 chapter 2 we have 


= 'Q'-— TEN d PI wi) NE 
Eee, TAK son [e+ Pi ] Wa cor STOOD 


` l=1 
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* w' k—1 R w— w’ = 
SAP) =r E bu ] gc SiG) p SAAMSOSIED 
1-1 


= diz $1(5) + ba Sa (Wa) + +--+ im Sm (5) 


7} m 
where b; =>, [ purus + 1; i; (w — :')— (t — w’) Zazi | 


tel 


whence $;(2,) = S Si 5; (Ps) 
: pn 
$21,2,...,m 


ALA kP, w—w' 
aa rer + Tegar AaS a S 


1-1 j-1 


Hence if the s and s'-th treatments are i-associates, the variance of the 
difference of their best estimates is e 


2 " 
Fee Dewy |*- 9-2 


1-1 
EXERCISES TO CHAPTER III 


1. Ina randomised block experiment with v varieties and r replications, 
variety 2 in Block 1 is interchanged with variety 1 in Block 2. Indicate the 
structure of analysis of the modified design and give estimates and variances 
of different treatment comparisons. What is the loss in efficiency due to this 
interchange ? 

2. In a' randomised block experiment originally planned with v varieties 
and r replications, it is found later on that there is not enough material of 
variety 1 and excess of material of variety 2. Variety 2 is then applied twice 
to blocks 1, 2, ..., 7; and once to the remaining blocks. Variety 1 is applied 
only in blocks (7, +1) (Ta +2), oT. Write down the C matrix of the 
design and obtain variances of different treatment comparisons.‘ Indicate 
how you will carry out the analysis of variance. 

8. In a randomised block experiment, variety 1 was used twice by mistake 
in the first two blocks and as a result variety 2 did not occur in these two 
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blocks. How will you analyse such an experiment and test the overall differ- 
€nces among varieties and also the difference between individual pairs? 

4. Inan s x s Latin square, two treatments in one column get interchanged 
by accident. Obtain estimates of different treatment comparisons and their 
variances. 

5. There are v new varieties to be tested. An experimenter divides his field 
into v blocks and each block into tv +1 plots. Each new variety is replicated t 
times in each block and a standard variety is assigned to the (tv +1) th 
plot of the block. Give the analysis appropriate to this design and compare 
it with the arrangement in tv + 1 randomised blocks of v plots each. 

6. Let N be the incidence matrix of a symmetrical BIBD. Consider the 
matrix 


ics TET V=) ci 

MN En N 

Show that MM’ = M'M = (r —3) I,,, and hence NN’ = N'N. 
7. Let.N be the incidence matrix of a BIBD, 


(i) Show that det N'N = 0 when the BIBD is non-symmetrical 


(H) Show that the eigen values of NN’ are rk and r—23 with 
multiplicities 1 and v — 1 respectively 


(iii) Obtain a non-singular G such that 
G (NN’)G = diag (rk, r —3, r — Aye, T —2). 

8. Denoting by N the incidence matrix of a BIBD with parameters v,b,r, k 
and A show that (i) Ej, N'N — rk Ej, (4) (N'N) (N'N) —(r—2)N'N-- 
k?n Ej Use these results to prove that if a BIBD can be constructed with 
v= $h(k +1), b — 3k +1) (k-- 2), k, r=k 4+2, x= 2, then between any 
two blocks there are either 1 or 2 varieties in common. Show further that 
there are 2% blocks having 1 variety in common with a given block and 
kik — 1)/2 blocks having 2 varieties in common with a given block, 

9. If N, is obtained from N, the incidence matrix of a BIBD by replacing 
each zero by — 1, show that N, N', = A(r — NI, + [b — 4 (r —23)] Ex. 


Hence N; NG =b I, if b = 4 (r —2) i.e. 2k =v +vv. < 


If N, and Ng be matrices obtained from the i 
possible BIBD satisfying these relations, show that the kronecker product 
N, X N; will lead to another BIBD. Obtain the parameters of this new BIBD. | 
10. Prove that a symmetrical BIBD can always be construeted with para- 
meters 9 = b = 2^ — 1, r = k = 27-1, } — oma Utilise this for the construc- 
tion of a new BIBD with parameters o — gm-1 _ 1, b=2 (gn-1—3) 
jih 2 eE and Na 9m-s. 4 


neidence matrices of two such 


2 
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11. Show that the Lattice Design of Section 7 is a special case of a PBIBD 
with two associate classes. Obtain its primary and secondary parameters. 

12. There are m k treatments arranged in m sets of k each and treatments 
of a set are assigned to a block and there are r replications. Show that the 
resulting design is a PBIBD with two associate classes. Obtain the para- 
meters of the design and show that the design is not connected. Demonstrate 
also that by choosing r properly, we can make the number of blocks equal 
to, less than or greater than v, the number of treatments. 

18. An incomplete block design is called a Linked Block design if (i) each 
block has the same number of treatments k (ii) each treatment occurs in r 
blocks (iii) any two blocks have the same number of treatments A* in 
common. Show that by interchanging the roles of blocks and treatments in 
a BIBD, one can obtain a Linked Block Design. Obtain the matrix D given by 
(2.1.15) for this design and show that that its rank is (b — 1). Show further 
that the eigen values of D are o and 1/r [k (r —1) — A*] having multipli- 
cities 1 and (b — 1) respectively. 

14. By taking v(v — 1)/2 combinations of v treatments taken two at a time 
we get a BIBD with parameters b — v(v — 1)/2, k =2,r =v —1 and a = 1. 
Taking blocks for treatments and treatments for blocks in this design, obtain 
a new design in which any two blocks have a single treatment as a com- 
mon link. Show that the new design is partially balanced and obtain the 
parameters of the design. : 

15. Show that in the case of PBIBD, the eigen values of NN' are r and 
the eigen values of A with appropriate multiplicities where A is the matrix 


m 


(aj) and aj = x Api; —n Nj xd 
Um 

m 
ag =r + 2 a pi — nihi, 0j =1, 2, ..., m. If NN’ has a zero eigen value 

Im 
with EX A u show that b > v— u and when the design is resolvable 
b>0+(r—1)—u. In the case of PBIBD with two associate classes, deter- 
mine the values and multiplicities of the eigen values of NN’ in terms of the 


parameters of the design. 
16. Show that in a PBIBD with m associate classes, the determinant | ly| 
m 


Klenk 


/ when k > r has the value zero where ly = r èy + 23 Me Pis — Aj ny and 
^. 


35 —1 if i.—j and 0 ifi Aj. 
17. Show that the efficiency factor (ratio of the variance of a varietal com- 
jawi in a design in ordinary randomised bloeks to the average variance of 
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the particular design occupying the same number of plots and having the 
same error variance per plot) for a p X q lattice is (pg — q)/(P 4 +p 4- q—8). 

18. Obtain the efficiency of the Lattice design in section 7 when analysis 
is made (i) utilising intrablock information only and (ii) with recovery of 
interblock information. 


19. p Xp xp varieties are arranged in a three way table i.e. Spatially in 
the form of a cube and treatments along lines parallel to the edges of the 
cube are then assigned to blocks. Obtain the variances of different treat- 
ment comparisons. Calculate the efficiency of the design relative to the 
corresponding randomised block design. 

20. Obtain three 4 x 4 orthogonal Latin Squares to construct (i) a BIBD 
with parameters v =b —91, r=k=5,2 = 1, (ii) a set of 15 orthogonal 
comparisons from 16 quantities a, V»... Tig (Hi) a balanced Lattice 
design with 16 varieties and (iv) a PBIBD with v 7:16, 0 —10, r=4, 
k=4, 2, =1, n, —12, A» = 0, n = 3, pli —8. 
` 21. A BIBD with parameters v, b, v, k, A is said to be affine resolvable 
if the blocks can be separated into r sets each forming a complete repli- 
cation and further ifany two blocks of different sets have the same number 
of treatments in common. Show that for such a design (i) the parameters 
can be expressed in terms of two integers n and £ (n 22,07 0) as. 
v=nk =n? [(»—1)t--1] b=nr=n (1? -Ent--1), A 5 nt E 1 (ii) the 
number of treatments common to any two blocks of different sets is 
k?|o = (n —1) t+ 1. 

22. An incomplete block design with v treatments each replicated 7 times 
in blocks of size k (< v) is said to be group divisible if the treatments can be 
divided into c groups each of size d so that any two treatments belonging 
to the same group occur together in A, blocks whereas those coming from 
different groups occur together in As (> 0) blocks. Show that (i)a necessary 
condition for the existence of such a design is rit —Az,0>0 and (ii) the 
following relations hold between the parameters: v — cd, bk =vr, 

r(k —1) = (d —1) X, + d (e —1) dg, T> oT mM. 
23. Denote the design of Ex. 21 by A(n, t) and the following Group 
Divisible Design by G(n, t) : 
v=b=(n—1) (nt +n +1), r =k =n [(n—1)t +1]. 
47 =0, às =(n—1)t+1 where c — n*t 444 
and d =n —1, n and t being integers, n > 2, 1— 0. 
Show that A(n, t) and G(n,t) either both exist or both do not exist. 

24. Defining Efficiency as in Ex, 17, obtain an expression for it in the case 
of an incomplete block design and obtain a criterion so that an incomplete 
block design may have maximum efficiency. 
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Establish that in the case of incomplete block designs having an equal 
number of plots per block, the BIBD is most efficient and in the case of two- 
way classified incomplete blocks of equal block size, the Youden Square is the 
most efficient. 

25. To each block of a BIBD with parameters v, b, r, k and Aa control 
treatment is added to each block, so that the block size is now (Je + 1). Give 
the analysis for the design, estimates of different varietal comparisons and 


their variances. 


CHAPTER IV 
STANDARD DESIGNS—Continued. 


1. Factorial Designs. Often instead of studying factors singly, it is more 
advantageous to study several factors in all combinations. Let the treatments 
consist of combinations of factors F}, Fo, ..., Fp, the ith factor being avail- 
able at Oth, 1st, 2nd, ..., (s; — 1)th levels. Thus à typieal treatment combi- 
nation in which the ith factor is at z;th level (i — 1, 2, ..., m) may be denoted 


by TUM ufi or by « (23, 25, ..., Vm). A linear function of the treatments 


al Wa, css Wm) T( Wy, Los noes Mm) (4-1-1) 
will be called a contrast if all c’s are not zero and 


Delay Rus 44) = 0 (4-1-2) 


the summation extending over all values of a, 25, ..., a, 


A treatment contrast will be,said to belong to the (k—1) th order 
interaction between the facto?s T.F ; Fi, if (i) C(x, Eb overs 5.) 
depends only on Zip ci, ..., Ci, ind (ii) the sum of C( aj, Za, ..., m, ) over 
any one of the EET Vis, Lig »-., Vi, iS zero. The 0-th order Uo OU 
are called Main Effects. If we write ra —ficTfiT.,fü, i—1,9,..., m, 
then a typical linear function of treatments belonging to (k —1)— th. TEES 


interaction between the factors Fi, TOR A Ey is 
$1 po. " 
DU TEUER Vh, Us o) ds d, fi, ia d i (4:1:8) 


where 0 2 (ain lin een zi) = SM Go ins ry Oi) = 
gi, 


EI Tiga ans Di. ) (4:1:4) 
One a ious way of obtaining these linear functions is to choose 
l(mi, Wi ea Ti) =h m la wig. Ik DA 
where L. is such that Sa D =0 
J i J Ti 


5 
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j =1, 2, ..., k. Thus there are (= HOY = 1h) eos (si, =r) 
independent contrasts 


$1 ba -hm l 2 1 a | l 44 
oot ESS A 253 . fäi "S Ress afl 
$i, Di, $i, m lei, i H" = 2 ti, fi 2 em k vi, m (4-1-5) 


belonging to the ( k — 1) =th order interaction between Fi, Fiy ..., Fi. 
From the definition, it is evident that (7) interactions of different orders are 
mutually orthogonal (i?) two interactions of the same order when the 
interacting factors differ at least in one element, are mutually orthogonal. 
Thus the total number of independent treatment contrasts sy. 9. ...... Sn — 1 
can be broken up into (1)s,—1, Sa — 1, ...... s $,— 1 sets of independent 
contrasts belonging to Main Effects (2) (s, — 1) (s —1), (s, — 1) (83 — 1), 
o (8m-1—1) (55 — 1) sets of independent contrasts belonging to lst order 
interactions, ...... (m) ($ — 1) (sa —1) ...... (54 — 1) independent contrasts 
belonging to (m — 1) —th order interaction. Note also that s, s,......: Sin — 1 


=[1+(%—1)] [1+(s2—-1)] - [14-(54—1)] — 1 


m^ 
-9s-1)4 $i (-16 -1)- + (=D) ( 71) (8,1). 
EN ET 1) 
Let the factorial experiment be carried out in r randomised blocks of 
($1 Sae- Sm) plots each. Then the best estimate of a treatment contrast is 
the corresponding lincar function of the mean yields. Let Tie in 
(2i, hy aoo ai) denote the total of treatments in which the 7, th factor is 


at a; —th level, i, th factor is at the aj, —th level, ... and ij the factor 
is at wa— th level. Then the sum of squares due to main effects of F, is 


$1—1 2 T? 
c A 
S Deia — ———__,, T denoting 
T S3 Sg... T $1 Sa -Sm 
%=0 


the grand total. The sum of squares due to interaction between F} and F, is 
4 = 


2 
9 S Tyo (ty; a) — Sum of squares due to 
dua, (fae) 


— Sum of squares due to main effects of F, 
T? 


T Sy Sam 


main effects of F; 


The caleulation of sum of squares due to higher order interactions is also 
obvious: When s, —54— “= sm = (say), the factorial design is said to 
f p 0$ e 


1 1 r SEGO - m 
be symmetrical and often it will be expressed by s". 
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2. Split-Plot Experiments. Here for experimental convenience, the plots 
of a block containing a main treatment are split up into a number of sub- 
plots to accommodate an equal number of subplot treatments. ‘The subplot 
yields on account of having the same main treatment common and on account 
of the contiguity of the plots, are expected to be correlated. Let yir denote 
the yield of the plot containing the k th sub treatment of the j th main 
treatment in the? th block (i =1, 2, ...,7; j —1,2,...2; k= 1, 2,...6) Then 


Vir SE + ait pj Tr + Sje t Egk (4.2.1) 
Where a; denotes the effect of the i th block, pj the effect of jth main treat- 


ment, 7; the effect of k th subtreatment and 8 ,,, the interaction between j th 
main treatment and k th subtreatment and tjj S are normal variates. 


E(Ej,)—0 
V (ejr) = V ( Yür) =a 


Cov (gio Yer ) =P 07 if i=’, j —j', k Fk’, 0 otherwise. We also assume 


E a; = Yo; = Uy, —X54 = Edy = 0 (4.2.2) 
i j k j k 
Let 
Vir EDS. bra Mba; big—1 U; 
v8 4 
Yijz ES NL MT —1 Zj 
= Ve 22 921 aß p (4.2.8) 
i Ls 
v bg: bee .. beg—1 Zij (g—1) 


Vijg ve 


where the first matrix on the right hand side of (4.2.8) is orthogonal, Tt is 
easy to verify s 


V(U4) = o} = o° E ve | = VEA Ts ac uad 
w 


and all covariances are zero; i= 1, 2, ..., r; 271,2, ... à; 


A=1,2, ..,B —1. The normal equations arè obtained by equating the 
differential coefficients of 
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T — r a 2 
oS S S | zi | to S [ 5-50] 
i = i=l j=1 
r X B-1 B ^ 2 


$5 [ za- zi he (e+ du) | 3 


l A=l1 k=1 


Uum ži Us — [VB (u + a + 225 


with respect to the parameters to zero. They are 
w ray... = w [rapu] 


w apy.. =w [app Hapa] 
io Saye 


Re eee DE 


^ j—1,2,. 
w[ra(y oj 
HT w'[ray...] Wes + wran 


) 
k=1, D BN 


w[r(yjg— Yj) Hw [ry4.] = wr(z 4-3 5.) +w rty) 
OR eos ra d SERES l (4.2.4) 


“From ( 4.2.4 ) we easily get the results in Table 4.2.1. below. 


d TABLE: 4.2.1 


Best Linear Estimates of Some Linear Function of Parameters with Variances 


Linear Function of Best Linear Estimate Variance of best 


Parameters Linear Estimate 
ai — Oy Vu. —W.. 203/08 
Pj — Py y. Yy 203/rB 
Tp ^7 YY 2ci|ra. 
Sj, — Oye Yoja Yoja — (Yj. — Wy) = 1— à) 
Sin — Sie y Ve — (Y. y.-k’) ex 12) 


64 DESIGN AND ANALYSIS OF EXPERIMENTS 


5s ^ ^ ^ 

Lun Ceu EU EUST a Ys Land 

$2 =Yik — Mj — Y- y. isa set of solutions of the normal equations 
(4.2.4). Hence 


Regression sum of squares 


=w" [rag y*...-Fag Zi (y; — y. i *rBZ(ys —y...)? wiran (yny) 
TPES (Hi — Yj. — Yr ty)? ] 
Error sum of squares 


=w BED (yy. — Yi, — yy, ty)? + WEEE (Yi — Yi, = Yin yg.) 


Optimum tests are not available if w'/w or pis unknown. In the absence of 


knowledge of 9, we can, however, perform certain exact tests of significance, 
For if we set 


1 
NEUE yaar ist S Din Cin (4.2.5) 
A=] 


where the b’s are the same as in (4.2.8 ) we have V ( Cia) = Lw, V (ny) 
= 1/w', all covariances and expectations are zero. 


Now 


BBS a S tent 
Sa 2 25 (Cin — 6.5 y (4:2:6) 
DM REA 


T [A 
= Oy Waren. Trang =a)? EEE (ny — ne — ny tna)? 
S " a 


By Cochran’s rule, the terms on the tight hand side of (4.2.6 
independently as Z? o? with &—1, (x—1)(8—1 Janda (7-1) (B—1) degrees 
of freedom respectively. Similarly the terms on the right hand side of 
(4.2.7) are distributed independently as X e$ with df. Wen al 
(%—1) (r —1) respectively. Now 
TE (Hy ne)? =rBE(ej.—e...)2 = TEE (y.;.— 
Na sca = Vas a = Ng SO. 
ale —2...)* = «8X (ej... — e.. Y= aB (yr =y), 


if a =o =... ==, 


) are distributed 
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Z2(ny— ni =h +...) 2? = BZZ(ej.—s,.. LESSE)? 
= BEE (yy. — Yi — yy. y... )? 


B—1 B 
TX s Con ES S (s... — &...)* 
A=] k=1 
=rad(y...~ —y..), if t1 = Ta =...78 =0. 


TABEL =rdy (ej, —e;.—e..., +...) 
=7 ll (Yoja — Yo — Yr HY) if 85, —0, j=l, 2, 355/665 
k=l, 2, ..., B. 
DEX (Cy, — ja)? = VUE (ey, — ey. — Eje + ey.) 
—EXX(Wug Ug. — Yoja yu.) 
Hence we have the following analysis of Variance Table. 


TABLE 4.2.2 


Analysis of Variance of Split-Plot Design with Main Treatments, 
Sub-treatments in r Randomised Blocks. 


Dn — 


Source DF SS 
Blocks T —1l «pZ(gy..—y...)* 
Main Treatments x —1 TBZE(yj.—y...) 
Error (a) : 


Main Treatments x Blocks (x—1) (r—1) BZZ(yj.— y;...— Yg y.) 
a c 
Total TX—1 BEE (Yj — gy...) 


Subtreatments B=1 Tan( ya — Youn)? 


Interaction : 
Main Treatment x Sub- 


treatment («—1) (B—1) TERU Ya yy. — y by os 2 
Error (b) : «(r—1) (B—1) ZZZ(yg. — yy. — y. +y)? 
Total «r(B—1) EEE (Yj — Yy)? 


—————— aaa aaa aa — —À 


` Total raB—1 EEE (Yi — yY)? 


poe 
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For testing Block S S and Main Treatment S S, we have to use Error (a) 


and for testing subtreatments and Interaction between Subtreatments'and 
Main Treatments we have to use Error (b). 


Let us try to compare the efficiency of the split-plot design with a rando- 
mised block design having af plots per block. Writing Tjk Y; +t + Sz, 
we can obtain the following table of expectations : 


TABLE: 4.2.3. 


Expectations of sums of squares in Split-Plot Design and the Corresponding 
Expectations had the experiment been designed in r Randomised Blocks of «. 
plots each. 


Split-plot Design Randomised Block Design 
Source E(SS) Source E(SS) 
Main Treatments rg X (qj. —7..)? + 
(a—1)o3 
Error (a) (r—1)(«—1)c2 
Sub treatments ra(t.,.—7..)?+ Treatments r XX (zig — 7..)? + 
(B—1)ef («B — 1) oig 
Error (b) &(r—1)(B—1)oi Error («8—1) (r — 1) ofp 
Total  rEZ(vg—T..)*--r(«—1)o2 Total 


) rEXG—gh 
+ ar(B—1) oy 7 («8 —1)odg 
———————————— ka 


c?  (a—1)c2-d &(B—1)o2 


. = = —1 
a cxoWaBrerj te (eee) 


So the subtreatments and interactions between main treatment and subtreat- 
ment are measured with greater precision whereas the main treatments 


are measured with less precision and the overall precision is the same in 
both the cases. 


- M 
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3. Split-Plot in a Latin Square. Here the main treatments are arranged 
in a Latin Square and the plots are further split to accommodate an equal 
number of subplot treatments. Let Vir denote the yield of the sub-plot 
in the jj the plot having k th subtreatment, the main ‘treatment being l, i, j, 
1—1,2, ..., t; k = 1,2, ..., r. We can write 


Yih — UT 65 + B; + er + Sr T eg 
where E (ej, ) = 0, V(eiy) = o? and 


Cov (egp Ene) — pc? if i— i', j =), kk’ and zero otherwise.. 


Let ny) = [sii + Eja +... TF Ej — Sial VAA FI) 
G7] eb A NGING h coy oca 


Ug = Nie = (Cg tuas Heyr) | Vr sij =1, 2, unt 


Ther, expected values of the new variables are zero and covariance between 
any two of them is zero. Also 


Veni) Te e) =0f, V() =0 [1 +(r—1) 9] =o}. 


Now 


1 Li Lx t 
I E xXowp—Cuht.clit X (u —u..) +t X (uj — u..)* 
i-i j=1 


i=l j=l 


t t t 
TÉ EX(uwg—u.)?-zE E o(ug—us—u;— tF 2u. . )* 
te i=l j= 


T—l 
II S 2 ` Th =e 2x Watt E S Cna — 7...) 
1 AEST ZEN 
zb- » Si Si [5910] $ 
i J AN 


In I, by Cochran’s rule, the constituents on the right hand side are distri- 
buted as x? o3 with d.f. 1, £—1, t—1, i—1 and (#—1) (t—2) respectively. 
In II, again by application of Cochran’s rule, the constituents on the right 
hand side are distributed independently as y? of with d.f. r—1, (t—1) (r—1) 
andt (t—1) (r—1) respectively. Hence we have the following analysis 
of Variance Table. 
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TABLE 4.3.1 


Analysis of Variance for Split Plot int x t Latin Square 


———————————— 


Source DF S. S 
t 
Rows 1—1 : rd (Yi. —g...)2 
i=l 
Columns t—1 rd (yy.—y...)? 
j=1 
t 
Main Treatments  £—1 rd (Yem —y...)? 
1-1 
Error (a) (—1)(0—2) © XE (y p gg 3 
Total ?—1 TED (yy —y...)? 
F, T 
Subtreatments r—1 Bx (y.— 9) 
kel 
Interaction 
Main Treatment x 


Subtreatment (r—1)(#—1) PIE Yie — Ya ay — y. tyne)? 


Error (b) t(t—1)(r—1) - BEE (Yit — y — Yi FY.) 


4. Repeated Subdivision. In section 2 
into y sub-sub-plots to accommodate an e 
our set-up is 


; if the subplots are further divided 
qual number of sub-sub treatments, 


Yia =U + B; ty, + 48 jk + Orr + dg 
ee 


vino = Pi sedia Sin San Su Ba, 


kel 


-Ee-Es-Eu-fe-Ne- ume Situ 
E (egu) =0, V (egu) =0? 


Cov ( tip epp ) = 01 o? if i = Cek k lY 


and Sp ifii yk 


and — 0 otherwise 
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Let us consider the three sets of variates 
I. ning? =L 2, +9139 =1,2,...,0; k—1,2,..., B; A= 1, 2, ay] 
IL viy T= 1, 2, Goo ELSE ULP cete) be ec [e es 1l : 


mu. 15 gg BTE E haves OLS 


ijo 
Nijk1 03 Qor ++. AY Ek Eijk 
Nijk2 Ai Age «+» AYa Eijk? A Eijk2 
NijkÝ —1 QY—1---4YY_1 tükY EijkY 


[477 Oy ba bgi VY DT VA egi 
Cue d big Doe bga Vy Eija =B Vy Eija 
gi bip—1 02 g—1 --bpp—1 VY Sap. Vi sup. 


where Ei J s] is orthogonal 
al vy = BY Sij. j 
The new variables are all independent having expectations zero and 
variances of the different e are respectively Xj—o? (1—p, ), Z2 = 03(1—p’), 
= o {1 +(B-1)e 4); 


Y es 
l-c(y—1l)e 


where cå =o? {1 EE EE r= 


Also we have 


[A Y-1 


2 g 
Bop a in aren Sh. AARIN (5-4 = -E 
AS Se ees 

HBS ju MASS eA Su d 
RED TC. 

*221512i n — an 


n ggg Osor i t+ gi Si tamt 
i j y 


y=1 


2 2; 2, 2i (čij — Cy )? 
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III EXg—rav..-EaXE(v. —v.. E rZ (v —v..)* 
LE] t 2] 


+22 (05 —0;.—v.5+0..)? 
ij 


Straightforward application of Cochran's rule 


gives the following analysis 
of Variance Table. 


TABLE 4.4.1. 


Analysis of variance of Split Plot Design with r splits of subtreatments 


ee 


Source DFE SS 


Blocks r—1 aby X (yi... — y...) 


Main treatments « —1 TBYZ(y...—py....)2 
i 


Error (a) (r—1) («—1) BYZE (yy. yy... — s. . Ey... 
ij 


= : E 


ubtreatments 6-1 TAYE (Y. e =Y... JO 
k 
Interaction : 


Subtreatment x. 
Main treatment (&—1) (8—1) TYEE (Yik. —y.,. Y. e HY.) 
ik 


Error (b) &(r—1) (B—1) YEER (Yije — V. -Yje yg)? 
ijk 


Sub-sub treatment y—1 T«pX (y... — yore)? 
1 
Main treatment x (x —1) (Y —1) 


LESSE yen ae 
sub-sub treatment ji (Vi y 1 Yj. ty.) 
Sub-treatment x (B— 1) (y—1) BION Lum á : 
Sub-sub treatment ki Wa Ya e Yr ty.) 


Main treatment x (a—1) (8—1)(y—1) 


Sub-treatment TEENS — yn) ay 
ub-treatment x 

V5. d- y... Me MAYA 
sub-sub-treatment ja TY EY. pa — quus) 
Error (c) €B(r—1) (r—1) DUIS 


(ga — yg — Yir. — Yje) 


n NG ee 
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Y—i 
V (Y.a y.u) =V [ zx (an — 401.4 | 
2e 


ADEDE] 


pl 
WOY. aaa Ya) =V E (biy—bry) Y. -y ] =2 33 /ray 


v=1 


Similarly V(y.j.. —y.j-.) =232/aBy, V(y--m—y- om) 22Z2]ra 


PSI 
V(yji— Yz) =V Ls Qj (Mjaa | = 2i = Te CBE 
NET TBY 


y 
V(yji—Vjr) = TE (5. -annaa | =22;/rß 
V(Y. r> Yeri) — VI. — Y- -r 4 X40, (A. — naka) 
o _ 2X; d 227 (¥—1) 

ray ray 


^ 5, Split-Plot with sub-units in Strips. Here the subtreatments are arranged 
in strips across each replication. Let y,;, denote the yield of the plot having 
the k th subtreatment, j th main treatment in the? th replication (i = 1, 2,...7; 
(RES Riot nd 5 ems Pss B). Then the yields yj;, yyy, and gi, > gy 
are correlated. Our set-up is : 


Yije =H + t; + Y; + + Oe F Eije 
with La; —Zy; Te pe iQ 
i i k 


E (cgr) =0, V (ejr) = 05 Cov (eges ere) = p1 0° 
ifiai, j=j’, kk and =pa0° if i—i', k—k', j #j’ and =0 otherwise. 


Let (eg)  A(vBeg- to c ri) 


where A= 


1/8 agi a2 -agg 
is an orthogonal matrix. Then the two sets of variates 
L qu (i—L225j =L2. 4b—12.548—1) 


IL y= VE syo (FH LPF =l hes t) 
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are independent in the probability sense. "ija ' S have the same variance but 
Hy, and $54 are correlated. Similarly u; and Uy are correlated. Making use 
of a further transformation 


1/MVa bir 035... big — 1 Va NED Mania 
1 Boy Dog... bou — ti i 
ux) [Va bzr Ds... Da, 1 il, =B Ul 
1/MVa ba, baz... baa — 3 Gi(a—1). či (an) 
Va t 
where B is orthogonal and ( ug)—B vp 
Ui (x—1) 


We get the following sets of variates which are independent. 


I. vig = Va U; = Mah tj, DERI PAD 


a a 
IL. vi, = S bj, wj = VB > biy ej, $—1,2, T; y= 1,2. 


A — l. 
j-1 Sex 
Me sx Va ny em 3 or D Wa E oh 6 a 
k=l 
a o B 
IV. Gia => b jy Tij) = y 21 b jy aky Eijk 
Jæ j=1 k=1 
Kg Baka A lea a Gl 1,2,..., B—1 
Thus 
1 ihe) oe B—1 
ijk = ix + —— Doe 1 
Map 7% VB E jv ?iy + NE a [2527 
Ba gt A 
AÉ 2i aka bi Cun 


| 
| 
i 
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Now 
EXo?—rXwg-EEX(vn—v JB 
wv 
2 EEE? —rEES t EE(b: EQ» 
A QN 


EXXXUGQ-—rXàXX Cy? + DUD | CT C Y 


Straightforward application of Cochran's rule gives the following analysis 
of Variance : 


TABLE 4.5.1 


Analysis of Variance of Split-Plot Design with sub-unit treatment in strips 


Ln. —————— À— 


Source DF SS 
Main Treatments a—1 PREYS — y.-.)* 
Error (a) (r—1)(«—1) B XX(yg.— yi. .—J.4. + Y) 
Subtreatments Bp—1i r&X(gy..y — y...) 
Errox (b) (ri iun) LEE (yi Wi. «V. rty)? 


Interaction : 


° «Main Treatment x (x«—1)(8—1) TEN (yayi a Yy FY) 


= 


Sub-treatment 
Error (c) (r—1)(«—1)(8—1) DEE (Yi Yin Yj ja T V 
e HYS ct U- Youn) 
CTI MUS ULTRI eee 
We enn nlso establish very easily that 
V (gg. — Yj) E 201 [76 
VY — Yew) = 208 [ra 
V (9 = Vu) 2[(@—1)o§ 4+ oi ]/rp 
V (Yir — un) = 2[(8—1)o3 +05 ]/ra 


|| 


where o?, o$ and og can be calculated terms of o*, p, and ps. 
The analysis of sub-units in strips when main treatments are arranged in a 


Latin Square is exactly similar. 


6. Confounded Arrangement in Split-plot Designs. Suppose that we haye 
p main treatments and s" subtreatments. Then we can have ps whole plots 


and in a set of s whole plots having one particular main treatment, we can 
$7—! sub-plots and assign to them the subtreatments 


split each main plot into t : 
d.f. belonging to the highest order interaction of the 


in such a way that (s—1) 


74 DESIGN AND ANALYSIS OF EXPERIMENTS 


m Subtreatment factors are confounded with the whole plots. "The subdivi- 
Sion of the degrees of freedom is as follows : 


D.F. 
Between Blocks b —1 
Between Main Treatments p—1 
Highest order Interaction of Subtreatment factors s — 1 


Interaction : 


Main Treatment x Highest order interaction of 
Subtreatment factors (p —1) (s — 1) 


Error (a) (sp—1) (p — 1) 
Subtreatments s (s™—1 — 1) 


Main Treatment x Subtreatment s (s™-1_]) (p—1) 

Error (b) sp(b—1) (s™—1—1) 

Total bps™ — 1 e 
Another method would be to take s y blocks (of p whole plots) consisting of r- 
sets of s blocks. In any one set, the main plot is divided into sg"-1i subplots 
to accommodate 5"—1 of the s” subtreatments— different subsets of s”-1 of 
the s" subtreatments being used for different members Of a set. Here the 
subdivision of the s" subtreatments into 's sets of s™—1 is such that interest 


contrasts belong to the highest order interaction of the Subtreatment factors. 
The partitioning of the total degrees of freedom is as follows : 


Source D.F. 
Between blocks rs—1 
Main Treatment Wil 
Main Treatment x Highest order interaction 
of the subtreatment factors (p — 1) (s — 1) 


Bros (a) $(r—1) (p —a) 
Total between whole plots rsp — 1 
Subtreatments s (smi. 1) 


Main Treatment x Subtreatment s (571.1) (p—1) 
Error (b) 8p (r—1) (5713) 
Total m 


T$"5—1 
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EXERCISES TO CHAPTER 4 


= 1 1 =" Be aes: Ai Aia 
L tet = [1] [maa A “ile RE 


and X; — (X; 4, X; Xi}, i=2, 8, ..., m. Show that 


A; 


i 


xc Be Apa Xa 


(4—1) Ap; Xi, 


2 Hence, establish Yates' process for obtaining main effects and interactions 
from total yields of treatment combinations in a 2? design. Show that this 
can be used inversely to obtain the yields of treatment combinations from the 
total, main effects and interactions. 


2. Show that Yates' method can be extended to designs in which some or 
all factors are at more than two levels. 


| Extend in particular Yates’ method to 3" design, taking 
S9 


© 


n wei vii 
A, = —] 0 1] ;X,={1, z, a3} 
nce. 43 : 


A; = —AL a o Aya ; X={X_, a; X a2; Xi 4 b 
Ai —2A;, Aia #=2 8, ..., m and 
noting 
(14-2; + 27) Ay, X45 
A; X= (—1+27) Aii Kira 
Aia Xia 


How will you arrange the treatment yields in order to make Yates' method 
operative ? 
3. Deduce that in a 2” factorial design replicated r times in randomised 


blocks of 2" plots, each main effect and interaction has variance o?/r.2™—2 
| where o? is the per plot intrablock variance. 


i 
| z (1 — 2m; +2?) 
| 
| 
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4. Let 


$—1 


| KN ba ce 
442) 3s (2) X a (e) (41 (i021) X, s Goya) (i-a) 
be orthogonal polynomials for the set of values @ — 0, 1, 2,......... $; — 1. 
We ean define for ; = Zuo Mos 
5,—1 
Ea D ae Xr CO Sc M 
z=0 


as due to the linear, quadratic, cubic,.....effect of Fi. Show that the sum of 


squares due to the above effect in a randomised block experiment with 7 
replications is 


F (1) s (8:21)... (s 2—2) 
[zs Fu | T5555 (BOT (BE pI) I 


5. Three manures N, P, K each at two levels 
eight manures Ny) P, Ko =a, N, P, K, =>, N, P, K,—coN P,K3 — d, 
i4,K, =h. Give the 


6. Show that split-plot arrangements are Specialised types of designs 
involving the confounding of main effects, 

7. Show that when the main plot is Split into tw. 
can be substantially simplified by calculations of 
adjoining sub-plot yields. 


9 parts only, the analysis 
Sums and differences of 


8. Show that the main effects, 
in a 2" design are linear funetions 
(1), a, b, ab, c, ac, be,, abc, — 


two, three and more factor interactions 
of the treatments 


9. In the analysis of split plot design with « main treatments and f sub- 
treatments in r randomised. blocks, arrange the Observations Vas as Zy, Zo,- .Z 
where. A = (i — l)eg + (j — Ue + k: ijk EATA 
Write down the dispersion matrix Vand its inverse V—1i E 

- - Util " 
to obtain the normal equations. ise the results 
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10. In Ex. 9, obtain best unbiased linear estimates of main treatment 
contrasts, subtreatment contrasts and utilise the normal equations to obtain 


their variances and covariances. 


11. Obtain the expected values of sums of squares in Table 4.2.2 and 
Table 4.3.1. 


12. Suppose that there are no main plot treatments and there are p plots 
in an area. Divide each horizontally into b strips to accommodate b levels 
of treatment B. Independently of this divide the plots vartically into c 
strips to accommodate c levels of a treatment C. Write down the structure 
of analysis for this experiment, the partitioning of the total degrees of freedom 
and the variances of different treatment comparisons. 


18. Indicate the change in the analysis when a set of main plot treatments, 
4, say at a levels, are imposed on to the Experiment in Ex. 12. 


14. c, and c, are the standard errors associated with main-plot and sub- 
plot respectively, there being f sub-plots per main plot. Establish that the 
variance between any two treatment means not belonging to the same main- 
plot involves the weighted average of c; and c$, the weights being in the 


ratio 1: (B — 1). 


15. There are « main treatments where « is even, 2r blocks of « plots 
arranged in r sets of 2. Each main plot is split into 2"7-! parts. The 2” 
subtreatments are divided into two sets of 2"! such that the interest contrast 
belongs to the highest order interaction of the subtreatment factors. To 
«|2 main plots of the first block of a set, we assign the first set of subtreatments, 
and the second set of subtreatments to the remaining main plots. In the 
second block of the set, we assign the subtreatment sets in the reverse order. 
Give the computational procedure of the arrangement and the partitioning of 


the total degrees of freedom. 


16. Give the structure of the analysis of covariance for the split-plot design. 


Obtain the adjusted plot yields and the variances of different adjusted treat- 


ment comparisons. 


A8x2x8 design is applied to the whole plots in blocks of 6 whole 
partially confounded and a factor D at two levels is 
Give the computational procedure for the scheme 


17. 
plots, with AB and ABC 
applied to the subplots. 
that you would adopt. 
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18. Given a two-way whole-plot design with split-plot superimposed on 
one of the whole plot treatments. There are sı Varieties and they are assigned 
to the columns randomly ; in addition there are Sa fertiliser levels which are 
assigned to the rows randomly. sg rates of seeding are to be tested and they 
are split-plot divisions of the varieties. There are r replications, Obtain 
(2) the analysis of variance for the scheme (ii) the best unbiased linear estimates 
of defferent treatment comparisons together with their variances, 4 

19. In a split-plot experiment with r replications, a whole plots and B 
sub-plots, take the yield’ of the plot in the first replication having first main 
treatment and first sub-treatment to be zero. Let ?;, the value of a con- 
comitant variate, be —1 for this plot and zero elsewhere. 
of variance and covariance table and deduce the value of 
sion coefficient. 


Obtain the analysis 
the error (b) regres- 


20. Give the principle of construction and analysis of a design for a 4 x 22 
factorial where the units are in a 4 x 4 Latin Square and the plots can be 
divided into two halves only, smaller units being impractical, 9 


CHAPTER V 


APPLICATIONS OF GALOIS FIELDS AND FINITE 
GEOMETRY IN THE CONSTRUCTION OF DESIGNS 


1. Galois Fields. A field F is a set of more than one element for which 
there are defined operations of addition and multiplication which satisfy 


Commutative Law: a +b =b +a, ab = ba 
Associative Law: a+(b+c)=(a+b) 4- e a(bc) =(ab)c 
Distributive Law: a(b+c) = ab + ac. 


Laws of Inverse Existence : 
(i§ For every pair a, b there exists an æ such that œ + a = b 
(ii) For every pair a, b satisfying the condition a z& 0, there exists an y 


Se such that ya — b. 


A field containing.a finite number of elements s is called a Galois Field and 

is denoted by GF (s). A Galois Field can always be constructed when s = p” 

where p isa prime and n is a positive integer. When n = 1, the residue 

' elasses modulo p constitute GF (p). In the general case the residue classes 

modulo F (x) where F (a) is a ploynomial in æ of degree n which belongs to 

and is irreducible in GF (p), constitute GF (p"). These functions F (z) are 

known as minimum functions. Minimum functions likely to be of use in 
Design of Experiments are listed in Table 5.1.1 on page 80. 


The elements of GF ( p” ) can also be expressed by 0, 1, o, ad, ...o?" -? where 
0,1 and « denote respectively the null element, unit element and a primitive 
root of 2? -1 — 1 = 0. Thus the nine elements of GF (8%) may be identified 
as ug =0, €, =l, u5 — e, ti, =o — x +l, u, — 03 = 2x +41, uy = = 2, 
ug = «5 —2a,u, = «ê —2«--2 and u=a?=a+2. Now ts + us = 
2 -- 2x + 2 = 2« +1 =u, and ts. Ug = od. a = «9 = a? = tus. 


If F(a) = a, + ma +...+ a,” where ag; Qy; d5,..., a, are elements of GF(p) 
is a minimum function for GF(p") then it can be established that 


n 
—1 : S 
a! — (— 1)” ay mod F (v) where t = m . Thus knowing the series 0, 1, a, 


ahi = 2 up to of - 1 all others can be constructed by multiplication by 


ag and its powers. 
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80 


That Z—ap=,0 


Sg 
BHLG+ yo 


6—Ut—gU 


g— t=; 


6— NGU yi — gV = gv I+2=,v 
o+r7=,0 Gan 
Guha 1+ 2+0G+,2g= Qv 
Sg GUS v 


Graco 1 +ag=_v 


suogounyr wnunurpy fo ST 


TI'S rSv 


TUHU hetu 
Ihsan, 
Lra—,m 

I+o= gv 

Ie? ot 

I=? 

Lagu 


I+2=,2 


bz] 
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2. Finite Geometry. The symbols PG (k, p”) and EG (k, p”) stand respec- 
tively for finite k-dimensional projective geometry and finite k-dimensional 
Euclidean Geometry. 


In PG (k, p"), (o, 8i, Las ..., ay) where a, Lys Los ..., v; are elements of 
GF (p”) at least one of which is different from zero will stand for a point. 
If o 40 be an element of GF (p"), then (a, ay, 5, ..., %) and (pay, pay, ..-p25) 


E+ 
will stand for the same point. Thus there are a s ka 47 
DES 


points in our geometry, s being written for p". All the points which satisfy 
the (k — g) linearly independent equations. 


Qio Uy + Ay, 9, +... Ay % = 0, i= 1,2,...k — g... (5.2.1) 


will be said to form a g-dimensional subspace or a g-flat in PG (k,p”). 
Evidently a g-flat is generated by (g +1) points, the matrix of whose co- 
ordinates has rank (g +1). The number of points in a g-flat is evidently 
sh t1_y 
$—4 
It is now easy to find out the number of g-flats in PG (k, p^), g< k. A 
£.flat is generated by (g + 1) linearly independent points. The first point 
can be chosen in P, ways. The second point in P, — P, ways, the third 
point in Pj — P, ways and the (g +1) th point'in P, —P,., ways. Hence 
the total number of ways of choosing (g +1) independent points = 
P, (Pj — Po) (Py —P,)..(P, —P,..,) Each g-flat can be generated by any 
one of P, (P, —P))...(P, —P,_1),(g-+1) linearly independent points. 
Hence the number of g-flats is 
Bis Bie) «(Be PIE) 


= P, say. 


FH CERE CO TCSESS ONE OE UN 
(s** 1 — 1) (s* — 1)...(sF-4*1—1) 
C qgr*1—3) (FD)... — 1) (6.2.2) 
One can easily verify that 
Q (k, 5,8) = (sk f=") (5.2.8) 


An ordered set of k elements of GF (p") (v, a, ...2;,) will be called a point 
in EG(k, p"). Two points (a, a, ... z,) and (a, aa, ..., a) will be 
identical if and only if a; = e, í =1,2,..., k. Thus there are s* points in 
our geometry, s as usual standing for p". The set of points satisfying the 


consistent and independent k—g equations. 

Qio + dj; y+ litat ... + ay, %, =0, $—1,2,...k —g (5.2.4) 
will be said to constitute a g-flat of EG (k, p"). Evidently the number of 
points in a g-flat is s*. 
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From (5.2.1) and (5.2.4), it is obvious that EG (k, p") may be derived from 
PG (k, p^) if we exclude from PG (k, p") the points satisfying 7) = 0 which 
constitute a (k — 1)-flat in PG (k, p”). From (5.2.8) the number of g-flats in 
EG (k, p^) where g < k is thus 

> (558) —6(k—1s8g) (5.2.5) 


8. Orthogonal Latin Squares. In section 2 of Chapter 8 we have already 
defined a Latin Square. If two Latin Squares with the same letters be such 
that when the two squares are superimposed each letter of one square coin- 
cides exactly once with each letter of the other square, then they are said 
to be orthogonal. — — 


A set of Latin Squares, such that any two of them are mutually ortho- 
gonal will be said to constitute a set of orthogonal squares. If the Latin 
Squares be of side s, the number of squares in an orthogonal set can be at 
most (s— 1). To understand this, we first of all note that permutation of 
letters in the Latin Squares does not destroy the property of orthogonality. 
Thus we can re-write the squares of the set such that the first row in each 
is AB... Now between any two squares B coincides with itself in the first 
row, hence in the first column B cannot occupy any position twice or the first 
position. Thus there are at most (s—1) positions available for B. When 
a set of (s — 1) orthegonalised Latin Squares of side s exists, the set is «cid 
to be completely orthogonalised. 


When s = p" where p is a prime and n is a positive integer, a set of com- 
pletely orthogonalised Latin Squares can always be constructed. 


Let ug—0, t =l, U= &, ug — o3, ..., u, 1 = «^ ? be the elements of 
GF(s) Then in the z th row and y th column of the i th Latin Square, we 
put the subseript of the element 

U; Uz + Uy 
1—1,2,..,8—1; 7—0,1,2,..,5—1,9 —0,1,2, ...,8 — 1. 


A square so defined is evidently Latin. For if in the z th row, the same 
number occurs twice, we must have u;u,+u, = u; U_t+Uyr, whence u= re 
Similarly if the same number occurs twice in the y th column, we must 
have Us Ug tUy=Uy Mgr Huy whence Wi Uz=Uj Uy, and since U; #0, u,=Upr 
Again any two squares say the i th andi th are orthogonal. For if in the 


æ th row and y th column of 7 th square we have j and in the same position 
of the 7’ th square we have j', then 


Uj Ug + Uy = U; 
Uy us + Uy = Uy 
(U; — Uy) u, = Uy — Uy 


Since u; — uy + 0, there is a unique solution u, and hence one and only 
one u, Thus the combination (j, 7’) occurs once and only once. 
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The work of construction of the Latin Squares can be greatly simplified. 
The integer j; (x, y) to be put in the æ th row and y th column of the i th 
Latin Square is given by Uj (2, y= Wi Ust Uy Hence if z £0, 1 j; (a, y) = 
Ji+ı (2—1, y), and j; (1, y) = ji+ı (S —1, y). Thus to get the complete set, 
we have to merely get the fist Latin Square which is called the key Latin 
Square. Again it is easy to verify that for 1 < k < s — 2,1«k'«s—2. 

(i) jik +1, k +1) =0 when j,(k, k') = 0 
(ii) Jı(k+1, k’+-1)=144; (k, k') when j; (k, k') = 1, 2, ..., s—2 
(iii) ji (E +1, k' +1) = 1 when j, (k, k')=s— 1 


The 0-th row of the key Latin Square is 0, 1, 2, ..., s—1 and we have 
merely to construct the Ist row. Knowing this row, we can write down the 
top right half of the square from the properties enumerated above. The 
other half can be written down from the property of symmetry. 


Thus when s=p”, (s—1) orthogonal Latin Squares of side s always exists. 


When s is not a power of a prime, let p,* ita Di! be the canonical re- 


presentation of s as product of distinct primes. Let r= min (pi 1, 


p^ =l, ap ni! — 1). Then we can show that a set of r orthogonal Latin 
Squares of sides exists. The proof is not difficult. Let a, be a primitive 


element of GF (nf) and let uj —0, [ua — 1, ug — ap ...., "pfi = 1) 


py ; ; 
=a; * — 2 be the elements of GF (2,3). Consider the set of s elements u 


= (uni, wai, s un). If u= (uri, usir, ... uri ) and uP = urin, ..., uzin) 
be two elements of the set, we define u£-Eu9 (urin Tous s uis + uti” ) 
and ul. y? = uri, uit, aa, Ul. un^). The set constitutes a commu- 


tative ring and an element none of whose constituent coordinates is zero, 
possesses a unique inverse. Let us arrange the elements of the set so that 
the first (r + 1) elements of the set are Uj=(Uyj, Usj -- Uy), 0 <J <7, the 
rest being numbered arbitrarily. In the z th row and y th column of the 
Jj th square let us put the subscript of the element given by Uy Us + Uy, 
(«= 0,1, 2, ...,8 —1; y=0,1-2...,s —1; j=1,2,...,r). Note that for 
the considered values of j, u; possesses a unique inverse as also Uj— t; if j255' 
and both lie between 1 and r. It is easy to establish now that the Squares 
are Latin, and any two of the r Latin Squares so constructed are orthogonal. 


4. Construction of BIBD. The points of PG (N, p”) may be identified 
with varieties and the points lying on a t-flat, 1<i <N may be said to con- 

+H. N+1_ N—t+1 
SNH MESURES (s 1)...(s —1) 


: block. Thenv= WOLDIERO an 
stitute a bloc e mei (EI) (s—1) ' 
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siti_y 
k = number of treatments per block = =r 
tions = number of t-flats passing through a fixed point 
(Py — Po) Py —P)..(Py—P,.,) _ Hi | 
(Pa Py) (P,— P1) (P, =P) sh 
à= the number of times a pair of treatments occur together in a block 
= the number of ¢ flats passing through a line 
(Py — P4) (Py—P,) ... (Py — P, .) 5 ($+ — 1) (s! — 1) 6 (N, 5,1) 
(P, — P) (P, —P4).... (P, — P, 1) (s *1 — 1) (sN — 1) dc. 
BIBD's can also be constructed from EG (N, p^). We identify as before 
the points of the geometry with varieties, and points in a t-flat where 1 «UN 
with varieties in a block. The parameters of the design are v = sY, b = 
DIN, s, t) — à (N — 1, s, t), k= number of points in a t-flat—s', r = number 
of t-flats passing through a point 
(s¥ — 1) (sN-1 — 1)... (NE — 1) - 
(s—1) (571— 1)... (s —1) and à = number of ¢-flats pss 


,T = number of replica- 


o 


$ (N, s, t) and 


N=1 __ N—t+ 
S Wee 1 
through a fixed pair of points = ( = Jan 2) 
E. ($-1—1)...(s—1) 
From a given BIBD, it is alsa 
inthe following way. In the j t 
treatments which do not occur 


c 


b'—b,r—b—rk-—v—k 
d all treatments contained 


Let us, Ui, ..., Us; be the elements of GF (s) where s = 


n x 
can find k elements us, Uio, ..., us such that k (k p". Then if we 


— 1) differences Uy, — Us s lm 
m 
are all the non-null elements of GF(s) each occurring exactly à = k (k — 1) 


go 
times, then we can construct a symmetrical BIBD wi 
D wit pes 
Yt h parameters v=b=s, 
8 


= 1 * The proof is trivial and the blocks can be obtained 
by adding to the initial block (u 


r=k,\ = 
ip Mp e Wi) the non-null elements of GF(s) 


the clements 1 ^ 
172,8, 4 B, 6, 7, 8, 9 and 10 each occmrring twice, CHR om 
9.299.789 5 U Jg Vig m2 
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6, 10), (3, 5, 6, 7, 11), (4, 6; 7, 8, 1), (5, 7, 8, 9, 2), (6, 8, 9, 10, 3), (7, 9, 10, 
11, 4), (8, 10, 11, 1, 5), (9, 11, 1, 2, 6), (10, 1, 2, 8, 7), (11, 2, 8, 4, 8) constitute 
thé symmetrical BIBD: v = b = 11, r =k = 5, A= 2. This is à particular 
case of a more general result: If 44 +3 = p", then a symmetrical BIBD 
! exists with parameters v = b = 44 + 8, r = k = 21 + land A. Letabea 
primitive element of GF (v). Then the elements of GF (v) may be expressed 
as 0, 1, a, «°, ..., a? *. Consider the set (1, a3, «*, ..., &33). The differences 
are + (a — a), + (at? — a3, ..., + (aiti GAN) o a ME AMER e 
— land writing a? — a? = qfi, these differences 


membering that «Atı = 
can be rewritten as 
q. dq. 4; +47, 
(CELESTE E a teres 
Mem jenem i oir 
These are the 4A + 2 non-null elements of GF (4) + 8). Hence every 

non-null element occurs exactly A times. Thus to construct the design 
v = b = 19, r = k = 9, à = 4 we have merely to find the initial block. Now 
a primitive element of GF(19) is 2. Hence the complete design can be con- 
structed by adding 1, 2, 3, ..., 18 successively to the initial set (1, 4, 16, 7, 9, 
17,11, 6, 5) and taking Residues (mod. 19). By deleting this block and all 
‘tréatments contained in it, we can construct the BIBI): v = 10, b = 18, 


T —9, k = 5 and A = 4. ? 
5. Construction of PBIBD. EG(N, s) and PG (N, s) where s = p” may be 
employed for the construction of PBIBD. From EG(N, s) we may omit the 


origin and all m-flats passing through it. If we take the retained m-flats as 
our blocks and the retained points as our varieties, we get a PBIBD with the 


| , following parameters. 


v= sN—1 
= ¢ (N, s,m) — $ (N — 1, s, m) — $ (N — 1, s, m — 1)* 
d r= ¢(N—1,s,m—1)—¢(N —2, s, m — 2) 
M Ecol 
1 mn—sN—s 4 =¢(N—2, s,m—2) —$ (N —8, s, m —8) 
| Ay = $ —2 Ag = 0 
aCe a 


PG(N, s) can also be employed for the construction of PBIBD. We can 


omit the point (1 0 ... 0) and all m-flats passing through this point. Regarding 


*For the sake of convenience we define ġ (N, s, m)—1 if m=—1 and ¢ (N, s, m)=0 


ifm < — 2. 
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the retained m-flats as blocks and retained points as varieties, we get a PBIBD 
with the following parameters. 


s(s¥—1) 
(er EINE 
b —6$(N,s m) —& (N —1, s, m —1) 
k —(sn*!-—1)[(s—1) 
T= $(N—1, s,m—1)—$6(N—2, 5,m—2) 
ù = >(N--2, $ m —2) —6(N —8,5m—8) 
à —0, m = 8? (sN71—1)](s—1) n —s—1 


1.) —({%—%—1 n, na (atts) 10 
pc = a) 


Many PBIBD can be constructed by considering geometrical configura- 


tions. Thus we can get a PBIBD by taking the corners of a rectangular 
parallelopiped as varieties and faces as blocks : 


D 


Thus the blocks are (1, 2, 8, 4), (1, 4, 5, 6), (1, 6, 2, T), 
and (8, 4, 5,8). The parameters of the design are v — 8, 
n; = 1,2, = 1, n; — 8, As = 2, ng = 8 and 


(5, 6, 7, 8), (2, 8, 8, 7) 
b=6,r=8,2,=0, 


( £9 fm ‘ODOM: 0 1 0 
(py)= [0 0 8 (pi) = |o 2 0 (py) = |1 0 2 

0 8 0 i o 042 
If A = (ag) is an m x n matrix and Bapxq matrix, then the Kronecker 


product of A and B is defined as 


o 


a B a B du a,, B 
AxB= | a,B a» B as, B 
AmB am:B Gan B 


which is an mp x ng matrix. If N, and N, be the incidence matrices of two 
BIB designs with parameters 0, by, Tas ky, A, and Uas be, To, ke, X, then the 
design with incidence matrix N — N, x N, is a PBIBD. In fact when 
Ty kg — T; 2, £0, the new design is a PIBD with three associate classes 


[ 
» 
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having parameters v' = 9, os, b' = b, by, r' = r, 72, k' = ky ka n'y = 0, — 1, 
: 2 
Ty = 0, — 1, n3 = (v; — 1) (va — 1), 2^, = Ti àa, A'o = To Ap 23 = A43, and 


2,—2 0 0 
(phy) =/0 0 (v, — 1) > 
0 (va —1) (vı —1) (a — 2) 
0 0 va — 1 
(py) = 0 €w—2 0 5 
vg—1 0 (vı — 2) (v, — 1) 
0 1 7, —2 
(py) = 1t 0 71—2 


7,—2 0,—2 (v—2) (v,—2) 


Many interesting PBIB designs have been constructed in recent years by 
taking the Kronecker product of incidence matrices of known designs. 
e 


6. Construction of Confounded Designs. - Let us consider the symmetrical 
Yastorial design in which there are m factors F,, F,,  .,."; each at s levels 
where s — p", a positive integral power of a-prime. A treatment can be 
represented by F, 9 Fi" cdd 9" where Xj, Wo, ... 7, are respectively 
the levels of the different factors and may be identified with the point 
(Xis Los «+e, Vm) of EG (m, s). The effect of this treatment may be represented 


by «(X Las -es Vm) OF simply by (vj 2; ... m) also. Let the sum of the 


effects of the treatments satisfying 45; T1 +--- Haim %m=%, 1—0, 1, 2, ...,8—1 
where do dg, ... % 1 ate the elements of GF(s) and ay, Mo, ..., Qim are m 
fixed elements of the same field, will be denoted by =; The treatment 
contrast [Eroh En b haa E; saa Where lo + hi... +h ,1=0 
will be said to be a contrast carried by the pencil P(a@, aya; ..., Cim). 


Evidently, there are (s—1) such independent contrasts. Ifonly ar, ai, een Org, 


are different from zero while all other a’s are zero, then in the contrasts be- 


longing to the pencil there is symmetry with respect to the factors other than 
tis ias D te Let Plas; Aga) ees dom) be another pencil such that the rank of 


re Ayo on | 
1 (55 ane om 


+ dig 2a F + -F lim Tm — % and aa @, + ... + aom toy, 


is 2. Then a3, 2. 
nie we consider the contrasts belonging to 


have 57-2 points common. Hence if 
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these independent pencils lio Zio + ha Zt F heta Sy aw and. lao 25g 


za Xa d... Toi Xa will have the sum of the product of their 
corresponding coefficients equal to 


s? at (To Fh bo "3351 ) Oa gp se 1,1) — 0. 

Thus the contrasts of any pencil are orthogonal to the contrasts of any 
other independent pencil. It is, therefore, apparent that the (s—1) degrees 
of freedom carried by the pencil P (315 455, ..., Qim) in which only the elements 
iy +++» Gy are different from zero, belong to the (k—1) th order interaction 


of the factors Fi; Fiy no Ei è 


which these elements are different from zero and the rest are zero. These 
carry (s — 1)* degrees of freedom belonging to the (k—1)th order interaction 
of Fi, Fiy e$ Fi. Let ja;, s be elements of GF (s) such that the rank of 


Evidently there are (s —1)-1 pencils in 


dun aiz Lim 
Gy dos Aom 
aja [n oed lkm 


D 


isk. Let us assign the treatments satisfying a; Taj wak... "T jg, By = Ch, 


P= ly 2, 57 E fora" Dick and let us designate this block as well 
of the treatment effects beloiging to it by Ex.» Xp 
d 2 


is 
as the sim 
ty Ge Evidently this 
block contains smk treatments and there are s* 

Dy L z DEM =0 


- Oy, [1 ras 
H poe». a. poe. t iy t 1 í 
hn Cin Ak. ty d tk X ja 2 3 ^ 


such blocks. The contrasts 


are confounded with the blocks. 


They are carrie by the AL pencils 
$— 
Bhar, T Ag 3 + ob Me djs very Ay 


where X4, Ag, ..., ^, are element: 
degrees of freedom. 


Uy, + g Gam +... ry Gyn) 


s of GF(s), not all zero, each having (s — 1) 


For a representative contrast carried by the pencil 


POY a4 +r, 5 F e ty Biv «y À 


can be written down as 


b Both X^ 4... Boa MS oy lot + uiu l,4=0 (5.6.1) 
where 2/; is the sum of the effects of all treatments satis 
Hee Ey Ops) OE os Oia a F Ak Grm) d, 
A (31 83 + 44 Ga co ae, Tm) + 


1 dim +... d- A di) 


fied by (A, 433 + Ag Gay 
= «ji i.e. satisfied by 
ES (ag me su H Oin Cpa) Vege s 
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setting a, @ + aya % +... + Tim Tm = jy ee Uy Vy pve F dp, py = ys 


we get ry ai, F As di, de s H Ae %, = Wo, Hence D; — X (ai, IM jb 
the summation X being over all Oi iy se. A, Such that Ay i pee EA 0, ma a 


Henee (5.6.1) is an interset contrast and is confounded with the blocks, 
The actual effects confounded can only be ascertained after knowing the 
number and positions of the non-null elements in P (A, aj + As doy +... + 
Ag Gg, ss Oa, F -e F Ak Apm). 

7. Hypercubes of Strength d. Let (ay, a5, ..., p) be an ordered set of m 
elements such that each a; can assume s values. There are, therefore, s" 
such sets. A subset of sf of these sets will be said to constitute a hypercube 
of strength d and will be denoted by (m, s, t, d) if all combinations of any d of 
the coordinates occur equally often, s'77 times. When s= p", a hyper- 
cube of strength d can be constructed when we can find sets 
(dis aj, «++, Qu), i= 1, 2, ..., m where each a, is an elements of GF(s) such 
that any d of mao are independent. Let yj, We ++ p, each belong to GF(s) 


and let a; = z ayj yj Then (vj Lo, ..., %,) constitutes an (m, s, t, d). For 
take the 7, th, i th, ..., ig th coordinates. Then the equations 
1 Gia Ya F Gia Da Toy d gg 
Eno? an 0 


M = 


will be satisfied by st” distinct sets of values of (yi ys ..., y) and thus 
this particular combination of the i-th, ith, ... i-th coordinates will 
occur st? times. 


An (m, s, t, d) can be utilised to construct an (s", sê) design — an s" 
design in which a complete replication consists of s"—' blocks with sf treat- 
ments to each block— such that no main effect, first, second, ..., (d—1)-th 
order interaction is confounded. For the array will contain the sets 


(4,5 C AL) am)» A= (eee aa 


and the contrasts carried by the pencil P(b,, ba» -++ bm) Will be confounded where 


93 a21 Umi dy 
1» agg “ma ese 
a i 
bm 
ay an Amt 


Supposing that among bs biy ..., D there is one null element, while the 


rest are zero, it will mean 
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etg ai Ga bi, 
saa BEN bi, =0 

NG aaa AA ; 

Gt aint Uo Giat bia 


ie. the i-th, i-th, ++» dgth sets are dependent, contrary to hypothesis. 
Example: Construct a (34, 3?) design such that the main effects and first 
order interactions are not confounded. 


In GF(8), the four sets (1, 0), (0, 1), (1, 1), and (2, 1) are such that any two 
of them are independent.’ Hence the contrasts carried by pencils P(b,, ba, bs, b4) 
which are confounded are orthogonal to (1, 0, 1, 2) and (0, 1, 1, 1) and are 
therefore P(2, 2, 1, 0), P(1, 2, 0, 1) and their generalised interactions P(2, 0, 1, 2) 
and P(2,1,1,0). The initial block is (A, u, A+ u, 22 H u), à u = 0, 1, 2 and 
the other blocks can easily be constructed. 


sk 


8. Balancing. In (sm, s*) designs, the ic Sets of (s — 1) degrees of 
freedom which are confounded in any one replication can be chosen at will. 
When, however, more than one replication is available, a better plan is not to 
confound the same degrees of freedom in each replication but to Spread the 
loss of information as uniformly as possible among the interactions to he” 
confounded. If each of the (s — 1)*-1 sets of (s — 1) degrees of freedom 
carried by the (k — 1) th order interaction of the factors E qs F, is 
confounded in 7, replications and remains unconfounded in 5 replications 
then we say that the interaction Fi Fi Lon F,, has been bal 


anced (there 
being uniform loss Ti | (Tara) on every degree of freedom belonging to this 
interaction). 


Now there are Te, (k—1)-th order interaction, carrying in all mo 
sets of (s—1) degrees of freedom, If each of these sets js confou 
replications and remains inconfounded in Ta replications, we say 


y (s—1)*-1 


We can achieve complete balance over the highest 
a (s™, s) design by taking r — (s — 1)"-1 replications 
distinct pencils TE Cry, ata dg), à; 5 0, i = ihy 
key block the treatments satisfying— 


0121 + E + am c, = 


Thus in (8%, 8) design, an arrangement achieving complete balance over second 
order interactions is given in the next page. 


order interaction in 
corresponding to the 
+, m and assign to the 


OTI 010 ore GOL 200 60% SII 210 SIG 616 SIL 610 


666 SZI $60 SIG SII 610 LOL L00 roc Tor T00 


Ila TETEL Ito 165 IGI TZ0 150 L3G IZI OIL 010 orc 


060 026 061 III IIo Iro 010 0165 OII GEL 660 


100 105 LOL 100 L0G LOL $00 5065 GOT 050 066 051 


ISI Ic0 166 010 015 OII $56 ZBL $50 $00 $06 901 


610 GIZ SII 061 060 066 IIZ EET Ito Ito IIG CT 
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va 
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To balance all interactions for designs of the class (s7, 57-1), we take 
(£y, Bos ~., Xm) where each z; belongs to GF(s) and is non-null. We construct 
the design with the help of the (m — 1) independent pencils P(a;, ..., i) 
POA no Mm — l where ay ai + ap a +... + lim Tm = 0. Evidently 
(2, Los ..., m) and p (21, Tos ..., 2,) where p is a non-null element of GF (s) 
will give rise to the same set of pencils. Hence there will be (s — 1)n-1 
replications and the constituents of the key-block will be Os (2, Uy, sss ay), 
Where ao, &, Qs, ..., %,_, are the elements of GF(s). No main effect is con- 
founded. How many times will the (s — 1) degrees of freedom carried. by the 
pencil P (0..., aj, 0, edi, s Orany tip 0, ...) be confounded ? Eviently this 


number 7, is equal to (s — 1)?-* times the number of solutions of 
[S T Qj Ri +... + di, on = 0 


by distinct sets Tip Big ory v; none of which is zero. This latter number is 


Pal JEN ei HEN ees M1 eah M NOS 
$—1 (8) sl a) nee ge) = 


Ter Gp) ra Gs) a] 


TE aa yg} T 


1 (ti p (e 
Seed e. CU ecc eC Cas) 


SUE (aU Sy) Uo a yd Eee 


$(s—1) s 
Hence 
n= B-D- , 
8 
Thus the loss of information on every (s —1) degrees of freedom belonging 
to every (k — 1) -th order interaction is 


(s 1) 1— (yea 
$(s—1)-1 


EXERCISES ro CHAPTER V 


1 Construct Addition and Multiplication Tables for GF (33) and GF (43). 
Let w be a primitive element of GF (32), Writing the elements of the field 
3&5 Gah hi i= w i i i 
zn 0; ag where URN SEQ E 2, 8, ..., 8, (i) obtain the value of 
Ag+ (da aaa s] (ii) solve the equations 43 9-035 and s Co 9 —a. 
2 sU HA7 Y= org. 
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2. A non-null element of a field is called a quadratic residue if it is the 
square of some element of the field, otherwise it is called a quadratic 
non-residue. When s = p”, p> 2 show that the number of residues is 
equal to the number of non-residues in GF (s). Defining x(0), =0, x(a)—--1 


or —1 according as a is a quadratic residue or not, show that X y (v? —1) 
v 


= — 1, where v runs over all elements of GF (s). 


8. Let2u-+ l= p" where pisa prime, and æ be a primitive element of 
GF(p"). Show that among the u(u — 1) differences of the elements a9, 2?, 
ae, 42 "-? (i) every non-null clement occurs exactly (u — 1)/2 times if u be 
odd (ii) every non-square element occurs u/2 times aid every non-null square 
element occurs (u — 2)/2 times when wis even. 


4. Construct two orthogonal Latin Squares of side 12. 


5. Number the rows and columns as 0, 1, 2, ..., 5—1. Obtain row 1 of 
the Key Latin Square from the Addition Table of GF (s). To fill up the 
Take any element in row 1 and proceed 


remaining rows proceed as follows : 
diagonal. Ifthe initial number 


by 1 single step in the direction of the leading 
jn row"! is 0, fill each successive cell by zero; if it is other than 0, put in each 


successive cell one greater than the number in the preceding cell, remembering 
that when (s — 1) is reached, the next cell is to be filled up by 1. The remain- 
ing portion of the key Latin Square can be filled up easily as the square is 
symmetrical about the leading diagonal. The other squares are obtained 


by cyclic interchange of columns leading 1, 2, ...) (8 — 1). 
— 4, 8, 9, 16, 25 and 27. 


6. How may c married couples play a (c — 1)-round bridge tournament, 


if each man plays one round with each lady except his wife and everybody 
se. Show that a necessary 


plays with everybody else except his or her spou 
condition for a solution to exist is that there are three c X c Latin Squares, 
orthogonal in pairs. Show further that when c = 2", n2 1, a solution 


always exists. 


Do the actual construction for s 


ares of side s where s is a power of a prime, 


7. Using orthogonal Latin Squ: : 
block designs : 


construct the following incomplete 
(2) pus basttsrastiLk=sand 


(zt) 


8. Associate varieties wi 


A=1 


v=b=st+s4+1, r=k=s+4+1, 4=1. 
th residue classes mod v. Let the elements of the 
the sets (aji, doi +++ ay), t= 1, 2, <-> t be residue classes mod v such that 
(i) the members of each set are distinct (ii) among the tk (k — 1) intraset 
differences, each non-zero clement occurs A times. Then (a3; + 0, ag; + 0, 
e, dg + 0) where à = 1, 2, ..., t and 0 runs lover all residue classes mod v, 


constitute a BIBD with parameters b = vt, v, r = kt, k and à. 
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9. Use Ex. 8 to construct (i) the BIBD: v —9, b = 18, r— 8; — 4 
and A = 8 from the initial sets (0, 1, 2, 4), (0, 3, 4, 7) ; (ii) the BIBD with 
v varieties where v is a prime power and of the form 6¢ + 1, b = vt, r — 8t, 
= 8 and à = 1 from the initial set (25-1, tel A EID TD <. and 
z is a primitive element of GF(v). Construct in particular the BIBD : v = 18, 
b —26,r = 6, k = 8 and à = 1. 
10. Use Finite Geometry to construct the following BIBD : 
(i) v=b=15,r=k=7, 1 =8 
(ii) v= 8, b = l4, r=7,k=4, A =3 
(iii) 0 =15,b =85,r=7T,k=3,)=1. 


11. Show that a necessary condition for a symmetrical BIBD with an odd 
number of varietics to exist is that 2? = mod P is solvable when v = 4t+1 
anda? + à = 0 mod p is solvable whe 


Nv = 4t +3 where p is a prime dividing 
the square free part of (r — 3). Hence show that it is not possible to con- 
struct the following BIBD: 


(i) v=b=48,r=k=15,.=5 p 
(ü) v=b=2,r=k=8, 12 
(iii) u= b= 67,r=k=12,,=9 
(2) v=b=77, r =k = 20, 4 s. 


have two companions. 


days no girl will walk with any one of her school 
than once. 


elements are residue 
i 1) intraset differences, the non-nuil 
54 G9) each occur ^; times, i=1, 2, ...,m. (ii) among the 
() D a =t (u, w=1,2, o My Uw), each of «0, ODE | 
Zn; occur pi; times (iii) among the n; nj difference, Oy ad, (u—1, 2, ..., np 

w=1, 2, ...,n;) the numbers ey as, s cau) each occur p! times, [—1. 9, ... m. 

Show that (a,;--0, Tai 40, ..., a,;--0), i—1,2,.., 1 Ad 9 running Praten | 
elements mod o will constitute a PBIB with parameters vt, o, tk, k, X 
Ja: Am and (pj), ij, 1— 1, Drees aM. eer eo 


14. Use Ex. 18 to construct the PBIB . 
Mm = n3 = 8, 4 = 1, A9 (Dy 
blocks (0, 2, 9, 11), (0, 1, 4, 5). 


v=] 


7, b = 88, T—8, k = 4, 
= (6 4) and (p 


^) = (1 $) from the initial 
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15. p?xq varieties are arranged as a three dimensional lattice of points, 
p along x and y-axes and q along z-axis. If blocks are formed consisting of 
all, treatments represented by points lying in planes parallel to the a z and 
y x coordinate planes, show that the resulting design is a PBIBD. Obtain the 
paramaters of the design and show that its efficiency factor is 


(pg — 1) ] [Ga — 1) -2( — 2]. 


16. Use Finite Geometry to construct the following PBIBD : 
(i) v=14,b=28,r=6,k =8 
Ay = 1, my = 12, Ay = 0, m = 1, pin = 10 
(ii) v=b=12,r=k = 06,4, = 3, n, = 8,24 = 2, n = 8, pty = 4. 
(iii) v =b = 36, r= k =5, 4, —1, nj —20, 4,—0, n= 15, phi 10. 


17. Associate two varieties to each element 0, 1, 2, 2°, ..., v~1 of GF(s) 
where s = 2u + 1 = p^, p>8 and v is a primitive element, by means of 
lower suffixes 1 and 2. Show that by adding the elements of GF(s) in succession 
to the constituents of the initial blocks 


4 2u—2 2 2u—2 

s (Or 1,, 22, 2, a, La, Ua, e Te”) 
4 2u—2 2 —2 

(eg Ch OR coy CHER OT GAY oo 8 -) 


we obtain a PBIBD with parameters v = b = 2s, r = k = s; 


A = (8 1)/2, a = 5 — 1, n; = 28 — È, n, = 1, py = 28 — 4. 


18. In a 2" factorial experiment in 2* blocks of 2"-* plots each, if 
the main effects are denoted by 4, B, C, ..., the two factors interactions by 
AB, BC, ..., and so on, show that if any two interactions are chosen to con- 
found, the third is automatically determined by throwing together theletters 
of the two already chosen and suppressing any letters they have in common. 
Show further that the constituents of the key block are treatments having 
an even number of letters in common with the confounded effects and the 
constituents of the other blocks can be obtained by multiplying the consti- 
tuents of the initial block by treatments not present in it and suppressing 


any letters they have in common. 


19. Construct (i) a 25 design in blocks of 8 plots confounding ABC, ADE 
and BCDE (ii) a 2° design in blocks of 8 plots confounding ABC, CDE, ADF, 
BEF, ABDE, BCDF and ACEF (iii) a 25 design in blocks of 16 plots confound- 
ing ABCD, ABEF, CDEF. 


Establish that in s" design where s isa power of a prime and the 
he maximum number of factors which can be accommodated 


ffect or two faetors interactions are confounded is (s* — 1)/ 


20. 
blogk size is s*, t 
so that no main e 
(s — 1). 
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21. Construct a 215 design in blocks of 16 plots so that main effects or 
two factors interactions are not confounded. 


22. Show that using blocks of s” plots (s — p", where p is a prime), it is 
possible to test all combinations of So many as N, (but not more) factors in 
such a way that any interaction confounded does not involve less than k 
factors, if N, denotes the maximum number of distinct points which can be 
chosen in PG (r — 1, p") such that no (k — 1) of the chosen points is conjoint, 


Show in particular that for an experiment in which each factor is at two 
levels using blocks of 27 plots it is possible to accommodate a maximum 


number of 27-1 factors, if confounded interactions do not involve less than 
four factors, 


28. Construct (i) a confounded 33 design in 9 blocks of 3 plots in which 
four of the degrees of freedom confounded are carried by the pencils P(1, 1, 0) 
and P(1, 1, 1). (ii) a confounded 43 design in 16 blocks of 4 plots each in 
whith six of the degrees of freedom are carried by the pencils P(1, 1, 1, 1) and 
P(1, 1, 1, «) where « is a primitive clement of GF(4). 


24. For a 89 design in nine blocks of three plots each, the constitution of 
the blocks are as follows : (000, 011, 022), (101, 112, 120), (202, 210, 221), 
(001, 012, 020), (102, 110, 121), (200, 211, 222), (002, 010, 021), (100, 111, 122) 


3 
(201, 212, 220). Which 5 degrees of freedom are confounded ? 


25. Show that the total nu mber of systems of confounding for a s 
ment in blocks of s*es being a power of prime, is 


CAEI ET (Gar Sai) 
(s"-5 — 7) (er =o (s—1) 


experi- 


26. Show that in the case of 5* design in s blocks of s plots each, 
4 complete balance over the first and second order inter 
replications are sufficient, 


to achieve 
actions (s — 1)? 
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29. For the s? design in s? blocks of s plots each, show that in order to 
achieve balance over first order interactions only, the number of replications 
necessary is (s — 1) when s is the power of an odd prime and 2(s — 1) when 
s i$ a power of 2. 

30. Show that if we take (s — 1)? replications of a s4 design in s? blocks 
of s? plots each in which the key block of a particular replication contains 
the treatments 

v + xw + Be, = 0 

Wy Bo lay + da 1», = 0 
where c, d are fixed non-null elements of GF(s), s = p" and « and f are all 
possible non-null values, a complete balance over the second order inter- 
action is achieved. 


31. Construct 3 replications of a 4? design in 16 blocks of 4 plots each 
achieving complete balance over the first order interaction with a loss of 
1/8 of information. 


2 


CHAPTER VI 


SOME SELECTED TOPICS IN DESIGNS OF 
EXPERIMENTS 


l. Missing Plot Technique. Often in experimental work, one or more 
experimental units may be missing through accident. The experimental design 
with which we started thus gets modified. The usual Least Square method 
can be applied to meet the new situation. The missing plot technique, however, 
seeks to provide numbers in the missing plots so that the ordinary 


procedures 
of the original design can be used. 


Suppose that y, yo, ..., Yn are the observed yields in the existent plots and 
Jp Ws --.. Vp are the yields that might have come from the missing plots. 
Usually, the expectations of these yields are linear functions of some para- 
meters 0,, Os, ..., Om- We give below the ordinary least Square analysis and 
the analysis with the substitution Vi, Vos «.., 0, In the missing plots. 


(a) Analysis with existent observations 


IE 31033 verre y 0, 
22. 6 

E(y)=E | IJ = | 1» —A0.. (6.1.12) 
Ya Ani Ig sasas Anm On 


Sum of Squares to be minimised — 


Y— A0 (y — A8) — y'y — WA y --0'A'AQ (6.1.22) 
Differentiating (6.1.22) 


with respect to 8 and equating to zero, we get the 
normal equations 


A'y — A A 8 A (0.1.32) 


The residual sum of Squares is equal to 


yy—6'A’y (6.1.42) 
Suppose that we want to test the hypothesis DS 05,54... — 0, — 0. 
Then we have to ealeulate ; ith 6, +1=9,, A RI 
76,—0. squares and (6.1 tan) 
against (6.1.4a). To 
d as follows : — 


=A, $ (6.1,5a) 


War Ange ] 0; 
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Sum of squares to be minimised is equal to 
(y — Ay Y (y — Ai 9) (6.1.6a) 


The normal equations are 
^ 
Any- ALA > (6.1.72) 
Conditional residual sum of squares is 
^ 
yy—Q A,y (6.1.82) 
, 2), in the missing plots. 


a 


(b) Analysing with substitutions of a. %4 


[51 B b11 b12- Vim 6. 
E(y —A0,Ex)—E|: | =] - ?* | —B6 
& bra bre- br AT 
k J Ia Vj: Orm 8, J 
A (6.1.1b) 
Sum of squares to be minimised is equal to p 
(y — A9) (y — A 0) + (x — B 0)' (x — B 0) (6.1.2b) 
Normal equations are 
A’y + B'x = (A A + B' B)8 (6.1.3b) 
i y A+ x! B —0' (A'A + B' B) s z (6-1-4b) 
B= de (A' A +4- B' B) ^ (6.1.5b) 
dx 
Residual Sum of squares is equal to $? where 
Sp =y y +x’x —0'(A’ y + B’x) (6.1.6b) 
Differentiating S? with respect to x and equating the derivative to zero, we get 
$ EDEN Cien A 
2x — —— (yd B’x)— B@=0 (6.1.7b) 


Remembering (6.1.5b), x = B 6 (6.1.8b) 


Substituting (6.1.8b) in (6.1.8b), we get the same normal equations as in 
(6.1.32) “and substituting (6.1.8b) in (6.1.6b), we get the same residual sum 
a). Hence we have Fisher's rule : 

If variates vj. Wos ... 2 be missing, we can proceed as if they are present 
vided that we put for them values which would minimise the residual 
The residual sum of squares so obtained is the actual residual 


of squares as in (6.1.4 


pro 
sum of squares. 
sum of squares. 

Let us find out the conditional residual sum of squares. Here the sum of 


squares to be minimised is 
(y — As 9)! (y — Ar 6) + (x — Bi WW (x — Bi 6) (6.1.9b) 
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The normal equations in the case, are 


(A^ y + BY x) = (A', A, + B', B) $ (61.105) 

whenee 
p, 29 (A^, A; + B’, Bj) (6.1.11) 
dx 

The residual sum of. squares $7? in this case is 

yy+x'x — 6 (A^ y B’, x) (6.1.12b) 
Equating the derivative of (6.1.12b) with respect to x to Zero, we get 

iuc CU POUR 

2X — dx (ALY + BY x) — 10—0 (6.1.13b) 


Remembering (6.1.11b), we finally get 


E 
K=B 90 (6.1.14b) 
Substituting (6.1.14b) in (6.1.10b), we get the same normal 


equations as in 
(6.1.7a) and substituting (6.1.14b) in (6.1.12b) we get the same conditional 
residual sum of Squares as in (6.1.8a). Hence 


Sum of Squares due to hypothesis = Min (conditional Residual Sum 6f 


Squares) — Min (Residual Sum of 
s Squares) 


The sum of Squares due to hypothesis caleulated from the 


completed table 
will, therefore, always have a positive bias, 


is u ka; 1L Tj» itis well known 
Tj are respectively the grand mean, 
and difference between 


; i» the actual ith block total by T;, 
7, actual grand total. By (6.1.8b) 


whence» = "Bit tT; ~6 l 5 
Ec S (6.1.15) 
If the yields of the jth treatment in the itb block 


i-th block be missing and if v and w 
we have if 7 zd 253 


and j'-th treatment in the 


are substituted for the missing yields, 
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— sp On. Upon Ene aru 

t T tr 
AA Bem. Tw G 4-2 4-0 

t T tr 


c 


pu 


Whenee we get easily 
(r — 1) (t — 1) [r B; +47, —G] —[r By + tT, — 6G] 


v= (r—1pg(—19—1 
wo SHV C= DPB, 17; —G] — [rB; 17; — G] 
(FD n ME 


The estimate of the missing plot can also be obtained by utilising the method 
of covariance analysis. Let x = 0 and y = the actual yield for the existing 
units and æ = — 1 and y = 0 for the missing unit. The best estimate of the 
missing yield is simply the error b regression coefficient. Thus for the random- 
ised block with one missing unit, we have the covariance analysis D 


S (a?) S (xy) b 
1 1 B; G 
locks Ss re cens 
Blocks ; E : + z 
"Tr E 1 mon G 
r 'Treatments mie zi = S 2 
To hls G er BU KANG 
Error —1 — 1)jrt, — + 2 - =] i i 
et (DG cour te A eme 
"Total jue e 
K TA rt 


If the yield of the plot in the ith row and jth column containing the kth treat- 
ment in the r x r Latin Square be missing and if Rj, Cj, T; and G denote 
respectively the actual row, column, treatment and grand totals we have w, 


the missing yield given by 


bles Rita o Cw ty Te +æ 2(G + a) 
T 7 T 72 
whence 
f Qo TT C; + Ti) — 26 
(r-)(-—2 

The covariance technique would also give the same result : 

S(e?) S(xy), b 
Rows (i/r — 1/7?) — Rilr + Gr? ve 
Columns — (1/r — 1/7) — Cj[r + G[r* 
Treatments (1/r — 1/7?) — Tyr + Ghr? at. 
Error (r—3)(—23)P, [r(R + C+.) 26], Us tO an aa 

Dir 2 


G[r* 


Total 1—1/r° 
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When performing the analysis of variance, it is also unecessary to obtain 
the conditional error sum of squares. Thus if we obtain the actual block sum of 
Squares and total sum of squares, then by subtracting from the total sim of 
Squares, the block sum of squares and error sum of squares, we obtain the 
actual treatment sum of squares. In the randomised block experiment with 


one missing yield, we have seen that (6.1.1 5) provides an estimate of the missing 
value. Evidently, 


Bla) — + y 


2 
E zr — 1) EP(r—1)- rt —1) — 266 1) — ao —1) ^ 
d: (r—1) (1 E 

OT — E r 
—(—318-1* (6.1.16) 


The variance of the estimate of a difference between a treatment with one 


mussmg unit and a treatment without any missing unit = J” Tika =) 
T r 


r+t—1 247 —1) — 9 E 
L a pr HET) aq) 
D (r—1)(t—1) e 
t 
1 ————— 1 
y T 26) i] (6.1.17) 
Whereas the variance between two treatment means without any missing 
2g? E A 5 
unit is—— . 


T E 


Average variance of treatment differences 


202 na a (£—1) 4 2% (t=1) (t —2) 
UT 3( —1)(—1) DL ERUNT DER 
t(t — 1)/2 
20? 1 
EXE HT PN 1 
= + E (6.1.18) 
Hence the loss of efficiency is 1/[(r — 1) (t — 1) + i]. 
The actual treatment sum of Squares 
UA [t+ LB E BY uu. KAB Laj oa. +B 
t 
T5 TA e. +(Tj+@)? 4p... apt (Ge =| 
= - 
tr 
Bh ee +BY BE IB 
i E Biti) 
E [ T9, n TE... TT, (G La)? 
T zn = ] — Bias. 
Bisse Di 4 Ee — 2Ba t=1 B, 7 
1(1— 1) t CEDE caes | (oo = 


t 
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The result is quite general. In fact, the sum of squares due to hypothesis in 
the analysis of variance table with substitutions for the missing yields will 
always be greater than the actual sum of squares due to hypothesis and if this 
turnf out to be non-significant compared to the error sum of squares with 
appropriate degrees of freedom. the actual sum of squares due to hypothesis 


will also turn out to be non-significant. 


2. Fractional Replication. In factorial experiments, when the number of 
factors is large or one or more of the factors include a large number of levels, 
combinations. Thus with 6 factors each 
729 treatments. In situations like this, 
ay not be sufficient to have even one 
if properly employed, is useful 


we have a great number of treatment 
at three levels, we have as many 2s 
the resources of the experimenter m 
replication. The device of fractional replication, 
in such situations. 


Let us consider the symmetrical factorial design 
h i=l, 2 o k 


s" where s is a power of a 
prime. Let P(aj, di» -+++ lim (6.2.1) 
where each a; is an clement of GF(s) denote k independent peneils. Then 
in a (1/s*) replicate of the s” design, we shall take only the smk treatments 
satisfyiie 
jg Vy H Aga to Poe H lim Um = Gp j= 1, 2, -..,k where each ZA belongs to 
GF (e). Usually cach dj, on the right hand side is taken equal to %, the 
null element. Thus of the s” treatments 
smk treatments satisfying g 
t an B= 0 t=] 2; es k (6.2.2) 


yy Vy E lito H E lim Tm 


4 E 
(ag, Uat- Un) We retam only the 


Let P(b;, b, ... + Ym) be a peneil independent of (6.2.1). Let 


> denote the sum of treatments satisfying 


H aab Ug Un — X 
jo dig I ese ljm Um = ij J = Vs 2, oes k 


rried by the pencil P(b;, bos «+ bm) may be written 


Then any contrast ca 
in the form 


aes ly S > 
e, Doa; Pigs Iya SEWA pesca fas Vip Kig n Mi, ar 


(6.2.3) 


) and S stands for summation. In terms of the 


ee C 


s 


where f, : 
retained treatments, (6.2.3) can be written as 


Euh Diag 0... 0 (6.2.4) 
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The contrasts carried by any one of the s* pencils 


S 
A 
Tiya «sey Ob DNE nh RE iis ts ta, (6.2.5) 


and they become in terms of the retained treatments exact] 


y (6.2.4). Summing 
up (6.2.5) over all values of 24, Ao, ..., ^, we get . 


SS [exe Bo, uus, Za,_,0...0 
ap ke te parere Ea +... + Dey — (300 ...0 +... 4 Ba, 00 ... 0)] 


+ hf TE c ) 
= sf L X4 l Xo, oote +0, Xa oo | 
TÓC EN. I LJ)[ ] 


ms X .-- 15200...0 4L ... HLE 9,1 0...0] 


Thus each pencil is one of a set of s* alias set of pencils such that the sum 
of similar contrasts carried by all of them is estimable. The result of using 
a fractional replication is that the contrasts are no longer Separately estimable, 
but the sum of the contrasts belonging to any complete alias set is estimable. 


| /—1)-flats in the (m—1)-flat »,—0 in PG(m, s) and 
is equal to $ (m — 1, m — fp — 1,5) =$ (m—1, k— 1, s) which wirtten out 
in full gives 


(s*— 1) (ska ~1)...( 


In an s" design in 1/5" replicate, 
confounding for effective eliminatio: 
quite straightforward. We choose gi 


it is sometimes 


n of soil heterogeneity. The process is 
ndependent pencils 


necessary to resort to 


D (Cs, eps, +13 Cim) Y = 1, 2g (6.2.6) 
which are also independent of (6.2.1) and assign to the 58 blocks the treatments 
m 
Z -6i-ne. pk 
De 
m 
E FG mud =1, 2, d 


ee 


DESIGNS OF EXPERIMENTS 105 


The confounded effects belong to së pencils and their alias sets. The sum 
of the contrasts carried by any pencil independent of these k + g pencils 
and its aliases, is estimable. Let us illustrate this by the following example. 


TABLE 


Half Replicate of a 2° Experiment with Confounding to Reduce Block Size to 


Eight 

Block I Block II Block III Block IV 

Ly ay V3 
Ua Us Te Wa Cg Ve Q4* Tg Ue Ty Vg Te 
o 9 O ( 9 0 TOEIC IO Um mn m oto 
(a d. IM T (yy al (y 31. 4 i0 © 
( i | G Wi ib 04410 sk MI Oo @ i 
() dou og il O u d (y 707. f) io gh ® 
i (0 0 i a 0 ( xt 9 0 0 1 beg 3l 
TO al @ i m i qo imi 10 0-0 
Ee 0 EO Sh) Oo G0 (- 3, 3H ibo rit 
dox EET n cw nu ib «9.0 Ola O 


te of a 29 design involving six factors 
]s 0 and 1 respectively. The retained 
= 0. Confounding was further 
ded pencils being P(0, 0, 1, 0, 


This is an example of a half replica 
A, B, C, D, E and F each at two leve 
treatments are 2, + @ + Us + 94 + a@5 + Te 


resorted to for reducing block size, the confoun 
0, 1) and P(0, 1, 1, 0, 1, 0). Thus CR =ABDE,BCE=A DFand BEF = 


A C D are confounded with the blocks. The alias of a main effect is a five 
factor interaction and that of a first order interaction is a four-factor inter- 


action. 


8. Confounded Asymmetrical Facto: 


factorial design, the m factors Fi, Fo, .. 
number of levels. In such a case the experiment is known as asymmetrical. 


Problems of reducing block size by confounding only high order interactions 
do not always lead to solutions and even in favourable situations, the problem 
of construction: of such designs is more difficult than the symmetrical case. 
Thus in the case of s, X Sa design, except in special cases, it is not possible 
to arrange the sı X Sz treatments in sz blocks of 8; treatments each, without 
confounding some contrast belonging to the main effect of the second factor. 
Hence we try to have a number of replications such that the loss of informa- 
tion on every degree of free effect of the second factor 


dom belonging to mam. 
is uniform and the loss of information on every degree of freedom belonging 
to the two-factor interaction is al 


]so uniform. This leads to a balanced arrange- 
ment. A few special cases are described below : 


x s, The arrangement in Table 6-1 


rial Designs. In the general case of a 
., Fm usually do not all have the same 


(i) The case 2 
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« 


sts. 
2 
(z—'s) 


Li Sc. 
(g—s)” 


tst 
(1)” 
25 


Zs - gt 5 r— 9s) rd $5 (T — 8s) eee 355 c nd Bs t+ e zo 339 A I 3po[g 
2-9 TIY, 
ts tsy se "sap Loe toe [acr] E Bon Bln I'$5 ee Igo TI% g 
By AT [^ Bp E z S z á I emn I I z 
35 es. z I ene ze z I Loa wee I I 
- o BA Jo [9^9'T Tg Jo PAT 
Beery Aes (T—*s) pte (r— 9s) v — "s SS T Pt este 5s d Se ed I POL 
2-9 TIV, 
———À——————————M——M——MÁÉ—Ó ——— ——d—— nae 
[^ ee Z one it ose E^ op A I Bs "m z I [4 
L^ ae. Bg m L^ aoe z aod [4 PA I zm I m I 
LLL DRM LI AA D M C M I M pu i 
3A Jo PAT Tq jo PAYT 
m ————Á—Á————————————Á———À ÁÓHÓOÓ€——ÀÓ 
fgg -E*s(r— s)'r-- 5s(1— 95) c7 Boyt d fs p H- US sg Eo [4 t POT 
[-9 @Iavy, 


$$... 86 
(2—i5)? (as) 
BSy... By 
(1) (1) 
zc z 


[E [4 


ZJ JO paT 


TASE 15. I 
£ (9)? (gs)? t $ 
rs. ag r 
S (D^ (9^ | 8 
Z I qe iP ik [4 
I Es Yn oT I 


ly jo [oae 
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will lead to a confounded arrangement having parameters v = 2 X Sa b = GA 
k = 2, r=s, Here the main effect of F; is unconfounded in each of the sz 
replications and there is a uniform loss of information 1/s, on every degree of 


freedom belonging to either F, or F, Fe 


(i) Thecase 8X 555278. A LatinSquare of side s, can be constructed with 
(1, 2, 8, ..., 8) in the leading diagonal. Let 


O11 “12 m Als, 
«21 «22 . €25s 
sul [v1 [E 


be such a Latin Square. Then the arrangement in Table 6:2 


will lead to a confounded arrangement having parameters V = 85; b = Sh 
k =8, r = ss. Here the "main effects of F, are preserved in each of the s 
replications but there is a uniform loss of 1/s, on every degree of freedom 


belonging to F, and Fi Fa 
o 

(iii) 0 = 5, X5 where s; is a power of a prime and sı < sa +1. 

f orthogonal Latin Squares whose first rows 


are l, 2, «++» Sa for constructing confounded arrangements with v = $1 X Sa 
b = s, k = Sp 7 = Sy preserving the main effect of F, but with uniform 
loss of information 1 [sı on every degree of freedom belonging to F; and F4 Fa 


In this case we can make use 0 


2nd Square (s, — 2)th Square 


1st Square 
1 Darana S2 1 PAN Ss 1 2 83 
2 2 (2) (—2)— 59 (s;—2) 
Q) (tes «(D 2 (2) teg 1 “os, 
2) (Urt l aa bakang ah (51—2) 
QE cen "EM a alte 


The balanced arrangement would be as given in Table 6-8 


bor s. design, if 3,41 = DLL $778 (say) we can 
In thes, X X X fo ee -1, The s™~! combinations of the last 


m 
write it as 8; X 9 Xe MH XS 
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(m — t) factors can be divided into s! groups of 5"-- 
these s' groups as the s' levels of a fictitious factor Fu. By arranging the 
division of the s™—' sets into s! groups of 5"-!-! each such that there is 
uniform sacrifice of information on interactions of the factors F 
F, we can obtain balanced confounded arrangements in the 


X 8 X s"-'design. For example in the 2 x 32 design, 
division of the 3? design involving F, and F}. ` 


each. We can imagin 


t» Fipo a 
SEEE 
each of the following 


5 


I 


II 
(1,2), (2,1), (8,3) (1,1), (2,2, (8,3) 4| 
(2,2) (8,1, (1, 3) ; (8,2, (2,1) (1, 8) : | 
(8,2, (1,1), (23) (2), (8,1), (2,8) | 


into 8 sets of 8 each will confound 2 d.f. belonging to F,F;. Associating thes | 
sets with the levels of a fictitious factor Fa, we can from case (i) of section 2, 
obtain the design given the Table 6-4 Riad 


a 


TABLE 6-4 

Blocks 128 456 7.8/9. 101112: 18-14 15/9 16 iy as 
Levels 
of Level of F, 

F,F, 

Sh ZERO a LS oem CIC RM NARI Lt 
21 TTE, 2:29 8838 222 TET, 333 
TET 222 888 eat al 111 838 222 
81 838 Fatal 222 833 222 PTT 
22 128 1238 128 218 218 218 
12 281 231 231 132 182 182 
82 812 812 312 


<=> 
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Omitting the blocks, 1,5, 9, 10, 14 and 18 we get another good design with 
F, F, and F, F, Fs partially confounded. Many papers have appeared 
in recent years on construction of balanced confounded designs in the asym- 
metrical case but a fuller discussion is not possible in a book of this size. 


4. Weighing Designs. When we want to find out the weights of a number 
of light objects, it is possible to increase the precision of weighings by taking a 
suitable combination of the objects instead of weighing them separately which 
is the usual practice. Let o? be the variance of a single weighing. Then weighed 
separately, each of two objects will be estimated with variance af, If however, 
both objects are placed on one pan and the estimated weight is y, and the 
two objects are placed on different pans and y, is the estimated weight placed 
on the pan containing the lighterof thetwo objects to obtain balance, then 
the estimates of the weights of the two objetes are (y, + y;)|2 and (y, — y;)/2 
and the variance of each estimate is c*?/2. Evidently the second procedure is 
superior to the ordinary method of weighing the two objects separately. 

Let Tiy take the value 1, if, the ith object be included in the ath weighing 
in the right pan, the value — 1 if the ith object be included in the ath weighing 
in the:left pan and the value 0 if it is not included in the ath weighing. Let 
there be p objects having true weights £;, Ba... Bp and let there be N weighings, 
the observed results being yz, ys ..., Yy. Now 


o 


Yı Tu Tg. Bi 
E()-E|^|-|?n?w"» | | P Lx... (6.4.1) 
UN. TIN Cy. Bp 


In order to find estimates of B the usual method is to obtain those B which 
will minimise | 

(y — XB) (y — XB) = yy — 2 8’ X'y + g' x’ xp (6.4.2) 
Equating the derivative of (6.4.2) with respect to B to zero, we have estimates 
[3 of B given by ; 

X’y =X’ X 8 (6.4.8) 
Assuming x'x to be non- Singular, we can write B= (KA) IK y. Evidently 
E (Ê)= B, E (8—8) (B — BY =2(X’X)-1. Let Cadet (X'X) = ley 1 
i,j = 1, 2, ..., p and let C;; denote the co-factor of cy in C. Further, we denote 
(x'x)^! by- (e). Now : 


= ĉn Cy — > C; e; Gij 
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The second constituent on the right hand side is a positive definite quadratic 
form ; hence the minimum value of ci! will occur when Cry = 4g =.. 


5 = Cip 
= 0 and the absolute minimum value 1/N will occur when further aj = 
+1, j= 1, 2,..., N. Thus when KK — NI 


p» all the objects will be weighed 
with maximum precision. 


This introduces the desirability of constructing Nxp Hadamard matrices X 
with elements + 1 or — 1 such that XX=N Ip. Since 


OFANA ATTA: M ICA! (A A — RAA O 
6. — a) A—AJ — [A'—A' 14—A] 10 2AA 
and l 


1-1 is a Hadamard matrix, a Hadamard matrix of order 2* 


always be constructed.. 


can 


1 A 
X X)y-1— P siis ———— 

(E X) SN 3 T—) L (r — 3)rk Es o 

Hence the variance of the estimated weight of each of the v objects is 
o o (rk — A) | (r—2) rk 

The efficiency of a weighing design is often defined as Min det (X’ X)-1 

det (X' X)-1 

det (X' X) 


= Max det (X Xj When a Hadamard matrix H N41, Of order (N --1) exists, 


gn of maximum efficiency 
constructed. We.have merely to “sweep out” 
omit the first row and ist column and replace the 
resulting matrix of order N'by+1. The deter- 
minant of such a matrix is CN 1) 2 V in absolute value. Accord- 


ing to this definition of efficiency 


a spring balance weighing desi 
ings of N objects can be easily 
the Ist column of Hy, 
non-zero elements of the 


involving N weigh- 


1 0 1 
M 1 1 
1 1 0 


is the best weighing arrangement with three weighin 
a spring balance. The arrangement can also be identi 
treatments in a half replicate of a 23 design: a, 
combination 0 0 0. 


Es of three objects in 
fied with the retained 
+ ® + e= 0, barring the 
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-5. In section 6 of Chapter 3, we have shown that a symmetrical BIBD can 
always be converted into a Youden Square. Another method which is more 
helpful in actual construction is to write down the symmetrical BIBD as 


follows : 


Positions LE Loe aa a AN b 


As a measure of “ non-Youdenness " we shall take M = = M(i, j), where 
hJ 
M (i, j) is 0 if the jth variety occurs in the ith position and 1 otherwise, i —1, 2, 
+.) ks J = 1, 2, ..., 0. Fora Youden Square M = 0. If M > 0, a particular 
variet$, say j, will not occur in the ith position and will occur twice in some 
other position i’. The situation is represented below 


e 


à Block 
Position a y Q 
T E 
i J Jj 


where both treatments j, and jẹ are different fromj. Interchange of j, and j in 
block y will normally lead to a reduction in the value of M except whenj, 
already occurred in position 2’ and j; did not occur in position 7 except the one 
place indicated. In this case, after the interchange, ja occurs in two places in 
position i' and is absent in position i. We may repeat the procedure with j, and 
it can be shown that ultimately M will be reduced by at least 1. Thus by 
repeated applications of this system of interchanges, we can reduce M to zero. 


By a slight modification of the above arguments, one can prove easily 
that a balanced incomplete block design, in which the number of blocks is 
v m, an integral multiple of v, the number of treatments, can be converted 
into a design in which every treatment will occur m times exactly in each of 
the k positions. In the notations of section 2 of chapter 2, L =m Ey, NN’ = 
(r —d) I, + AE, Hence F of (2.2.4) is = I,-— = E,,, the eigen values 


of which are 0 and x having multiplicities 1 and v — 1 respectively. Hence, 
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c? 


all elementary treatment contrast will be estimated with variance 


Further (2.2.11) is satisfied and thus the estimates of row contrasts are 
orthogonal to estimate: of treatment contrasts. 


EXERCISES To CHAPTER yi 


1. In the following 5 x 5 Latin S 


quare Experiment, certain plot yields 
marked with an asterisk are missing 


AE CEDA E 
B CD EA 
CD Pd AUS 

D E A' Bc 
EAEI GD 


Obtain the variances of different treatm 


3 s ent comparisons and calculate 
the loss of efficiency due to the non-availability of yields of the three plots. 


2. Give the analysis appropriate to an s x s Latin Square when a single 
row, column or treatment is missing. 

8. In a balanced incomplete block design, caleulate estimates of different 
treatment comparisons with their variances when (i) the yield of one plot is 
missing and (it) the yields of all plots in one particular block are missing. 
Calculate the loss of efficiency in each case. 


4. In a split-plot experiment arranged in r randomised blocks, having 


levels: of subtreatments B, the yield 
d error of the estimate of (i) difference 


plot yields are not Separately available, 
known. Obtain the best unbiased linear estimates of the Separate yields. 
Compute the variances of different treatm, 


ent comparisons and find the loss 
of efficiency due to the mixing up of the plots, b S 
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6. In the following 4x4 Latin Squares experiment, the yields of the two 
plots marked with asterisk got mixed up 


E A B c D 
c * D B 
* c B oe ak 
B D A c 


Obtain best unbiased linear estimates of the separate yields. Obtain 
formulae for the estimates when the mixed up plots belong to the same row, 


column-or treatment. 
7. Below is given an one-quarter replicate of a 28 design where confounding 
has also been resorted to for reducing block size : 


Block I: ab, acf, beg, aeh, cefg, cd, abcdgh, bdf, abdefg, abcefh, gh, 
? bcfgh, adfgh, acdeg, defh, bedch. 


Block II: bef, aeg, (1), beh, abcefg, abcd, cdgh, adf, defg, cefh, abgh, 
H acfgh, bdfgh, bedeg, abdefh, acdeh. 


Block III: ce, bef, acg, beh, abfg, abde, degh, acdef, cdfg, fh, abcegh, 


aefgh, bedefgh, bdg, abedfh, bdh. 
Block IV: de, bcdef, adg, bdh, abcdfg, abce, cegh, aef, fg, cdfh, abdegh, 
acdefgh, befgh, beg, abfh, ach. 


Identify the alias sub-group of pencils and the confounded interactions. 
Write down the aliases of the main effects and two-factor interactions, Assuming 
interactions of 8 or more factors to be negligible, give the partitioning of the 
total degrees of freedom in the analysis of variance table. 


8. Ina cotton mill, one of 5 spindles was found to be winding defective 
weft. In detecting the cause of the defect, the four component parts of the 
spindle designated by the letters A, B, C, Dare interchanged. The reconstructed 
spindles are then tested five at time over five periods P. You have the techno- 
logist's assurance that there is no interaction between the factors P, 4, B, C 
and D. Design a 1/5? replicate of a 55 design choosing a suitable alias sub- 
group for testing the effects of the various factors. 

9. Construct a 45 design in 1/4 replicate with confounding to reduce the 
size of. the block to 16 plots in which only 3 degrees of freedom from one first 
order interaction-are confounded. Give the constitution of the principal block 


and the partition of the total degrees of freedom in the analysis of variance. 
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10. Show that if the alias subgroup involves only interactions of the 
(t+-k—1)-th and higher orders (17K), the alias of an interaction of order (k—1) 
would be (? — 1)-th and higher order interactions, Hence or otherwise show 
that if an alias subgroup involves interactions of only the fourth and higher 
orders, the main effects and first order interactions are estimable, when 
interactions of second and higher orders are negligible. 

11. Show that if we adopt a hypercube (m, s, t, d) of strength d as x 
replicate of a s" design, interactions upto order k — 1 where 2k <d+ are 
estimable when interactions of order equal to or greater than d — k are 
negligible. 


12. Six factors are available each at three levels. Show that it is possible 
to arrange the experiment so that (i) it is in 9 blocks of 27 plots each and 
forms (1/3) replicate of the 8° design (ii) the aliases of main effects and two- 
factor interactions are interactions of at least four factors and (iti) two of the 
degrees of freedom between blocks may be interpreted as an interaction of not 


less than three factors, while the remaining six degrees of freedom between 
blocks are interactions of at least four factors. 


Write down the constitution of the key block and indicate how others may 
be found. If residual error is estimated from four-factor interactions, ascertain 
the number of degrees of freedom that will be available for the purpose. ^ 
. 18. Construct a confounded q X 22 design in blocks o 
replications such that the loss of informati, 
(q — 1) d.f. belonging to 4 B C is 4(q — 1 
be generalised so as to give q X 2" designs 


f q X 2 plots in q 
on on B C is (q — 2)*[g? and in 
)/q?. Indicate how the design can 
in blocks of q x 2-1, 

14. Plan a confounded g x p? design where p is either an o 
power thereof in blocks of q X p plots in q(p — 1)/2 replications such that 
the loss of information on the (p — 1) degrees of freedom belonging to B C 
is [(p — 1) q? — 2pq + 2g]/g* and that on the (4 — 1) (p — 1) degrees of 
freedom for A B C is 2p(q — 1)/g*. Indieate how this design ean be extended 
to q X p" design in blocks of g x p^-iplots, ^ 

15. Let p be a prime or power of a prime and 
than p. Write the p — 1 orthogonal Latin S i 

quares of side d 

number denote level of the first factor, If we omit QE E Pee NM 


numbers q + 1, q + 2, ..., e the fa 
re j indi 
levels of the second factor, Lane numbers indicate 


in (p — 1) replications. Show that the losses of information on A, Band AB are 
0 and p/q(p — 1), ; 
16. In Ex. 15, suppose q = p and we tak icati 
uppo e k replicati 
(p—1). Show that in this case only interaction AB eR AE ro 
Caleulete the loss of information on AB. ie d 


dd prime or a 


q be a positive integer less 
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17. Use Ex. 15 to construct a p x r X s confounded design where p has 
the same meaning as in Ex. 15 andr X $s =q p by the substitution of the 
TXS combinations for the q levels in that example. Illustrate your method by 
constructing the 4 X 2 X 2,5 X 2 X 2 designs in blocks of 4 plots or 7 X 3 x 
| 2 design in blocks of 6 plots. Show that the loss of information on A is (p — rs)/ 

.  r$(p — 1) and that on each of AB, AC, ABC is p[rs (p — 1). 


| > 18. Given m sets (ay, az, ..., a5) where each a; can assume q values 0, 1, 2, 
I .. q — 1. A subset N of these sets is called an orthogonal array of strength 2 

oa if all combinations of every pair of coordinates occur an equal number of 
p=N/¢ times. Show that these orthogonal arrays can be utilised to con- 
struct designs for sı X S2 experiment in blocks of s, or s; plots preserving the 
2 main effects of A or B respectively. 


ariance of the estimates of p objects in N 
1 and in the most efficient design it is g*/N. 
e taken as 


19. Show that the average v: 
weighings is (o?/p). Trace (X' X)~ 
Hence the efficiency of a weighing design can b 


PIN Trace (X' X)-1—(1/N) Harmonie mean of the eigen values of (X' X) 
Obtaií the efficiencies of the following weighing designs : 
1 1 1 1 AA aO 
R tg oun x " 5" 2 
= = (Bee 

| a)X= |1 1 (2)X= ql (8) DECLA 

1 —Il LESI : 
1 —1 —1 —1 (ab gb on 


L 20. If Hy. isa matrix of (N — 1) rows and p columns consisting of +1 
and — 1 such that H’y-1 Hy_1= (N — 11) show that the efficiency of 


the weighing design 
E Ei 
£s [ Hy. ] 


o 


is (N—1) (Np — DI NN +P —9* 


¢ E 
à X= [ "a ] 
EnO 
O x= [ " ] 


21. Show that the most eficient weighing design for N = 2" can be 

written down by writing (1) @ b, “=> ab, ..., abe, along columns and A, B, ie 

$ AB ABC... along rows and filling up the body of the table by 1 or —1 
Pom Eos iss 


Obtain the efficieney when 


TEL 


116 DESIGN AND ANALYSIS OF EXPERIMENTS 


according as the rows have an even or odd number of letters in common with 
the columns and finally adjoining a row of + 15, 


22. If N=2 (p"--1)—4 k, where p is an odd prime or zero, show that an 
NXN matrix A can be constructed with + 1 and —1 such that A^ A=NIy. 


23. We are interested in the total weight of v ob 
balance. Naturally ell the objects can not be weighed simultaneously. We 
take as the weighing design X — N’ where N is the incidence matrix of a 


BIBD with parameters v, b, r, k, à. Show that the variance of the best unbiased 
linear estimate of the total weight is c?v[rk. 


jects by means of a spring 


24. A weighing design has the following X and y matrices : 


= ear esi 


20 
1 1 1 1 80 
X= —1 —1 1 1 y= <10 
—1 1 —1 1 —'5 
—1 1 1 —1 n Dy 


(i) Obtain the best estimates of individual Weights 
(ii) Calculate the efficiency of the design 


(iii) Indicate how you would proceed if you feel that there is bias in 
the balance. 4 


(iv) Can you estimate the bias from the above material? If not, state : 
what more information you need. 


25. A square matrix witha single entr 


y of 1 in each tow and in each column 
and with all other entries zero is called a 


: permutation matrix. Show that the 
incidence matrix of a symmetrical BIBD with i plots per block can be ex- 
pressed as a sum of k permutation matrices, 


26. A certain dance is attended b itis 
2 . oy has 
previously introduced to exactly k gi P en 
introduced to k boys. No one wants to m. 
these assumptions, will it be always possi 


with each other in such a way that no further intr 


such a way that the 
the number of times 
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in Rectangle. Show that anecessary and sufficient 


the integer i occurs in the Lat 
may be extended to an m X n 


condition in order that the Latin Rectangle 
Latin Square is that Nü)&rts—-micLhA: 

28. Suppose that we are given vr elements made up of v varieties of objects 
each repeated r times and that the set is arbitrarily arranged in a two-way 
classification of r rows and v columns. Show that it is always possible to 
rearrange the elements in each column so that each row will contain one and 


only one of each variety. 


min. 


29. A pack of playing cards contain the numbers Ace, King, Queen, Jack, 
13 tricks of 4 cards each 


10, 9, 8, ..., 2 of four suits. Arrange the 52 cards into 

in such a way that (i) any two tricks have one number in common (ii) any 
two numbers occur together in some trick and (iii) any trick contains allfour 
suits, 

t have each six leaves grow- 
f the 81 treatments to the 
e on the same plant and 


sen for an experimen 
Plan an allocation o 
eatments occurs one 
h leaf position. 


80. Thirty-one plants cho 
ing serially along the stem. 
leaves such that each pair of tr 
each freatment occurs once on eac 
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